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1. INTRODUCTION 



Previous studies of fast metal combustion such as exploding wire or under water 
shaped charge have demonstrated that a theoretical description of the process calls for 
constructing a transient model rather than existing quasi-steady models. In some particular 
cases of volatile metal combustion of large droplets (larger than 80 micron ), radiation heat 
exchange between the droplet surface and the surrounding oxide cloud cannot be ignored. 
Special care was taken in Ref. 1 to properly account for the oxide cloud and the metal sur- 
face radiation. Oxide cloud formation, convection and deposition can be described by the 
transport equation of the second phase. The metai is treated as a single phase whether it 
is liquid or vapor. This is accomplished through a judiciuos choice of enthalpy which 
accounts for vapor quality and latent heat of vaporization. Metal-vapor interface location is 
found indirectly by searching for a point where the temperature is equal to the saturated 
vaDor temperature corresponding to the local metal partial pressure, and metal vapor qual- 
ity is zero. 

In Ref. 1, we described the equations for soherical, one dimensional, transient com- 
bustion an including oxide as a second pnase. in oraer to imbed these equations in a 
computer program, we will use the general purpose code PHOENICS , Ref. 2. Since this 
code is not restricted to one spatial dimension, once the feasibility of the 1-D model in this 
phase of the work has been demonstrated, it may be extended to 2-D (axisymmetric using r 
and 0 as spatial variables) in a straightforward manner. 

2. THE GENERAL TRANSPORT EQUATION IN PHOENICS 

The underlying idea of PHOENICS difference scheme is that all transport equations are 
re-cast in a uniform format- the "general transport equation". The formating includes the 
energy conservation equation. The governing equations for our model are conventionally 
written as written as follows: 



Energy equation phase 1: 

piDh,/Dt = V(KtVTi) + Q"i 



Energy equation phase 2: 

p 2 Dh 2 /Dt = V(K&T 2 ) + q " 2 



( 1 ) 



( 2 ) 
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Species mass fraction: 



p.Dc./Dt = V(p,D j0 Vc.) + w". (3) 

Momentum phase 1 : 

piDVi/Dt = V(mVv,) + Mu (4) 

Momentum phase 2 : 

p 2 DV 2 /Dt = ^(nz’tfVz) + M 12 (5) 



The continuity equation for each phase: 





(piV,i%i) = m 2 i 


(6) 


d(p2^2)/ct + 


(p 2 V 2 ^? 2 ) = -m 2 i 


(7) 



As shown by Spalding (Ref. 2), these equations can be presented in th following gen- 

tt 

eralized form, utilizing in addition to the property O, source terms S ^ and exchange coef- 
ficients 



p.DO./Dt =7(r 0 7o.) + s " 0 ( 8 ) 

Where: J 

The variable <D. = hi,h 2 ,c.,Vi,V 2 (8a) 

The source s" = Q\ t Q W 2 ,w\,Mi 2 ,M 2 i (8b) 

The exchange coefficient: 

r (p = K/c p .PiD j0 ,fM.H 2 (8c) 
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This generalization enables the use of a single numerical scheme for all transport 
equations. PHOENICS is built in three main parts described in Fig. 1. 





Fig. 1 The various subprograms in PHOENICS. 



EARTH 

SATLIT 



FLDDAT 

GROUND 



Include the numeric finite differences scheme and written as 
binary code. 

The management program, where the problem is declared: 
steady/unsteady, parabolic or eliptic , boundary conditions, 
relaxation factors etc. 

Initial values definition for each variable O^being used. 

For variables updating and sources definition. 
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3. CORRECTION CONDUCTION TERM IN ENERGY EQUATION 



Special care should be excercised in recasting the energy equation into the general 
format. The problem arises from the heat conduction term which should be a divergence of 
K T. Since the generalizated energy variable is enthalpy rather than temperature, the heat 
flux is related indirectly to the enthalpy. In PHOENICS this difficulty is resolved by assum- 
ing a constant ( temperature and composition independent) value of c . 

This approximation is inacceptable for our model. In the sequel, we dwell on the rea- 
sons for that and on the proposed remedy. 

There are three major factors that prohibit the assumption of constant c . These are: 

(a) Large temperature gradients. 

(b) Large composition gradients (especially in the combustion 
region and around the liquid metal/vapor interface). 

(c) The conductive heat flux constitutes a non-negligible fraction 
of the total heat flux (conduction + convection + radiation). 

The remedy for this problem is to introduce a suitable source term into the PHOENICS- 
generalized energy equation. This source term is just the difference between the correction 
heat flux divergence and the enthalpy flux divergence. For our model, the appropriate 
energy source correction term is: 



S corr = <3[(K/c p )[c p (dc^dy) + c p (dcjdy) + c p {dcjdy)} (T-298)]/<?y 

- <3[(K/c p )h° wv ((3c wv /^y)]/ <?y 

- ^[(K/c p )[(1-x)h Mfg (<? CM /<?y) - h Mfg c M (dx/<?y)]]/(?y (9) 



Note that last term in Eq. 9 vanishes except in the transition region between metal liquid 
and vapor, where x is a fraction number. 

In the case of a 1-D configuration, finite- difference gradient approximation is available 
in PHOENICS and described in Ref. 3. 
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c[(D^h/ay) = (h E -h p ) / [ 5y # _/r p + 5y ef /r E ] 

- ( h P- h H ) > 1 5 v M - /r P + 



( 10 ) 



Where p is the center node and E and W are its sides neighbors, as described below: 



5 VW 
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7 w- 
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p 



5 y„ 



ov 




e 



We used this difference scheme in computing the corrective energy terms. 
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4. VARIABLE NAMES AND FLOWCHART 



In this chapter we wil present the metal combustion variables used in PHOENICS and 
the implementation of the model equations as they appear in chapter 6 of Ref. 1.. The defi- 
nition of the problem and its boundary condition was declared in SATLIT. The initial condi- 
tions were set in FLDDAT. The various sources, the iterations procedures of surface and 
oxide cloud location together with the other combustion model were inserted in GROUND. 
The reader can follow the flow chart together with using the FORTRAN list in the appendix. 



NAME IN 


NAME OF VARIABLE 


APPEAR IN 


* 


THE PROGRAM 


IN EQUATION 


EQUATION NO. 




GAI 


s. 

i 


6-27 


C 


GAO 


s 

o 


6-23 


C 


GAP 


a 

p 


6-26,30 




GAQ 


E/R 


6-34 




GCDT 


rt 

w 

m 


6-34 


C 


GCP 


C 

P 3 


6-51 c 


C 


GCPI 


C 

p g 


6-26C 




GCPOX 


c 

P ox 






GCPO 


c 0 

p 


6-51 c 




GCPI 


c 1 

p 


6-51 c 




GCP2 


c 2 

p 


6-5 1c 




GCP3 


c 3 

p 


6-5 1c 




GCP4 


c 4 

p 


6-51 c 




GC1 


C M 


6-39 


C 


GC2 


C wv 


6-39 


C 


GC3 


C H 


6-39 


C 
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GC1MIN 


N . 

nun 


parameter for 
radiation beginning 


GD 


D.o 


6-49 C 


GDD 


Dpm jk 


6-48 C 


GDI 


Pi 


6-55 C 


GD2 


P2 




GEG 


e 

g 


6-26,30 


GEP 


z 

p 


6-25 C 


GEW 


£ 

w 


6-26,30 


GF 


cjc 

n wv 


6-35a,b 


GFCW 


F 

C-W 


6-29 


GFG 


*H 


6-37a 


GFOX 


j ox 


6-37 b 


GGAMH1 


k/c 

p g 




GHFGM 




6-4 


GHFOX 


h ^ 
oxfg 


6-32 


GHFU 


H , 

fu 


6-5 


GHOXS 


tt 


6-6a C 


GHOXV 


Q ,3 


6-6b C 


GHOXO 


0 

h 

ox 


6-33a 


GH1M 


h Ml 


6-9 C 


GH3M 


h M3 


6-10 C 


GK 


k. 

3 


6-51 a C 


GKA 


B 


6-34 


GKMF 


K Mf 
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GKO 


k io 


6-5 la 


GK1 


k 3i 


6-5 la 


GK2 


k j2 


6-5 la 


GK3 


k 33 


6-5 la 


GK4 


k * 


6-5 la 


GL 


1 


6-25 


GLSOX 


^oxsf 


6-6b 


GM 


**3 


6-5 la 


GMDT 


t* 

w 

ox 


6-37a,b 


GMM 


M./M k 


6-52 


GMO 


^30 


6-51, b 


GM1 


•*31 


n 


GM2 


^32 


H 


GM3 


^33 


ft 


GM4 


^3* 


tf 


GN 


N. 

1 


6-23 


GNMI 


N i 










GNN 


N 


6-24 


GPAL 


Ph 


6-18 


GPT 


P 




GRC 


r 

C 


6-29 


GRCI 


r . 

Cl 


6-27,28 


GRCo 


r co 


// 


GRFC 


q ox + 
not 


6-26 


GROF 


P Mf 


6-21 
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GROG 


^Mg 


6-21 


GROO 


Pox 


6-23 


GROU 


P1U1 


6-1,39,45 


GROUCM 


piUidc^r 


6-39 


GROUOX 


piUi^c /dr 

r wv 


6-39 


GRW 


^wnet 


6-30 C 


GR2 




6-42.46 


GSIG 


<7 


6-26,30 


GSUMD 




6-49 


GSUMF 


***** 


6-53,54 


GSUMK 


K 


6-54 C 


GSUMN 


EN.AV. 


6-24 


GSUMU 




6-53 C 


GSUMV 


EAV. 


6-24 


GSXM2 




term in Eq. 9 


GSXM1 




" 


GTEMP 


T. 


6-17 C 


GTEMPO 


Tz 


6-33a 


GTG 


T 

g 


at steam 


GTOX 


T 

ox 


at cloud 


GTW 


T 

w 


at x = 0 


GVISC 


tl 




GXMET 


X 


6-19 C 


GXMETA 


X' 


assumed for iteration 



’ Calculated in that equation. 
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Fig. 2 General flow-chart of the model in PHOENICS 
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Fig. 3 Chapter 5 in GROUND for the various sources in 
the model 
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* Gas Phase * 




Fig. 4 Metal vapor distribution and surface location 
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Fig. 5 Oxide cloud location and cloud emissivity 
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5. SUMMARY 



We presented a prescription for embedding the 1-D droplet combustion model that had 
been proposed in a Ref. 1 into the PHOENICS code. Extensive programming modifications 
are required, primarily due to the introduction of special constructed (finite- difference) 
source terms. These terms are of two distinct types. The corrective type was inserted into 
the energy equation in order to offset errors introduced by the general equation format of 
PHOENICS, which is based on an assumption of constant c , that does not hold in the 

p 

present model. 

The second type is true mass and energy source terms. The mass source is due to the 
progress of chemical reactions. The energy source is due to radiative heat transfer. 

The debugging phase has just about been completed at the time this report was writ- 
ten. We regard the 1-D model as a test case for droplet combustion in steam. Most- but not 
all - the physical factors are represented in the model. Specially, oxide deposition on a 
moving droplet can be modeled in a 2-D (axisymmetric) computation. 

Once the 1-D test computations prove statisfactory, the road is clear to proceed with the 
2-D model. 
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In this appendix we will pesent the FORTRAN listing of PHOENICS which includes our 
model for metal combustion in steam. 



CSDIRECTIVEXXMAIN 

C XFILE NAME: MODGRD.FTN 

C XINCLUDED SUBROUTINES: THE MODELS OF MAIN, GROUND & STRIDE. 

C XDOCUMENTATION: PHOENICS INSTRUCTION MANUAL (SPRING 1983). 

C XSATEL LITE FILE NAME: MODSTL.FTN 

COMMON/ 1 SH I FT/ 1 1 1 ( 57 ) , NFMAX 

C SET F-ARRAY DIMENSION AS NEEDED, & SET NFMAX ACCORDINGLY. 

COMMON F(25000) 

NFMAX=25000 
CALL MAIN1 
STOP 
END 

C$DIRECTIVE*XGRCUND 

SUBROUTINE GROUNDCIRN, ICHAP, ISTP, ISWP, I2ED, INDVAR ) 

CSINCLUDE 9 , CMNGUSSI . FTN/G 
LOGICAL LOGIC1, LOGIC 
DIMENSION LOGIC(IOO) 

COMMON/ L DATA/ LOGIC1 (309) 

EQUIVALENCE (LOGIC(l), LOGIC1 (210)) 

DIMENSION INTGR ( 1 0 0 ) 

COMMON/I DATA/INTGRI ( 199 ) 

EQUIVALENCE ( INTGR( 1 ) , INTGR1 ( 95) ) 

DIMENSION RE(100) 

COMMON/RDAT A/REK421 ) 

EQUIVALENCE ( RE( 1 ) , RE1 ( 322 ) ) 

COMMON/BOUND/LOCREG(60) , 

8TR1,CP1RK7),VP1R1(7),CP2RK5),VP2RK5),CPNRK5),VPNRK5), 
8TR2,CP1R2(7),VP1R2(7 ),CP2R2(5) , VP2R2(5) , CPNR2( 5 ) , VPNR2( 5 ) , 

8TR3 »CP1R3(7),VP1R3(7)»CP2R5(5),VP2R3(5),CPNR3(5),VPNR3(5), 
aTR4,CPlR4(7),VPlR4(7 ) , CP2R4( 5 ) , VP2R4( 5 ) , CPNR4( 5 ) ,VPNR4(5), 
8TR5,CP1R5(7),VP1R5(7),CP2R5(5),VP2R5(5),CPNR5(5),VPNR5(5), 
8TR6,CP1R6(7 ),VP1R6(7 ),CP2R6(5) , VP2R6 ( 5 ) , CPNR6 ( 5 ) , VPNR6 ( 5 ) , 
&TR7,CP1R7(7 ), VPIR7(7),CP2R7(5) , VP2R7 (5 ) , CPNR7 ( 5 ) , VPNR7(5), 
8TR8,CP1R8(7),VP1R8(7),CP2R8(5),VP2R8(5),CPNR8(5),VPNR8(5), 
8TR9,CP1R9(7),VP1R9(7),CP2R9(5),VP2R9(5),CPNR9(5),VPNR9(5), 
8TR10,CP1R10(7),VP1R10(7),CP2R10(5),VP2R10(5),CPNR10(5),VPNR10(5) 
DIMENSION TCVREG(350) 

EQUIVALENCE ( TCVREG( 1 ) , TR1 ) 

CSINCLUDE 9,GUSSEQUI . FTN/G 

LOGICAL CARTES, POLAR, SPDATA, SKEW, TWODYZ, ONEDZ, STOVAR( 25), 

& S0LVAR(25), PRINT (25), RESID(25),CMPRSS, CON EMU, LSP1, 

S CONRHO, EMDOT,ONEPHS,INCORE(10),ISAVED,SAVEI,SAVEM, 

5 RESTRT , XCYCL E, MONITR, REGION( 1 0 ) , STEADY, WHOL EP , SL ABPP , 

6 RAIN, BLOCKZ, PWST AG, RUN( 30 ) , PLOT , RESMAP, FLAG, BLOCK, 

& TEST, C0NCK4), DISTIL 

LOGICAL SPRESS , PARAB, DON ACC, OVERLY, SACC, 

& GUSSIE, CATHY, CONNIE, CORA, ESTER, FLASH, FLORA, 

8 FOCS,GENMIX, HESTER, PICALO, PL ANT, SPLASH, HELGA, 

a TACT ,TIBALT,TOPSI, PAMELA, TABLES, WSTAG, 
a CONMOD,GROSTA,SUBPST , SUBWGR 
EQUIVALENCE ( CARTES , L OGIC1 ( 1 )),( POLAR , L OGIC1 ( 2 )) , 
a ( SPDATA, LOGICK 3)), (SKEW, LOGIC1 (4)), (TWODYZ, L0GICK5)), 
a (ONEDZ, LOGICK 6 )), ( STOVAR( 1 ) , LOGIC1 (7 ) ) , 

a ( SOL VAR( 1 ), LOGIC1 ( 32 ) ) , (PRINT( 1 ) , LOGICK 57)), 
a (RESID(1),L0GIC1(82)), (CMPRSS, LOGICK107)), 
a ( CON EMU ,LOGIC1(108)),(LSP1,LOGIC1(109)), 

a (CONRHO, LOGIC1(110)),( EM DOT, LOGICl(lll)), 

a ( ONEPHS, LOGIC1 ( 112) ) , ( I NCORE( 1 ) , LOGICK 113)) 

EQUI VALENCEdSAVED, LOGICK 123) ), (SAVEI, LOGICK 124)) , 
a (SAVEM, LOGICK 125)), (RESTRT, LOGICK 126)), 

a ( XCYCL E, LOGICK 127 ) ), (MONITR, LOGICK 128) ), 

a ( REGI ON( 1 ), LOGICK 129) ),( STEADY, LOGICK 139)), 

a (WHOLEP, LOGICK 140) ),( SL ABPP, LOGICK 141 ) ), 

a ( RAIN, LOGICK 142) ), (BLOCKZ, LOGIC1 (143)), 
a ( PWSTAG, LOGICK 144) ),(RUN(1), LOGICK 145)), 

a ( PLOT, LOGICK 17 5 )), (RESMAP, LOGICK 176)) , 

a (FLAG, LOGICK 177)), (BLOCK, LOGICK 178)), (TEST, LOGICK 179)) 

EQUI VAL ENCE(C0NC1( 1 ), LOGICK 206 ) ) 

EQUIVALENCE (SPRESS, LOGIC( 100 )),( PARAB , L0GIC(99 )) , 
a (DON ACC, L0GIC(98) ) 

EQUIVALENCE ( OVERLY, LOGIC( 96 )),( SACC, LOGIC( 95) ) , 
a (WSTAG, L0GIC(94) ) 



MEL00010 
MEL00020 
MEL00030 
MEL00040 
MEL00050 
MEL00060 
MEL 0 0 070 
MEL00080 
MEL00090 
MEL00100 
MEL 00110 
MEL 00120 
MEL 0 0130 
MEL00140 
MEL00150 
MEL 0016 0 
MEL 0017 0 
MEL00180 
MEL00190 
MEL 00200 
MEL 00210 
MEL 00220 
MEL 00230 
MEL 00240 
MEL 00250 
MEL00260 
MEL 0 0270 
MEL 3023 0 
MEL 00290 
MEL 00 300 
MEL 0031 0 
MEL 00 320 
MEL 00330 
MEL00340 
MEL00350 
MEL00360 
MEL0037 0 
MEL00380 
MEL00390 
MEL00400 
MEL00410 
MEL00420 
MEL00430 
MEL00440 
MEL00450 
MEL 0046 0 
MEL00470 
MEL00480 
MEL00490 
MEL00500 
MEL00510 
MEL 0 0 520 
MEL 00530 
MEL00540 
MEL 00 550 
MEL00560 
MEL00570 
MEL00580 
MEL00590 
MEL 0 06 0 0 
MEL00610 
MEL 00620 
MEL00630 
MEL00640 
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EQUIVALENCE ( DISTI L , L OGIC( 90 ) ) 

EQUIVALENCE (PAMELA, LOGIC( 67 )), (TOPSI, LOGIC( 66) ) , 

8 (TIBALT,L0GIC(65)),(TACT,L0GIC(64)),(HELGA,L0GIC(63)), 

& (SPLASH, L0GICC62) ) , (PLANT, L0GIC(61) ) , ( PICALO, L0GICC6 0 ) ) , 
a ( HESTER, LOGIC( 59) ) , (GENMIX, LOGIC( 58 ) ) , ( FOCS , LOGIC( 57 ) ) , 
a ( FLORA , LOGIC( 56 )),( FLASH, LOGIC ( 55) ), ( ESTER, L OGIC( 54) ) , 

a ( CORA ,L0GIC(53)), (CONNIE, L0GICC52)), (CATHY >LOGIC(51)), 

a (GUSSIE,LOGIC(50)) 

EQUIVALENCE (TABLES , LOGIC( 48 )),( CONMOD, LOGIC( 47 )) , 
a ( GROSTA, LOGIC! 46 )), (SUBPST, LOGIC! 45)) 

EQUIVALENCE ( SUBWGR, LOGIC( 43 ) ) 

INTEGER FSTEP,FSWEEP,TSTSWP,ITAB(8),MTABVR(8) 

DIMENSION ISPCSO( 25 ) ,LITER(25) 

EQUIVALENCE ( NX, INTGR1 ( 1 ) ) , ( NY, INTGR1 ( 2 ) ) , 
a (NZ,INTGR1(3)),(ISPCS0(1),INTGR1(4)), 
a (NREGN,INTGR1(29) ), (NPHI, INTGRK30) ), 
a ( LITKE, INTGR1 ( 31 ) ) , ( L ITHYD, INTGR1 ( 32 ) ) , 

a (LITH,INTGR1(33)), ( LITCNC, INTGR1 ( 34 ) ) , 
a (LITSLB, INTGRK35)), ( NRUN, INTGR1 ( 36 ) ) , 

a (LITER(l), INTGRK37)) , ( FSTEP, I NTGR1 ( 62 ) ), 

a ( FSWEEP, INTGR1 (63) ) , ( LSTEP, INTGR1 ( 64 ) ) 

EQUIVALENCE (LSWEEP, INTGRK65) ) , (NPRINT, INTGRK66) ) , 
a (IERRP,INTGR1(67)),(IMAXP,INTGR1(68)), 
a ( I EMU1 , INTGR1(69)),(IXMON,INTGR1(70)), 
a ( IYMON, INTGRM71)), (IZM0N,INTGR1(72)), 
a ( KEMU, INTGR1(73)),( KMAIN, INTGR1 (74 ) ) , 
a (KINDEX, INTGRK75) ), (KGEOM, INTGRK76) ), 
a (KINPUT,INTGR1(77) ) , ( KSODAT , INTGR1 ( 78 ) ) 

EQUIVAL ENCE( KFLAG, INTGR1 ( 79 ) ) , (KCOMPF, INTGR1 (80 ) ) , 

8 ( KSORCE, INTGRK81) ), (KLES1D, INTGR1 (82) ), 

a (KLES2D,INTGR1(83) ), (KLES3D, INTGR1 ( 84) ) , 
a (KCOMPP, INTGR1 (85) ) , (KADJST, INTGRK36 ) ) , 
a ( KFLUX, INTGRK87) ) , ( KSHI FT , INTGR1 ( 88 ) ) , 
a (KOUTPT, INTGRK89) ), (KDIF, INTGR1 ( 90 ) ) , 
a ( KCOMPU, INTGR1 ( 91 ) ) , ( KCOMPV, INTGR1 ( 92 ) ) , 
a (KCOMPW, INTGR1 (93) ) , (KCOMPR, INTGR1 (94) ) 

EQUIVALENCE ( KELIN, INTGR( 100) ) , (MEANDF, INTGR( 99 ) ) , 
a ( NUMCLS, INTGR(98) ) , 

a (IRH01,INTGR(97)),(IRH02,INTGR(96)), 
a (IZW1,INTGR(95)),( IZW2, INTGR(94)), 
a ( MGRID, INTGR(93)),( KWAL L , INTGR( 92) ) , 
a ( IZPR1 , INTGR( 91 ) ) , ( IZPR2, INTGR( 90 ) ) , 

a ( ISTPR1 ,INTGR(89)), ( ISTPR2, INTGR(88)) , 
a ( NTPRIN, INTGR( 87 ) ) 

EQUIVAL ENCE ( IMDOT , INTGR(86 ) ) , ( IHSAT , INTGR( 85) ) , 
a (ICFIP,INTGR(84)),(NTABLE,INTGR(83)),( NTABVR, INTGR( 82) ) , 

a (LINTAB,INTGR(81)),(NPRTAB,INTGR(80)),( NMON, INTGR(79)), 
a (ITAB(l), INTGR (71)), ( MTABVR(l), I NTGR (63)), 

a (NZPRIN, INTGR(62) ) , 

a ( NXPRIN, INTGR (61)), (NYPRIN, INTGR (60)), 

a ( LDISTL , INTGR(59 ) ) , (TSTSWP, INTGR(58) ) , 

a ( KDBEXP, INTGR( 57 ) ) , ( KDBRHO , INTGR( 56 ) ) , 
a ( KDBMDT , INTGR( 55 ) ) 

EQUIVALENCE ( I LOOP1 , I NTGR (53)), ( I LOOPN, INTGR (52)), 
a (IPBP,INTGR(51)),(IZPRF,INTGR(50)), 
a (IZPRL,INTGR(49)),(ISTPRF,INTGR(48)), 
a ( ISTPRL , INTGR( 47 ) ) , (KDBGEN, I NTGR (46 ) ) , 
a (IVELF, INTGR (45)), (IVELL,INTGR(44)), 

a (IKEF,INTGR(43)),(IKEL,INTGR(42)), 
a ( IENTF, INTGR (41)), ( IENTL, INTGR ( 40) ) , 

a (ICNCF,INTGR(39)),(ICNCL, INTGR (38)), 
a (NVEL,INTGR(37)),(NKE,INTGR(36) ), 
a ( NENT , INTGR (35)), ( NCNC, INTGR ( 34 )) , 

a ( NZSTP , INTGR( 33 ) ) , ( NPRMNT , INTGR( 32) ) 

DIMENSION SPARE1(20),XUDIST(30),YVDIST(30), ZWDIST ( 50 ) , 
a SIGMA(25),CRIT(25) , DTFALS( 25) , RESREF( 25) , 
a TITLE(25),FIINIT(25) ,TFRAC( 30 ) 

EQUIVALENCE ( SPARE1 ( 1 ) , RE1 ( 1 ) ) , 
a (TFRAC(1),RE1(21)), (XULAST, RE1 ( 51 ) ) , 
a (YVLAST,RE1(52)),(ZWLAST,RE1(53)), 
a (XUDIST(l) ,RE1(54)), ( YVDIST( 1 ) , RE1 ( 84 ) ) , 

8 ( ZWDIST( 1 ), RE1 ( 1 14 )), (SIGMA (1),RE1( 164) ) , 



MEL 00730 
MEL00740 
MEL 00750 
MEL00760 
MEL 0077 0 
MEL00780 
MEL00790 
MEL00800 
MEL00810 
MEL00820 
MEL 00830 
MEL 00840 
MEL 00850 
MEL 0086 0 
MEL0087 0 
MEL 00880 
MEL 00890 
MEL00900 
MEL 009 1 0 
MEL 00920 
MEL00930 
MEL 00940 
MEL 00950 
MEL 0096 0 
MEL 0097 0 
MEL 00980 
MEL 00990 
MEL 01000 
MEL01010 
MEL 01 020 
MEL 01030 
MEL01040 
MEL 01050 
MEL 01 06 0 
MEL 01 07 0 
MEL01080 
MEL 01090 
MEL Oil 00 
MEL01110 
MEL01120 
MEL 01130 
MEL01140 
MEL 01150 
MEL01160 
MEL0117 0 
MEL01180 
MEL01190 
MEL 01200 
MEL 0121 0 
MEL01220 
MEL 01230 
MEL01240 
MEL 01250 
MEL 0126 0 
MEL 0127 0 
MEL 01280 
MEL 01290 
MEL01300 
MEL01310 
MEL 01320 
MEL 01330 
MEL01340 
MEL01350 
MEL 0136 0 
MEL 0137 0 
MEL01380 
MEL 01390 
MEL 01400 
MEL 0141 0 
MEL01420 
MEL 01430 
MEL 01440 
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& (CRIT(1),RE1(189)), (DTFALS(1),RE1(219)), 

& (RESREF(1),RE1(239)), ( EMU1 , RE1 ( 269 ) ) , 

a (RH01,RE1(265)), C RH02 , RE1 ( 266 ) ) , 

8 (TLAST,RE1(267)),(CFIPS,RE1(268)) 

EQUIVALENCE ( AMDOT , RE1 C 269 ) ) , ( FI INITU ) , RE1 ( 27 0 ) ) , 
a (RELAXP,RE1(295)), (TITL EC 1 ) , RE1 ( 296 ) ) , 
a (DT,RE1C321)), 

a (RINNER,ZWDIST(50)) , ( SNAL FA, ZWDISTC 99 ) ) , 

a ( PBAR , ZWDIST ( 98 ) ) 

EQUIVALENCE ( SLOEMU , REC 10 0 ) ) , ( SLORHO, RE( 99 ) ) , ( RLXRHO, RE( 98 ) ) . 
a ( RHOMAX, RE( 97 ) ) , C RHOMIN, RE( 96 ) ), 

a (EMUMAX,RE(95)),(EMUMIN,RE(99)), 
a (TKEMAX,RE(93)),(TKEMIN,RE(92)), 
a (RLXPZ, RE( 91 ) ) , 

a (ABUOY,REC90)),(HREF,REC89)), 

8 (AGRAVX, REC88 ) ) , ( AGRAVY , RE ( 37 ) ) , 
a (AGRAVZ,RE(36)),(ATIME,RE(85)) 

EQUIVALENCE (RLXPXY, REC89) ) , 
a ( EPSMAX, REC 83 ) ) , (EPSMIN,RE(32) ) . 
a ( ARHOl , REC 81 ) ), ( BRHOl, REC30) ) , ( CRHOl , RE(79 ) ) 

EQUIVALENCE ( AZW2, REC 73 ) ) , C BZW2, RE( 79 ) ) , 
a (PINT,REC75)),(CZW2,REC76)), 
a (DZ,RE(78) ) , (ZWL,RE(77 ) ) , ( ZW2M1T, RE(72 ) ) , 
a (VELMIN, REC70)) , ( VELMAX, REC 6 9 ) ), 
a ( FALCOM, REC 68 ) ) , (CMDOT,RE(67 ) ), 

a (RLXMDT,RE(66)), ( AMDTMX, REC 65 ) ) , ( AMDTMN , RE(64 ) ) , 
a ( RADMAX, RE(63) ) , (RADMIN,RE(62) ), 
a (PRESS, RE(61) ), CCP1,RE(60) ), 
a ( CP2 , REC 59 ) ) , (CP3,RE(58) ) , ( FSTOIC , REC 57 ) ) , 

a C ARRCON , REC 56 ) ) , ( PREEXP ,RE(55)),(CA1I,RE(5 ( 4)) 

EQUIVALENCE (CA2I,RE(53) ) , C CA1 3 , REC 52 ) ) , ( CA23 , REC 51 ) ) , 
a CH1SAT,RE(50)),(H2SAT,RE(‘49)), 
a ( EMULAM, REC 48 ) ) 

EQUIVALENCE ( TWPRCN , REC 92 ) ) , ( PSATEX , REC 91 ) ) 

EQUIVALENCES EXP, REC 39) ), C EWAL L , REC 38 ) ) , 
a (SIGE,RE(37)),(SIGK,RE(36)),(TAUDK,RE(35) ) 

DIMENSION L0C1(6),L0C2(6),L0C3(6),L0C4(6),L0C5(6), 
a L0C6C6), LOC7(6),LOC8(6),LOC9(6),LOC10(6) 

EQUIVALENCE ( LOCREGC 1 ) , LOCI ( 1 ) ) , ( L0CREGC7 ) , L0C2C 1 ) ) 

(L0CREG(13),L0C3(1)),(L0CREG(19),L0C‘4(1)) 

(LOCREGC 25 ),L0C5(1)), (LOCREGC 31), L0C6C1)) 

(LOCREGC 37 ), LOC7 ( 1 )),( LOCREGC 93) , L0C8 ( 1 ) ) 

(LOCREGC 99), L0C9C1)), ( LOCREGC 55 ), L0C10 ( 1 ) ) 

XUFRACC 30 ) , YVFRACC 30),ZWFRAC(30), ENDITC25) 
XFRAC(30),YFRAC(30),ZFRAC(30) 



EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
DIMENSION 
DIMENSION 



EQUIVAL ENCECXUDISTC 1 ) , XUFRACC 1 ) , XFRACC 1 ) ) , 
a (YVDIST(l) , YVFRACC 1 ) , YFRACC 1 ) ) , 

a ( ZWDIST ( 1 ) , ZWFRACC 1 ) , ZFRACC 1 ) ) , ( L ITCNC, LITC) , 

a (ENDIT(1),CRIT(1)),(RLXP,RELAXP) 

EQUIVALENCE ( NPRMON, LITER CD) 

CLUDE 9 , NMLIST . FTN/G 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 1 STARTS 



++MEANING OF SUBROUTINE ARGUMENTS-. 

IRN=RUN NUMBER; ICHAP=CHAPTER CALLED; ISTP=TIME STEP; 
ISWP=SOLUTION SWEEP; IZED=Z-SLAB; INDVAR: SEE CHAPTERS BELOW. 
++USER- INTRODUCED VARIABLES a ARRAYS-. 

TO AVOID CONFLICT WITH VARIABLE NAMES USED IN COMMON, ALL 
VARIABLES INTRODUCED BY THE USER SHOULD HAVE NAMES STARTING 
WITH 'G* IF REAL, 'J' IF INTEGER, AND ’G* OR 'J 1 IF LOGICAL. 
THUS GDZCIZ) MIGHT BE A Z-INTERVAL ARRAY; 

GW1 ( IY, IX) A 2-D ARRAY FOR AXIAL VELOCITY; ETC. 

USER-GENERATED SUBROUTINES SHOULD BE NAMED CORRESPONDINGLY, EG 
SUBROUTINE GVISCC GTEMP ,GCNC, GVSC) , FOR COMPUTING VISCOSITY 
FROM CONCENTRATION a TEMPERATURE. 

++GROUND-TO- EARTH CONNECTING SUBROUTINES: 

XUSE GETC NAME, G ARRAY , NY, NX) TO PUT VALUES OF VARIABLE NAMED 
’NAME' INTO ARRAY 'GARRAY' DIMENSIONED GARRAYCNY, NX) . 

XUSE SET(NAME,IXF,IXL,IYF,IYL, GARRAY, NY, NX) TO SET VARIABLE 
’NAME’ TO GARRAYC IY, IX) OVER THE REGION: IXF-IXL a IYF-IYL . 
XUSE PRNSLBC NAME) TO PRINT VARIABLE 'NAME' OVER X-Y PLANE. 

XUSE ADDC NAME, IXF, IXL , IYF, IYL , TYPE, CM, VM, CVAR, WAR, NY, NX) 



MEL 01950 
MEL 0196 0 
MEL 0197 0 
MEL01980 
MEL 01990 
MEL01500 
MEL01510 
MEL01520 
MEL 01530 
MEL01590 
MEL 01550 
MEL 01560 
MEL0157 0 
MEL01580 
MEL01590 
MEL01600 
MEL 01610 
MEL 01620 
MEL 016 30 
MEL 01690 
MEL 01650 
MEL01660 
MEL 0167 0 
MEL01680 
MEL 016 90 
MEL 017 00 
MEL01710 
MEL 01720 
MEL01730 
MEL01790 
MEL 017 50 
MEL01760 
MEL 0177 0 
MEL 01780 
MEL01790 
MEL01800 
MEL 01310 
MEL 01820 
MEL 01830 
MEL01890 
MEL 01850 
MEL 0186 0 
MEL 0 187 0 
MEL01880 
MEL01890 
MEL01900 
MEL01910 
MEL01920 
MEL01930 
MEL01990 
MEL 01950 
MEL 01960 
MEL 0197 0 
MEL01980 
MEL01990 
MEL02000 
MEL02010 
MEL 02020 
MEL 02030 
MEL 02090 
MEL 02050 
MEL 0206 0 
MEL 0207 0 
MEL02080 
MEL02090 
MEL 021 00 
MEL02110 
MEL 02120 
MEL 02130 
MEL02190 
MEL 02150 
MEL02160 



17 



FILE: MEL3GRD FORTRAN A1 



C TO ADD SOURCE TO VARIABLE NAMED 'NAME' (SEE CHAPTER 5). 

C XUSE READIZ(IZED) IN CHAPTERS 1, 2, 8, & 9 TO ACCESS PI,... DM 

C & VOL , . . . AHDZ . (SEE FOOTNOTE TO LEGALITY TABLE) 

C XUSE GET1 D( NAME, GARRAY, NDIM) TO PUT VARIABLE NAMED 'NAME' IN 
C ONE-D ARRAY 'GARRAY' DIMENSIONED NDIM, THUS: 

C CALL GET1 D( NAME, GNX, NX) FOR XG,...DXG & DIMENSION GNX(NX); 

C CALL GET1D( NAME, GNY , NY) FOR YG,...RV 8 DIMENSION GNY(NY) ; 

C CALL GET1D( NAME, GNZ, NZ) FOR ZG,... WGRID 8 DIMENSION GNZ(NZ). 

C+++++L EGAL ITY TABLE FOR USE OF EARTH-CONNECTING SUBROUTINES: 

C ENTRIES IN TABLE GIVE CHAPTERS IN WHICH SUBROUTINES CAN BE 

C USED FOR VARIABLES IN LEFT-HAND COLUMN. (SUBROUTINE 

C STRIDE IS REGARDED AS BEING IN CHAPTER 3) 



c 

c 

c 


: VARIABLE: i 


GET & : 

PRNSLB : 


SET 


i ADD 


: READIZ i 


: GET1D : 


L 

c 


: P 1 - RZ 




ALL 


6 8 7 


5 


1,2, 3, 9 


NONE : 


c 


: P 1 0 - RZH 




3-7, 10-16 


3 


NONE 


NONE 


NONE : 


c 


: VO L -AHDZ 




ALL 


3 


NONE 


1,2, 3, 9 


NONE : 


c 


: D1DP 




NONE 


10 


NONE 


NONE 


NONE : 


c 


:D2DP 




NONE 


11 


NONE 


NONE 


NONE : 


c 


: Mill , MUl H 




5,13-16 


12 


NONE 


NONE 


NONE : 


c 


:EXCO(L,H) 




NONE 


13 


NONE 


NONE 


NONE : 


c 


: C FP 




5 


14 


NONE 


NONE 


NONE : 


c 


: MDT 




5 


15 


NONE 


NONE 


> NONE : 


c 


: HST1 , HST2 




5 8 15 


16 


NONE 


NONE 


: NONE : 


c 

n 


:XG -WGRID 




NONE 


NONE 


NONE 


NONE 


ALL : 



C NOTES ON ABOVE TABLE: 

C XIN CHAPTERS 1, 2, 8, & 9 VARIABLES PI... DM 8 GEOMETRY 
C VOL... AHDZ CAN BE ACCESSED BUT ONLY IN CONJUNCTION WITH 
C USE OF READIZ, THUS: 

C DO 1 IZED=1 , NZ 

C CALL READIZ(IZED) 

C 1 CALL GET ( . . . AS REQUIRED..) 

C XGEOMETRY ACCESSED BY READIZ IS THAT AT INITIAL TIME. 

C XD1DP 8 D2DP ONLY ACCESSIBLE IN UNSTEADY FLOWS. 

C+++++GROUND SERVICE SUBROUTINES: 

C XUSE CONTUR(NAME,IPLANE,ILOC,NINT, II, 12, J1,J2, GARRAY, NDIM) FOR 
C LINE-PRINTER PLOTS OF CONTOURS. 'NAME* = U1,...C4; 

C ' IPLANE' = XPLANE, YPLANE, OR ZPLANE; ILOC SETS IX, IY, OR 

C IZ LOCATION OF IPLANE; II, 12, Jl, 8 J2 SET FIRST 8 LAST 

C CELLS IN HORIZ. 8 VERT. ON PLOT; GARRAY IS 1-D WORKING ARRAY 

C OF DIMENSION NXxNY , NXXNZ, OR NYXNZ DICTATED BY IPLANE; 8 
C NDIM SETS VALUE OF DIMENSION OF GARRAY. 

C XUSE FLD2DA(TITLE, GARRAY, NY, NX) TO PRINT ANY ARRAY DIMENSIONED 
C GARRAY(NY,NX) ; SET 'TITLE' TO REQUIRED NAME ( A HOLLERITH 

C CHARACTERS ONLY) . 

C XUSE FLD3DA(TITLE, GARRAY, NX, NY, NZ, IPLANE, ILOC) TO PRINT ANY 
C ARRAY DIMENSIONED GARRAY( NX, NY, NZ) IN PLANE SPECIFIED BY 

C 'IPLANE* 8 'ILOC' AS FOR CONTUR ABOVE; SET 'TITLE' AS FOR 

C FLD2DA . 

C VARIABLE NAMES FOR USE IN GROUND: 

COMMON/TYPE/CELL, EAST, WEST, NORTH, SOUTH, HIGH, LOW, VOLUME, WALL 
C0MM0N/VAR/P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS, 
8KE,EP,H1,H2,H3,C1,C2,C3,C4,RX,RY,RZ,S1,S2 
C0MM0N/VAR0LD/P10,PP0,U10,U20,V10,V20,W10,W20,R10,R20,RS0, 
8KE0,EP0,H10,H20,H30,C10,C20,C30,C40,RX0,RY0, RZO , SIO , S20 
COMMON/VARLOW/P1 L ,PPL,U1L,U2L,V1L,V2L,W1L,W2L,R1L,R2L,RSL, 
8KEL,EPL,H1L,H2L,H3L,C1L,C2L,C3L,C4L,RXL,RYL,RZL,S1L,S2L 
COMMON/VARHI/P1H , PPH,U1H,U2H,V1H,V2H,W1H,W2H,R1H,R2H,RSH, 
8KEH,EPH,H1H,H2H,H3H,C1H,C2H,C3H,C4H,RXH, RYH , RZH , SI H , S2H 
COMMON/ GMTRY/VOL,VOLO,AEAST,ANORTH,AHIGH,AEDX, ANDY, AHDZ 
COMMON/PROP/ D1 , D2,D1DP,D2DP,MU1,MU1LAM, EXCO, CFP, MDT, HST1 , HST2 
COMMON/PRPOLD/DIO, D20 
COMMON/ PRPLOW/D1 L , D2L , EXCOL 
COMMON/ PRPHI/D1H, D2H,MU1H, EXCOH 
COMMON/VARNX/XG, XU, DXU, DXG 
COMMON/VARNY/YG, YV, DYV, DYG, R, RV 
COMMON/VARNZ/ZG, ZW1 , DZW, DZG, WGRID 
COMMON/GDMSCI/XPLANE, YPLANE, ZPLANE, ITNO 
COMMON/GDMSCL/LSLAB , MSLAB,HSLAB, LAMMU 



MEL 0217 0 
MEL02180 
MEL02190 
MEL02200 
MEL02210 
MEL02220 
MEL 02230 
MEL02240 
MEL 02250 
MEL 0226 0 
MEL0227 0 
MEL02280 
MEL 02290 
MEL02300 
MEL 0231 0 
MEL 02320 
MEL 0233 0 
MEL02340 
MEL 02350 
MEL 0236 0 
MEL 0237 0 
MEL02380 
MEL 0239 0 
MEL02400 
MEL 02410 
MEL 02420 
MEL 02430 
MEL 02440 
MEL 02450 
MEL 0246 0 
MEL 0247 0 
MEL02480 
MEL 02490 
MEL 02500 
MEL 02510 
MEL 02520 
MEL02530 
MEL 02540 
MEL02550 
MEL 0256 0 
MEL 0257 0 
MEL02580 
MEL 02590 
MEL02600 
MEL02610 
MEL02620 
MEL02630 
MEL 02640 
MEL02650 
MEL 0266 0 
MEL02670 
MEL02680 
MEL 02690 
MEL02700 
MEL 027 1 0 
MEL 027 20 
MEL 027 30 
MEL02740 
MEL02750 
MEL 027 6 0 
MEL 0277 0 
MEL02780 
MEL02790 
MEL 028 00 
MEL 02810 
MEL 02820 
MEL02830 
MEL 02840 
MEL02850 
MEL 0286 0 
MEL 0287 0 
MEL02880 
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REAL NORTH, LOW 

INTEGER P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS, 

&EP,H1,H2,H3,C1,C2,C3,C4,RX,RY,RZ,S1,S2 
INTEGER P10,PP0,U10,U20,V10, V20 , WIO , W20 , RIO , R20 , RSO, 

&EP0,H10,H20,H30,C10,C20,C30,C40,RX0,RY0,RZ0,S10,S20 
INTEGER P1L,PPL,U1L,U2L,V1L,V2L,W1L,W2L,R1L,R2L,RSL, 

&EPL,H1L,H2L,H3L,C1L,C2L,C3L,C4L,RXL,RYL,RZL,S1L,S2L 
INTEGER P1H,PPH,U1H,U2H, V1H, V2H, W1H, W2H , R1H, R2H, RSH, 

&EPH,H1H,H2H,H3H,C1H,C2H,C3H,C4H,RXH, RYH, RZH,S1H,S2H 
INTEGER VOL,VOLO, AEAST , ANORTH , AHIGH , AEDX , ANDY, AHDZ 
INTEGER D1,D1DP,D2,D2DP, EXCO ,CFP,HST1,HST2 
INTEGER DIO, D20,D1L,D2L, EXCOL , D1H, D2H, EXCOH 
INTEGER XG, XU, DXU, DXG, YG, YV, DYV, DYG, R, RV, ZG, ZW1 , DZW, 

&DZG, WGRID 

INTEGER XP LANE, YP LANE, ZP LANE 
LOGICAL LSLAB,MSLAB,HSLAB, LAMMU, LSPDA 
EQUIVALENCE (Ml , R1 ) , ( M2 , R2 ) 

SATLIT-EQUIVALEN7 IRUN: 

EQUIVALENCE ( IRUN, INTGRC 11 ) ) 

(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 1 ENDS. 
(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 1 STARTS: 

ARRAYS ( DIMENSIONED NY, NX ) FOR USE WITH ’ADD': 

DIMENSION CVARC 1 0 0,1),VVAR(1 00,1), CMC 100,1),VM( 100,1), ZERO (1,1) 
SPECIAL-DATA ARRAYS DIMENSIONED & DIMENSION VALUES SET HERE: 
DIMENSION LSPDA(l), ISPDA( 1 ) , RSPDA ( 1 ) 

USER PLACES HIS VARIABLES, ARRAYS, EQUIVALENCES ETC. HERE. 
DIMENSION GH1 ( 1 00 , 1 ) , GC1( 1 00 , 1 ) , GC2( 100 , 1 ) ,GC3( 100 , 1 ) 

DIMENSION GH2(100, 1) ,GD1( 100,1 ) ,GD2( 100,1 ),GR2( 10 0,1) 

DIMENSION GTEMP (100,1),GK0(3),GX1(3),GK2(3), GK3(3), GK4( 3) 
DIMENSION GMO ( 3 ) , GM1 ( 3 ) , GM2 ( 3 ) , GM3 ( 3 ) , GM4 ( 3 ) 

DIMENSION GCP0(3) ,GCP1(3) , GCP2( 3 ) , GCP3( 3 ) , GCP4(4) 

DIMENSION GK( 3 ) , GM(3), GCC 3 ) , 3X( 3) , GW( 3) , GWI ( 100 , 1 ) 

DIMENSION GD(4, 100,1 J , GVOL ( 100 , 1 ) , GV1 ( 1 00 , 1 ) 

DIMENSION GN( 100 , 1 ) , GYP(100,1), GYNV(100,1), GAN(100,1) 

DIMENSION GX0:i00,l) , GXMET ( 1 0 0 , 1 ) , GGAMH1 ( 1 00 , 1 ) , GMDT(100,1) 
DIMENSION GCPI (100,1), GDD( 4 , 4 ) ,GVISC( 100, 1 ) 

DIMENSION GCP(3, 100) ,GTEMP0(100, 1) 

NAMELIST /MOL E/ AT I ME, GW, GWI 
NAMELIST /TEMP/ ATIME, GTEMP, GXMET 

NAMELIST /PROPM/ ATIME, GKMF , GCPMF , GROF , GFOX , GHFU , GHFGOX, GHFGM, GEW 
NAMELIST /HT/GH200,GHOXO,GCPOX,GTMOX,GROOX,GLSOX,GTOXV, GTMVA , GTMM 
NAMELIST /KSPEC/GKO , GK1,GK2,GK3,GK4 
NAMELIST /MUSPEC/GMO , GM1 , GM2 , GM3 , GM4 
NAMELIST/DDSPEC/ GDD 

NAMELIST /CPSPEC/ GCPO , GCPI , GCP2 , GCP3, GCP4 
NAMELIST /CPSPI/ ATIME, GCP , GCPI , GX 
NAMELIST /SUMKUD/ GSUMK, GSUMU , GSUMD, GSUMF 
NAMELIST/YPV/ GVOL, GAN 
NAMELIST /CINI/ ATIME, GC1 , GC2, GC3 
DATA NLSP,NISP, NRSP/1, 1, 1/ 

DATA CVAR, WAR , CM, VM, ZERO/ 401*0.0/ 

USER PLACES HIS DATA STATEMENTS HERE. 



(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 1 

(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

GROUND INPUT FOR METAL COMBUSTION MODEL 



ENDS. 



(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DATA JING,GTCR,GQOX,GBOX,GDOX,GNOX, GMOX/O, 1700 . , 3 . E4, 1 .,l.E-6, 

SO. 5,1/ 

DATA GAQ, GKA/3 . E5, 1 . 0/ 

DATA GFST , GFG, GEG, GNMI , GVP , GSIG, GAO , GAI , GC1MIN/ 

S 0.99843, 0.11 1898, 0. 0001, 4. E08, 3. 2E-8, 5. 668E-8, 0.5, 0.5,0. 001/ 
DATA GPT, GKMF, GCPMF, GROF, GFOX, GHFU , GHFGOX, GHFGM, GH1M, GH3M, GEW/ 

S 100000. ,121. 3, 1083. 65,2340. ,1.886532, 4. 45E07, 4. 7 28 E06,1.13E07, 

S 0.0973E07,1.378E07,0.09/ 

DATA GH200 , GHOXO , GCPOX, GTMOX, GROOX, GLSOX, GTOXV , GTMM, GTMVA/ 
S-1.31E07,1 .643E07,1 27 1.936, 2323. ,4000. ,1 . 067E06 , 3253 ., 900 ., 2323 ./ 
DATA GK0(1),GK1(1),GK2(1),GK3(1),GK4(1),GK0(2),GK1(2),GK2(2), 
SGK3(2),GK4(2),GK0(3) , GK1 ( 3 ) , GK2( 3) , GK3( 3) , GK4( 3 )/ 

52133. 33. 26. 54. 0. . 0. . 0. . 125490. . 300. 172. 0..0..0., -20915. ,134. 9,0., 

50.. 0./ 

DATA GMO ( 1 ) , GM1 ( 1 ) , GM2 ( 1 ) , GM3 ( 1 ) , GM4 ( 1 ) , GMO ( 2 ) , GM1 ( 2 ) , GM2 ( 2 ) , 



MEL 02890 
MEL02900 
MEL02910 
MEL 02920 
MEL 02930 
MEL 02940 
MEL 02950 
MEL02960 
MEL 0297 0 
MEL02980 
MEL 02990 
MEL 03000 
MEL03010 
MEL 03020 
MEL 03030 
MEL 03040 
MEL 03050 
MEL 03060 
MEL 0307 0 
MEL 0308 0 
MEL 03090 
MEL03100 
MEL 03 1 1 0 
MEL03120 
MEL 0 3130 
MEL 03140 
MEL 03150 
MEL 03160 
MEL 0317 0 
MEL 03180 
MEL 03190 
MEL 03200 
MEL 0321 0 
MEL 03220 
MEL 03230 
MEL 03240 
MEL 03250 
MEL 0326 0 
MEL 0327 0 
MEL03280 
MEL 03290 
MEL 03300 
MEL 0331 0 
MEL 03320 
MEL 03330 
MEL 03340 
MEL 03350 
MEL 0336 0 
MEL 0337 0 
MEL 0338 0 
MEL 03390 
MEL03400 
MEL 0341 0 
MEL 03420 
MEL 03430 
MEL 03440 
MEL 03450 



MEL03460 
MEL 0347 0 
MEL03480 
MEL03490 
MEL03500 
MEL 0351 0 
MEL 03520 
MEL 03530 
MEL03540 
MEL03550 
MEL03560 
MEL 0357 0 
MEL 0358 0 
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8GM3( 3) , GM4( 2) , GM0(3),GM1(3),GM2(3), GM3( 3 ) , GM4( 3)/ MEL0359 0 

833., 0. 23395, 0., 0., 0., 100., 0.1 0 08857, 0., 0 ., 0 ., 56 ., 0.297 5, 0., 0., 0./ MEL03600 
DATA GCPO ( 1 ) , GCP1 ( 1 ) , GCP2C 1 ) , GCP3C 1),GCP4(1), GCPOC 2) , GCP1 ( 2 ) , MEL 036 10 

&GCP2( 2),GCP3(2), GCP4( 2),GCP0(3), GCP1 ( 3),GCP2(3), GCP3( 3),GCP4(3)/ MEL 036 20 
8 770.45,0.,0.,0.,0.,2307.5865,0.,0.,0.,0.,16497.4,0.,0.,0.,0./ MEL 0 36 30 
DATA GW(1),GW(2),GW<3)/ 26.98,18.016,2.01594/ MEL03640 

DATA GDD(1,1),GDD(1,2),GDD(1,3),GDD(2,1),GDD(2,2),GDD(2,3), MEL 03650 

8 GDD(3,1),GDD(3,2),GDD(3,3)/ MEL03660 

a 0.418,1.49,0.664,1.49,3.85,4.35,0.664,4.35,0.57/ MEL03670 

DATA GTEMP/23X9 0 0 .0, 150 0 . 0,150 0 . 0, 1800 .0,3X2323.0,2500 .0, MEL 03680 

8 2700.0,2900.0,3100.0,3500.0,3200.0, 3000.0,64X500.0/ MEL03690 

DATA GTEMPO/ 25x900 ., 10x2000 ., 25X900 ., 40x500 ./ MEL03700 

DATA GXMET/26X0. 0,0. 5, 0.75, 72X1.0/ MEL03710 

IF CATIME .GT. 0.0) GOTO 2900 MEL03720 

PRINT PROPM MEL037 30 

PRINT HT MEL03740 

PRINT KSPEC MEL037 50 

PRINT CPSPEC MEL 037 6 0 

PRINT MUSPEC MEL03770 

PRINT DDSPEC MEL03780 

GTMJL=GTMVA MEL03790 

GTMJS=GTMVA MEL03800 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 2 STARTS: MEL03810 

C PLEASE DO NOT ALTER, OR RE-SET, ANY OF THE REMAINING MEL03820 

C STATEMENTS OF THIS SECTION. MEL03830 

2900 IFCSPDATA) MEL03840 

aCALL RDSPCCIRN, INTGR(12), LSPDA, NLSP, ISPDA, NISP, RSPDA, NRSP) MEL 0 3850 

CALL GRDUTY( IRN, I CHAP ,IZED,INDVAR) MEL03860 

IF( ICHAP . EQ . -5) GO TO 10 MEL03870 

I F C I CHAP . L E . 0 . OR . I CHAP . GT . 16 ) RETURN MEL03880 

GO TO (100,200,300,400,500,600,700,800,900,1000,1100,1200, MEL03890 

81300,140 0,1500, 16 00), ICHAP MEL 0390 0 

RETURN MEL 03910 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 2 ENDS. MEL03920 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 2 STARTS: MEL03930 

c MEL 03940 

C CHAPTER 0: MODIFY SATLIT DATA, AT START OF EACH IRN. MEL03950 

c MEL0396 0 

10 CONTINUE MEL03970 

C END OF INITIAL CONDITION MEL03980 

C IF( . NOT . NAMLST ) RETURN MEL03990 

C IF(IRN.EQ.NRUN) DATFI L =. FALSE . MEL04000 

C— READ SATLIT DATA NAMELIST HERE MEL04010 

C CALL WRIT40(40HENTER NAMELIST DATA FOR GROUPS 1 TO 24 ) MEL04020 

C READC 20 , G1G24 ) MEL04030 

C CALL HRIT40C 40HENTER NAMELIST DATA FOR GROUPS 25 TO 42 ) MEL04040 

C READ(20,G25G42) MEL04050 

RETURN MEL04060 

c MEL04070 

C CHAPTER 1: CALLED AT THE START OF EACH TIME STEP. MEL04080 

C SET ’DT’ HERE WHEN TLAST SET NEGATIVE IN BLOCK DATA. MEL04090 

C ’ATIME + DT' GIVES THE END TIME OF THE CURRENT TIME STEP. MEL04100 

C NOT ACCESSED IF STEADY, OR PARABOLIC. MEL04110 

c MEL04120 

100 CONTINUE MEL04130 

ATIME=ATIME+DT MEL04140 

RETURN MEL 04150 

c MEL04160 

C CHAPTER 2: CALLED AT THE START OF EACH SWEEP. MEL04170 

c MEL04180 

200 CONTINUE MEL04190 

RETURN MEL04200 

c MEL 0421 0 

C CHAPTER 3: CALLED AT THE START OF EACH SLAB; MEL04220 

C NOT ACCESSED IF PARABOLIC, BUT ’STRIDE' IS. MEL04230 

c MEL04240 

300 CONTINUE MEL04250 

RETURN MEL 0426 0 

c MEL0427 0 

C CHAPTER 4: CALLED AT THE START OF EACH RE-CALCULATION OF MEL04280 

C VARIABLES P1,...C4 AT CURRENT SLAB. ITNO= ITERATION NUMBER. MEL04290 

c MEL04300 
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MEL 0 A 31 0 
MEL 04320 
MEL 04330 
MEL04340 
MEL 04350 
MEL 0436 0 
MEL 0437 0 
MEL 0438 0 
MEL04390 
MEL 0440 0 
MEL 04410 
MEL 04420 
MEL04430 
MEL 0 4440 
MEL 04450 
MEL 0446 0 
MEL04470 
MEL 0443 0 
MEL 0449 0 
MEL04500 
MEL 0451 0 
MEL 0 4520 
MEL04530 
MEL 04540 
MEL 04550 
MEL 0456 0 
MEL 0457 0 
MEL 0458 0 

exxxxxxxxxxxxx CHAPTER 5 xxxxxxxxxxxx**x*x*xx**xxxxx*xxxx MEL04590 

MEL 0 46 0 0 

(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXMEL046 1 0 
IMPUTATION OF SMOKE OXIDE CLOUD AND METAL SURFCE LOCATION AT NEW TMEL04620 
*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXMEL04o30 



CALL GET (Cl, GC1, NY, NX) 


MEL 046 40 


CALL GETCC2 , GC2, NY, NX) 


MEL 046 50 


CALL GETCC3, GC3, NY, NX) 


MEL0466 0 


CALL GET (ANORTH, GAN, NY, NX) 


MEL 0467 0 


CALL GET(R2,GR2, NY, NX) 


MEL 046 8 0 


CALL GET (VOLO,GVOL,NY,NX) 


MEL 046 90 


CALL GET1DCYG, GYP, NY) 


MEL 047 0 0 


CALL GET1DC YV , GYNV, NY) 


MEL 047 1 0 


CALL GET ( MDT , GMDT , NY, NX) 


MEL 04720 


CALL GET ( D1 , GDI , NY , NX ) 


MEL 047 30 


CALL GET ( D2 , GD2 , NY , NX ) 


MEL 047 40 


CALL GET ( HI , GH1 , NY, NX) 


MEL 04750 


CALL GET (H2, GH2, NY, NX) 


MEL 047 6 0 


CALL GET ( VI , GV1 , NY, NX) 


MEL 0477 0 


PRINT YPV 


MEL 04780 


IF ( ATIME .GT. DT) GOTO 2700 


MEL 04790 


JL = 26 


MEL 048 0 0 


JS = 28 


MEL 0481 0 


S PHASE ZONE PROPERTIES: K,MU, CPI, CP, D 


MEL 0 4820 


DO 2800 1=1, JL 


MEL04830 


WRITEC6,*) 'xxxxxxxxxxxx GTEMP= 1 , GTEMPC I , 1 ) , I 


MEL 0 484 0 


GT=GTEMP( 1,1) 


MEL 0 48 50 


GT2=GTXX2. 


MEL 0 486 0 


GT3=GT2XGT 


MEL 0487 0 


GT4=GT3XGT 


MEL04880 


GCPI (1,1 )=GCPMF 


MEL 0489 0 


GWI (1,1) =GW( 1 ) 


MEL 0490 0 


DO 2850 I N = 1 , 3 


MEL 049 10 


GCP(IN,I)=GCP0(IN)+GCP1(IN)XGT+GCP2( IN )XGT2+GCP3( IN) XGT 3 


MEL 04920 


&+GCP4(IN)XGT4 


MEL 04930 


CONTINUE 


MEL 04940 


CONTINUE 


MEL04950 


DO 60 I=JL , NY 


MEL 0496 0 


GT=GTEMP( 1,1) 


MEL0497 0 


GT2=GTXX2 . 


MEL 04980 


GT3=GT2XGT 


MEL 0499 0 


GT4=GT3XGT 


MEL05000 


IT = 1 


MEL 050 1 0 


DO 61 IN=1 , 3 


MEL 05020 



DO CONTINUE 
RETURN 



CHAPTER 5: GROUND CALLED WHEN SOURCE TERM IS COMPUTED. 

INDVAR GIVES DEPENDENT VARIABLE IN QUESTION IE. U1,...C4. 

TO ADD SOURCE TO DEPENDENT VARIABLE Cl(SAY) FOR IX=IXF,IXL 
AND IY=IYF,IYL INSERT STATEMENT: 

I F< INDVAR . EQ . Cl ) 

8CALL ADD( INDVAR, IXF, IXL , IYF, IYL , TYPE, CM, VM, CVAR, WAR, NY, NX) 
NOTES ON 'ADD* : 

XS0URCE= (CVARCIY,IX)+AMAX1(0 . 0 , MASFLO) )x( VVARC IY, IX)-PHI ) , 
WHERE 1 PHI 1 IS IN-CELL VALUE OF VARIABLE IN QUESTION. 

X * MASFL 0 1 = CMC IY, IX)*CVM( IY, IX)-P), 

WHERE * P * IS THE IN-CELL PRESSURE. 

xFOR INDVAR= Ml, OR =M2, SOURCE ADDED IS 1 MASFLO 1 ONLY, 

EXCEPT FOR ONEPHS= . F . 3 MASFLO < 0.0 CIE. OUTFLOW) WHEN 
CMC I Y, IX ) IS MULTIPLIED BY R1*D1 C FOR Ml) 8 R2*D2 C FOR M2). 

XBOTH 'CVAR' 3 'CM' ARE MUTLIPLIED BY CELL-GEOMETRY QUANTITY 
DICTATED BY SETTING OF 'TYPE* (=CELL, EAST AREA .VOLUME) . 

XTYPE-SP EC I FI ED AREAS ARE CALCULATED AS IF BLOCKAGE ABSENT, 

BUT •VOLUME' WITH ACCOUNT FOR ITS PRESENCE. 

XFOR ALL SOLVED VARIABLES, INCLUDING Ml C 3 M2 WHEN ONEPHS=F), 
IF ' CM ' > 0.0 CALL 'ADD*; FOR Ml 3 M2 ALTHOUGH 'CVAR' 3 1 WAR 1 
HAVE NO SIGNIFICANCE THEY MUST BE ENTERED AS ARGUMENTS. 

X • CVAR 1 , ' WAR 1 , 'CM* 8 'VM* MUST BE DIMENSIONED NY, NX. 



)0 CONTINUE 
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61 



4000 

3000 



3100 



3200 



GK(IN)=GK0(IN)+GK1(IN)XGT+GK2(IN)XGT2+GK3(IN)XGT3+GK4(IN)XGT4 
GM(IN) = (GMO( IN )+GMl ( IN ) XGT+GM2C IN ) XGT2+GM3C IN ) XGT3+GM4C IN ) XGT 4 )X 
81 . E-6 

GCP(IN,n=GCP0(IN)+GCPl(IN)XGT+GCP2(IN)XGT2+GCP3CIN)XGT3 
8+GCP4C IN)XGT4 
CONTINUE 

GCPI(I,1)=GCP(1,I)XGC1(I,1)+GCP(2,I)XGC2(I,1)+GCP(3,I)XGC3(I,1) 
KK = 1 

IF( (GW( 1 )XGW( 2) XGW( 3 ) ) .GT. 1.0E-27) GOTO 3000 
PRINT 4000,KK,I,GW(1) ,GW(2) ,GW(3) 

F0RMATC2I5,3F12.3) 

GWT=GC1(I,1)/GW(1) +GC2CI, 1)/GW(2) + GC3( I , 1 )/GW( 3) 

GWI ( I , 1 ) =1 ./GWT 
KK = 2 

IF ( GWI (1,1) .GT. l.E-27) GOTO 3100 
PRINT 40 00 , KK , I , GWI (1,1), GWT 



1)XGW(1)/GWI(I,1) 

1)XGW(2)/GWI(I,1) 

1)XGW(3)/GWI(I,1) 



) GOTO 3200 



3300 



3400 



GX( 1 ) =GC1 ( I i 
GX( 2) =GC2( I , 

GX( 3) =GC3( I i 
KK = 3 

IF ( GTEMP (1,1) .GT. 1.0E-27 
PRINT 4000,KK,I,GTEMP(I,1) 

GDI ( I , 1 ) =GPTXGWI ( I , 1 )/ ( 8483 . XGTEMP ( I , 1 ) ) 

GD2( I , 1 ) =GROOX 
GSUMK=0 . 

GSUMU= 0 . 

DO 4 IN=1 , 3 
GSUMF= 0 . 

GSUMD = 0 . 

DO 5 IK=1,3 
KK = 4 

I FC GDD( IN , IK) .GT. l.E-27) GOTO 3300 
PRINT 4000,KK,I,GDD(IN,IK) 

GSUMD=GSUMD+GX(IK)/(GDD(IN,IK)XGTEMP(I, 1)X*1.5) 

KK = 5 

IF( (GM( IK) X GWI IK) ) . GT . 1.0E-27) GOTO 3400 
PRINT 4000, KK, I,GM(IK) ,GW(IK) 

GMM=GM(IN)/GMCIK) 

GWR=GW(IN)/GW(IK) 

GXNK=(1 ./SQRT(8X(1+GWR) ))X(1+SQRT(GMM)/SQRT(SQRT(GWR) ))XX2. 

GSUMF=GSUMF+GMMXGXNK 

CONTINUE 

KK = 6 

IF ( GSUMFXGSUMDXGCPI (1,1) .GT. l.E-27) GOTO 3500 



MEL05030 
MEL 05040 
MEL05050 
MEL 0506 0 
MEL 0507 0 
MEL 0508 0 
MEL05090 
MEL05100 
MEL 051 1 0 
MEL05120 
MEL 05130 
MEL 05140 
MEL 05150 
MEL 0516 0 
MEL 0517 0 
MEL05180 
MEL 05190 
MEL 05200 
MEL 0521 0 
MEL 05220 
MEL 05230 
MEL 0524 0 
MEL 05250 
MEL 0526 0 
MEL 0527 0 
MEL05280 
MEL 05290 
MEL 05300 
MEL 0531 0 
MEL 05320 
MEL 05330 
MEL 05340 
MEL 05350 
MEL 0536 0 
MEL 0537 0 
MEL 05380 
MEL 05390 
MEL05400 
MEL05410 
MEL 05420 
MEL 05430 
MEL 05440 
MEL 05450 
MEL05460 





PRINT 4000, KK, I , GSUMF, GSUMD, GCPI ( I , 1) 


MEL0547 0 


3500 


GSUMK=GSUMK+GX(IN)/GSUMF 


MEL05480 




GSUMU=GSUMU+GX(IN)/GSUMF 


MEL05490 




GDINI=GX( IN)/ GSUMD 


MEL05500 




GD(IN,I,1)=GD1(I,1)XGDINI 


MEL05510 




CONTINUE 


MEL 05520 




GGAMH1 (1,1) =GSUMK/GCP 1(1,1) 


MEL05530 




GVISC(I,1)=GSUMU 


MEL 05540 


60 


CONTINUE 


MEL 05550 


CD 


PRINT MOLE 


MEL05560 


CD 


PRINT CPSPI 


MEL 0557 0 


CD 


PRINT SUMKUD 


MEL05580 


CD 


PRINT CINI 


MEL05590 




JLL= JL-1 


MEL 056 00 




DO 62 1=1, JLL 


MEL05610 


C FLUID PHASE PROPERTIES: ROO, MU, GAMMA=K/CP 


MEL05620 




GDI ( I , 1 ) =GROFXGCPMF/GCPI ( JL+1 , 1 ) 


MEL 056 30 




GD2 ( I , 1 ) =1 . E-l 0 


MEL05640 




GVISCC I , 1 ) = 10000. 


MEL05650 




GD( 1 , I , 1 ) =1 . E~1 0 


MEL 0566 0 




GDC 2, 1 , 1 ) =1 . E-l 0 


MEL 05670 




GD(3,I,1)=1.E-10 


MEL05680 




KK = 7 


MEL 05690 




IF ( GCPMFXGTMVA .GT. l.E-27) GOTO 3600 


MEL 057 00 




PRINT 4000 , KK, I , GCPMF , GTMVA 


MEL05710 


36 0 0 


GGAMH1 ( 1,1) =GKMF/GCPMF 


MEL05720 


62 


CONTINUE 


MEL 057 30 


C DENSITY AT THE SUFACE JL 


MEL 05740 
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GR0G=GPT*GWICJL,1)XGC1(JL,1)/(GTMJLX3483. ) 

KK = 8 

IF ( GROFX GCPOXXGROG . GT . l.E-27) GOTO 3700 
PRINT 4000, KK, I, GROF, GCPOX, GROG 
10 GD1( JL, 1) = 1 ./( 1 ./GROF+GXMET (JL+l,l)x(l ./GROG-1 . /GROF) ) 
GD2(JL,1)=GR00X 
IECKING IF IGNITION STARTS 
I F( JIGN .NE. 0) GOTO 2750 
IF ( GTEMP ( J L , 1 ) .GE. GTCR ) GOTO 2710 
INITION YET NOT STARTED 

,LCULATTION OF HEATING SOURCE AT THE METAL SURFACE 
DO 2720 IY=1,NY 
CM(IY, 1)=0 . 

VM(IY,l)+0 . 

CVAR(IY,1)=1 .E-10 

IF( IY .NE. JL ) GOTO 2720 

GD0XT=GB0XX(GC2( JL+1 )x*GMOX/GDOXX*GNOX)xEXP(-GAQ/GTEMP( JL, 1 ) ) 
VVARC IY, 1 ) =GROOXXGQGXXGDOXTXl . El 0 
GDOX=GDOX + GDOXTXDT 
0 CONTINUE 

CALL ADD (HI, 1, NX, 1 , NY, VOLUME, CM, VM,CVAR, WAR, NY, NX) 

GOTO 105 
MBUSTI ON STAGE 
0 JIGN = 1 
JSFLAG=0 . 

JRFLAG =0. 

GSUMN=0 . 

GSUMV = 0 . 

DO 111 1=1, NY 

IY=NY-I+1 

GT=GTEMP(IY,1) 

GTEMP 0( IY,1)=(GH2(IY,1) -GHOXO-GCPOXX298 )/GCPOX 
GCPDT=GCPI(IY, l)x( GTEMP (IY, 1) -293) 

GCH20=GC2(IY,1)*GH200 

GCHFG=GC1(IY,1)X(GCP(1,IY)-GCPMF) 

KK = 9 

IF ( GVPXGCPMF .GT. l.E-27 ) GOTO 3300 
PRINT 4000,KK,I,GVP,GCPMF 

0 GN(IY, l)=GVOL(IY, 1 )XGR2( IY, 1 )/GVP 
IF ( GN( IY, 1 ) .LT. GNMI ) GOTO 11 

0. OF OXIDE PARTICLES BECOME SIGNIFICANT 
JCO= IY 
JRFLAG=1 
GOTO 12 

1 IF (JRFLAG .NE. 1) GOTO 13 

IF ( GN ( IY, 1 ) .GT. GNMI) GOTO 12 

.OF OXIDE PARTICLES BECAME INSIGNIFICANT AGAIN, GOING OUT OF CLOUD 
JCI =IY 
JRFLAG=0 . 

GOTO 13 

2 GSUMN=GSUMN+GN( IY , 1 )XGVOL ( IY, 1 ) 

GSUMV=GSUMV+GVOL ( IY, 1 ) 

3 I F( GC1 ( IY, 1 ) .LE. GC1MIN) G0T013 
GPAL = ( GC1 ( IY, 1 ) XGWI ( IY, 1 )/GW( 1 ) ) *GPT 
GTMVA= 36560 .8/(25. O16-ALOG( GPAL)) 

GXMETA=(GH1(IY,1)-(GH1M+ GCPMF^(GTMVA-GTMM) )*GC1 ( IY, 1 ) - 
8 ( GH200+ GCP( 2, IY)*(GTMVA-37 3.0) )^GC2( IY, 1 ) - 
& GCP(3,IY)X(GTMVA-298.0)*GC3(IY,1))/ ( GHFGMXGC1 ( IY, 1 ) ) 

I F( GXMETA .LT. 0.0) GXMETA=0.0 
IF ( IY .GT. JS ) GOTO 18 
IF ( JSFLAG .EQ. 0. ) GOTO 14 

GXMETA=(GH1( IY,1)-(GH1M+ GCPMFX(GTMJS-GTMM) )XGC1( IY, 1 ) - 
8 (GH2O0+ GCP(2,IY)X(GTMVA-373.0))XGC2(IY,1) - 
8 GCP(3,IY)X(GTMVA-298.0)XGC3(IY,1))/ ( GHFGMXGC1 ( I Y, 1) ) 

IF ( JSFLAG .GT. 1 ) GOTO 17 
FLAG = 1 WHERE 0<= X =>1 

IF ( GXMETA .GT. 0. ) GOTO 20 
=0 FOR THE FIRST TIME 
J L = IY 
JSFLAG=2 
GXMET ( IY, 1 ) = 0 . 0 
GTMJL=GTMVA 



MEL 0 57 50 
MEL 0576 0 
MEL 0577 0 
MEL 057 8 0 
MEL 05790 
MEL 05800 



MEL 058 1 0 
MEL 05820 
MEL 05830 
MEL 05840 
MEL 05850 
MEL 0586 0 
MEL 0587 0 
MEL 0588 0 
MEL 0589 0 
MEL 059 00 
MEL05910 
MEL 05920 
MEL05930 
MEL 05940 
MEL 05950 
MEL0596 0 
MEL 0597 0 
MEL05980 
MEL 0599 0 
MEL06000 
MEL0601 0 
MEL 06 020 
MEL06030 
MEL06040 
MEL06050 
MEL 06 06 0 
MEL 06 07 0 
MEL 06 080 
MEL 06 09 0 
MEL 06 1 0 0 
MEL 06 1 1 0 
MEL 06 120 
MEL 06 130 
MEL 06 140 
MEL 06 1 50 
MEL 06 16 0 
MEL 06 17 0 
MEL 06 18 0 
MEL 06 190 
MEL06200 
MEL 06 210 
MEL 06 220 
MEL 06 230 
MEL 06 240 
MEL 06 250 
MEL 06 26 0 
MEL 06 27 0 
MEL 0623 0 
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GTMJS=GTMVA 
GTEMPCIY, 1 )=GTMJL 
GOTO 20 

17 GTEMPCIY, 1)=CGH1CIY,1)-GH1M+GCPMF*GTMM)/GCPMF 

GXMET ( IY, 1 )=0 . 0 
GOTO 111 

CHECKING IF SATURATION OCCOURE 
1 A KK=10 

IF CGPTXGTMVAXGEW .GT. 1 . E-27 ) GOTO 3900 
PRINT 4000,KK,GPT, GTMVA, GEW 
3900 CONTINUE 

IF ( GXMET A .LE. 1 ) GOTO 19 
X>1, GAS PHASE 

18 GXMET ( IY, 1 ) =1 . 

GTEMPCIY, 1)=CGH1CIY,1)-CGH3M-GCPC1,IY)XGTMJS)XGC1CIY,1) - 
8 CGH200-GCPC2,IY)X373.0)XGC2CIY,1) + GCPC 3, IY)*298 . 0*GC3C IY, 1 ) ) / 
& GCP I C IY , 1 ) 

GOTO 111 

X<=1 , FOR THE FIRST TIME AT THIS SWEEP 

19 JS=IY 
GTMJS=GTMVA 
JSFLAG=1 

20 GCPI F = GCPI C IY, 1 )-C GCPI C IY, 1 )-GCPMF) XGC1 C IY, 1 ) 

IF C GCPI F .GT. l.E-10) GOTO 22 

WRITEC6,X)'XXXGCPIF,GCPI,GC1, IY', GCPI F, GCPI CIY,1),GC1CIY,1),IY 
22 GTEMPC IY, 1 ) =CGH1C IY, 1 ) - C GHlM-GCPMFxGTMM+GXMETA*GHFGM) XGC1 C IY, 1 ) 

& - CGH200-GCPC2,IY)X373.0)XGC2CIY,1)+ GCPC 3, IY)XGC3 CIY, 11X298.0) 

8 / GCPI F 

GXMET C IY, 1 ) = GXMET A 
111 CONTINUE 
PRINT TEMP 

WRITEC6 ,x) 'XXJL= JS = * , JL , JS, 'xxxTMJL TMJS= ' ,GTMJL,GTMJS 
GTW= GTEMPC JL, 1 ) 

GTW4 = GTWXXA 

IF C GSUMN .LE. l.E-10 ) GOTO 2950 
RADIATION LOSES FROM THE OXIDE CLOUD 
KK=1 05 

I FCGSUMV .GT. 1 . E-27 ) GOTO 5000 
PRINT 4000,KK, JCO, JCI,GSUMV 
5000 GNN=GSUMN/GSUMV 

GL = GYPC JCO, 1 )-GYPC JCI , 1 ) 

GRC=GYPCJCI,1)+0.5XGL 
GEP=CA.E-12)XGNN 
GRPW=GYPC J L , 1 ) 

GEP=1 .-EXPC-GEPXGL) 

GFCW=1 ./Cl./GEW+l ./GEP-1) 

JC=C JCO+JCI )/2 
GT OX=GTEMP C JC , 1 ) 

GT0X4= GTOXxxA 
JG=C JC0+NY1/2 
GTG=GTEMP C JG, 1 ) 

GTG4=GTGxx4 

GAP=GEP 

GSIG=5 . 6688E-8 
GAO=0 . 5 
GAI = 0 .5 

GRFC=CGAPXCGFCWXGAI*GEWXGTWA+GEG*GTG4XGA0+GEPXC1-GEW)XGFCWXGAIX 
a gtoxai-gepxgtoxa)xgsig 

KK=1 06 

IFCGL .GT. 1.0E-27 ) GOTO 5100 
PRINT A000,KK, JCI,GL 
5100 GRCI =-GRFC/GL 
2950 CONTINUE 

WRITEC 6 , X ) ' XXBEGIN SOURCE, INDVAR TIME SWEEP ', INDVAR, ATIME, LSWEEP 
GOTO Cl 05, 105, 105, 105, 105, 105, 105, 105, 105, 105, 105, 105, 105, 

8 401, 501, 105, 601, 701, 801, 105) INDVAR 
SOURCE TERM FOR C2 i OXIDIZER 
701 DO 710 IY=2, NY 
CMC IY, 1 ) = 0 . 

VMC IY, 1 ) = 0 . 

CVARCIY,1)=1.E-10 
GQ=GAQ/GTEMPC IY, 1 ) 



MEL06290 
MEL06300 
MEL06310 
MEL06320 
MEL06330 
MEL06340 
MEL 06 350 
MEL06360 
MEL 06 37 0 
MEL06380 
MEL06390 
MEL06400 
MEL06410 
MEL06420 
MEL06430 
MEL06440 
MEL06450 
MEL06460 
MEL 0647 0 
MEL06480 
MEL06490 
MEL06500 
MEL06510 
MEL06520 
MEL06530 
MEL06540 
MEL 06 550 
MEL 06 56 0 
MEL06570 
MEL06580 
MEL06590 
MEL06600 
MEL06610 
MEL06620 
MEL 06630 
MEL06640 
MEL06650 
MEL06660 
MEL 0667 0 
MEL 06680 
MEL06690 
MEL067 00 
MEL 0671 0 
MEL06720 
MEL 067 30 
MEL 067 40 
MEL 067 50 
MEL06760 
MEL 0677 0 
MEL 0678 0 
MEL06790 
MEL 068 0 0 
MEL06810 
MEL 06820 
MEL06830 
MEL06840 
MEL06850 
MEL06860 
MEL06870 
MEL06880 
MEL06890 
MEL06900 
MEL 06 9 1 0 
MEL 06920 
MEL06930 
MEL06940 
MEL 06950 
MEL06960 
MEL06970 
MEL06980 
MEL06990 
MEL07000 
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IF (GQ .GE. 5. ) GQ = 5. 

GKR=GKA*EXP(-GQ) 

GCDT = -3XGKR*GC2dY,l)XX3.XGCl(IY,l)XX2. 

IF( GV1 ( IY, 1 ) .LT. 0.0 ) GOTO 750 
GROU=(GD1(IY,1)+GD1(IY-1,1))*0.5*GV1(IY-1,1) 
GYN=GYP(IY,1)-GYP(IY-1, 1) 

GR0U0X=GR0UX(GC2(IY, 1 )-GC2(IY-l , 1 ) )/GYN 
GOTO 7 AO 

>0 GR0U=(GD1(IY, D+GDKIY+l, 1) )X0 .5*GV1(IY, 1) 

GYN=GYP(IY+1, 1)-GYP(IY, 1) 
GR0U0X=GR0UX(GC2(IY+l,l)-GC2dY,l))/GYN 
(0 IF( (ABS(GCDT)-ABS(GROUOX) ) . GT . l.E-5) GOTO 720 
THE REACTION IS CONTROLED BY THE RATE OF REACTION. 
VVAR(IY,l)=GCDTXl .E10 
PRINT 4000, GCDT 
100 FORMAT (5HGCDT=,F15.3) 

GMDTCIY, 1)=-(1 ./3. )XVVAR(IY, l)xl . E-10 
GOTO 710 

HE REACTION IS CONTROLED BY DIFFUSION FLUX. 

20 GF=GC1(IY,1)/GC2(IY,1) 

IF(GF .LE. GFST ) GOTO 730 

IIFFUSION FUEL RICH. 

GROUOX=-ABS(GROUOX) 

GMDT(IY,l)=-GROUOXXGFOX 
VVARCIY,l)=GROUOXXl .E10 
PRINT 4100, GROUOX 
00 FORMAT (7HGR0U0X=,F15.3) 

GOTO 710 
FUEL LEAN 

30 VVARdY, 1)=-ABS(GR0UX(GC1CIY, 1)-GC1(IY-1, l))/GYN)xl .E10 
GR0UC1=VVAR(IY,1)*1 .E-10 
PRINT 4200 , GROUC1 
! 0 0 FORMAT (7 HGROUC1 = , FI 5 . 3 ) 

10 CONTINUE 

VVAR( IY, 1 ) = 0 . 0 

CALL ADD(C2,1,NX,1,NY,V0LUME,CM,VM,CVAR,VVAR,NY,NX) 
WRITE(6,x) 'XXC2, 1,25, 28,35, 100 * , GC2 
WRITEC6 , *) ' VV C2 ' , WAR 

WRITEC6 ,30 'xxxxxx END SOURCE SINDVAR 

WRITEC6,*) 'XXENTHALPY' ,GH1 

RETURN 

iOURCE TERM FOR Cl : FUEL 
■01 DO 610 IY=2, NY 
CM( IY, 1 ) =0 . 

VM( IY, 1 ) =0 . 

CVARCIY, 1)=1 .E-10 
GQ=GAQ/GTEMP( IY, 1 ) 

IF (GQ .GE. 5. ) GQ=5 . 

GKR=GKAxEXP(-GQ) 

GCDT=-2xGKRxGC2( IY, 1 )XX3 . *GC1 ( IY, 1 ) XX2 . 

I F( GV1 ( IY , 1 ) .LT. 0.0 ) GOTO 650 

GROU = (GD1(IY,1)+GD1(IY-1,1))5<0.5*GV1(IY-1,1) 

GYN=GYP( IY , 1 ) -GYP( IY-1 , 1 ) 

GR0UCM=GR0U*(GC2( IY, 1 )-GC2( IY-1 , 1 ) )/GYN 
GOTO 640 

>0 GROU = (GD1(IY,1)+GD1(IY+1,1))3<0.5XGV1(IY,-1) 

GYN=GYP( IY+1 , 1 )-GYP( IY, 1 ) 

GR0UCM=GR0U*(GC2(IY+1 , 1 )-GC2C IY, 1 ) )/GYN 
iO IF( (ABS(GCDT)-ABS(GROUCM) ) . GT . l.E-5) GOTO 620 
THE REACTION IS CONTROLED BY THE RATE OF REACTION. 
VVAR(IY,1)=GCDT*1 .E10 
GMDT(IY,1)=-(1 ./3. )*VVAR(IY,1)5(1 .E-10 
GOTO 610 

rHE REACTION IS CONTROLED BY DIFFUSION FLUX. 

>20 GF=GC1(IY,1)/GC2(IY,1) 

IF(GF .LE. GFST) GOTO 630 

)IFFUSION FUEL RICH. 

GR0U0X=-ABS(GR0UX(GC2(IY,1)-GC2(IY-1,1) )/GYN) 

VVAR(IY, 1)=GR0U0X3(GFSTX1 .E10 
GOTO 610 
FUEL LEAN 
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MEL07070 
MEL070S0 
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MEL07100 
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MEL 07 120 
MEL07130 
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MEL07150 
MEL07160 
MEL07170 
MEL 07 13 0 
MEL07190 
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MEL 07 220 
MEL 07 230 
MEL07240 
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MEL 07 27 0 
MEL07230 
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MEL 07 320 
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MEL07590 
MEL07600 
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MEL07650 
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630 VVARC IY,1) = -ABSCGR0UCM)X1. El 0 
GMDTCIY, 1)=-GR0UCMXGF0X/GFST 
610 CONTINUE 

VVARC IY, 1 ) =0 . 0 

CALL ADDC Cl , 1 , NX, 1 , NY, VOLUME, CM, VM, CVAR, WAR, NY, NX) 
WRITEC6,x)' 3€XC 1 ,1,25,28,35,100' ,GC1 
WRITEC 6 , X ) ' VV Cl ' , WAR 

WRITEC6,x) 'xxxxxx END SOURCE ' , INDVAR 

WRITEC 6 , X) ' xxENTHALPY ' , GH1 

RETURN 

C SOURCE TERM FOR C3 : HYDROGEN GAS PRODUCT 
801 DO 810 IY=2, NY 
CMC IY, 1 ) =0 . 

VMC IY, 1 ) =0 . 

CVARC IY,1)=1.E-10 
GQ=GAQ/GTEMP( IY, 1 ) 

IF (GQ .GE. 5. ) GQ=5 . 

GKR=GKAXEXP(-GQ) 

GCDT=-2XGKRXGC2( IY, 1 )xx3 . XGC1 ( IY, 1 )XX2 . 

IFC GV1 C IY, 1 ) .LT. 0.0 ) GOTO 850 
GR0U=(GD1(IY,1)+GD1(IY-1, 1))X0.5XGV1(IY-1,1) 
GYN=GYPCIY,1)-GYP(IY-1, 1) 

GR0UCM=GR0UX(GC2( IY, 1 )-GC2( IY-1 , 1 ) )/GYN 
GOTO 8 AO 

850 GR0U=(GD1CIY,1)+GD1(IY+1, 1))X0.5XGV1(IY,1) 

GYN=GYP ( IY+1 , 1 )-GYP< IY, 1 ) 

GR0UCM=GR0UX(GC2( IY+1 , 1 )-GC2( IY, 1 ) )/GYN 
8 AO IFCCABS(GCDT)-ABSCGROUCM)) . GT . l.E-5) GOTO 820 
C THE REACTION IS CONTROLED BY THE RATE OF REACTION. 
WAR(IY,1)=GCDTX1.E10 
GOTO 810 

C THE REACTION IS CONTROLED BY DIFFUSION FLUX. 

820 GF=GC1CIY,1)/GC2(IY,1) 

IFCGF .LE. GFST) GOTO 830 

C DIFFUSION FUEL RICH. 

GR0U0X=-ABS(GR0UX(GC2(IY, D-GC2CIY-1, 1))/GYN) 

VVARC IY , 1 )=-GROUOXXGFGXl . E10 
GOTO 810 
C FUEL LEAN 

830 VVARC IY, 1 ) =-ABSC GROUCM) XC GFG/GFST ) XI . E10 
810 CONTINUE 

VVARC IY, 1 )=0 . 0 

CALL ADDCC3, 1 , NX, 1 , NY, VOLUME, CM, VM, C VAR, WAR, NY, NX) 
WRITEC6,x) *XXC3,1, 25, 28,35,100' ,GC3 
HRITEC6 , X) ' VV C3 ' , WAR 

WRITEC 6 > x) »x**x** END SOURCE ', INDVAR 

WRITEC 6 , X ) 'XXENTHALPY' ,GH1 

RETURN 

C HEAT SORCE TERM FOR HI 
A01 DO A10 I Y=1 , NY 
CMCIY, 1 )=0 . 

CVARCIY,1)=1.E-10 
GT =GTEMPC IY, 1 ) 

GPRP=GGAMH1 C IY, 1 ) 

GPRS=GGAMH1 C IY-1 , 1 ) 

GPRN=GGAMH1 C IY+1 , 1 ) 

GYNN=GYPC IY+1 , 1 ) -GYNVC IY, 1 ) 

KK = 12 

IFCGANC IY, 1) .GT. l.E-27)G0T0 A200 
PRINT A000,KK,IY,GANCIY,1) 

A200 GYN = GYNVCIY, D-GYPCIY, 1) 

GYS=GYPC IY, 1 )-GYNVC IY-1 , 1 ) 

GYSS=GYNVCIY-1, D-GYPCIY-l, 1) 

GPRNX=GPRNXC 1-GXMETC IY+1 , 1 ) ) 

GPRPX=GPRPXC 1-GXMETC IY, 1 ) ) 

GPRSX=GPRSXC 1-GXMETC IY-1, 1)) 

GPRNC=GPRNXGC1 C IY+1 , 1 ) 

GPRPC=GPRPXGC1CIY, 1) 

GPRSC=GPRSXGC1CIY-1, 1) 

KK = 13 

IFCGPRNXXGPRPXXGPRSX .GT. 1 . E-27 ) GOTO A300 
PRINT A000 , KK, IY, GPRNX, GPRPX, GPRSX 
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MEL08090 
MEL081 00 
MEL 081 1 0 
MEL 08120 
MEL08130 
MEL081A0 
MEL08150 
MEL 0816 0 
MEL08170 
MEL08180 
MEL08190 
MEL08200 
MEL 08 21 0 
MEL08220 
MEL08230 
MEL082A0 
MEL08250 
MEL08260 
MEL0827 0 
MEL08280 
MEL08290 
MEL08300 
MEL08310 
MEL08320 
MEL08330 
MEL083A0 
MEL08350 
MEL 08 36 0 
MEL 08 37 0 
MEL08380 
MEL08390 
MEL08A00 
MEL08A10 
MEL08A20 
MEL08A30 
MEL08AA0 



26 



MEL3GRD FORTRAN A1 



KK = 1 4 

IFCGPRNCxGPRPCXGYNN . GT . l.E-27) GOTO 4400 
PRINT 4000,KK,IY,GPRNC,GPRPC,GYNN 
KK = 15 

IF(GPRNXGPRPXGCP(1,IY+1) . GT . l.E-27) GOTO 4500 
PRINT 4000,KK, IY, GPRN, GPRP, GCPC 1 , IY+1 ) 

KK = 16 

IF(GPRSXGCP(1,IY-1)XGCP(2,IY) .GT. l.E-27) GOTO 4600 
PRINT 4000, KK,IY, GPRS, GCPC 1, IY-1) , GCPC 2, IY) 

KK= 17 

IFCGYSXGYSSXGCPC3, IY) . GT . l.E-27) GOTO 4700 
PRINT 4000,KK,IY,GYS,GYSS,GCP(3, IY) 

KK = 18 

IF(GCP(3,IY-1)XGPRNXGYN . GT . l.E-27) GOTO 4800 
PRINT 4G00,KK, IY, GCPC 3, IY-1 ) , GPRN, GYN 
IF ( IY .NE. JL ) GOTO 420 
OF LIQUID METAL 
IS THE VALUE OF RADIATION 
GRW=-GEWXGSIGXGTW4/GYN 
IF ( GSUMN .LE. 1.E-1Q ) GOTO 460 

GRW=(GEPxGT0X4XGFCW+GEGxGTEMP( ( JL+JCI )/2, 1)XX4)XGEW/GYN+GRW 
0 VVARCIY, 1)=GRWX1 . E10 
GOTO 430 

0 IFCCIY .GT. JS) .OR. (IY .LT. JL)) GOTO 430 
GION WHERE X<1 METAL SATURETI ON 

GSXM1=CGC1CIY+1,1)-GC1CIY,1))/CGYNN/GPRNX+GYN/GPRPX)-CGC1CIY,1)- 
SGClCIY-1, 1 ) )/ ( GYS/GPRNX+GYSS/GPRSX ) 

GSXM2=(GXMETCIY+1,1)-GXMET( IY, 1 ) )/ C GYNN/GPRNC+GYN/GPRPC) - 
8 ( GXMET ( I Y, 1 )-GXMETC IY-1 , 1 ) )/ ( GYS/GPRPC+GYSS/GPRSC ) 

GSXM=GHFGMx ( GSXM2-GSXM1 ) 

VVARCIY, 1)=GSXM 
0 IF ( IY .LT. JL ) GOTO 440 
EAT RELEASED 3Y METAL 
HOXV FOR HEAT RELEASED 



RFCE 
W 



ABSORBED AT THE SURFACE 



MEL 08450 
MEL 0846 0 
MEL 0847 0 
MEL 0848 0 
MEL 08490 
MEL 0850 0 
MEL 08510 
MEL 08520 
MEL 08530 
MEL 08540 
MEL 08 550 
MEL 0856 0 
MEL 0857 0 
MEL 08580 
MEL 0859 0 
MEL 036 0 0 
MEL 086 1 0 
MEL 086 20 
MEL 086 3 0 
MEL 086 40 
MEL 036 50 
MEL 0866 0 
MEL 086 7 0 
MEL 0863 0 
MEL 086 9 0 
MEL 037 CO 
MEL 037 1 0 
MEL08720 
MEL 087 30 
MEL 087 4 0 
MEL 087 50 



MEL 0876 0 

COMBUSTION + LATENT CV3S) HEAT OF OXIDE IY>JMEL08770 
BY VAPOR PHASE OXIDE CHANGE TO LIQUID OXIDE MEL08730 



RELEASED BY LIQUID PHASE OXIDE CHANGE TO SOLID OXIDE 



GTOXV ) 
GTOXV ) 

GTMOX) 

GTMOX) 



GHOXV=GMDT ( IY, 1 )XGHFGOX 
GHOXV=GHOXV+ GR00XxGR2(IY,l)x 

GH0XS=GMDT(IY,1)*GLS0X 
GHOXS=GHOXS+ GR00XxGR2(IY,l)x 



GHOXV 
THE OXIDE 



PARTICLES 



HOXS FOR HEAT 
GHOXV=0 . 0 
GHOXS=0 . 

IF C GTEMP C IY, 1 ) .LE. 

I FC GTEMPO C I Y , 1 ) .LE. 

&GV0L(IY,1)XGHFG0X 
IF C GTEMPC IY, 1 ) .LE. 

I FCGTEMPOC IY, 1 ) .LE. 

8GV0L(IY,1)*GLS0X 

VVARCIY, 1)=VVAR(IY,1)+GMFU*GHFU + GHOXS + 

AT TRANSFERED BY CONVECTION TO THE GAS FROM 
GAS=3.xGR2(IY,l)/2. E-6 

VVAR(IY,1)=VVAR(IY,1) + (5.23E06)XGASX(GTEMPO(IY, 1 ) -GTEMP ( IY, 1 ) ) 
RRECTION OF HEAT SOURCE DUE TO ENTHALPY CONDUTION 

GDHM=CGC1CIY+1,1)-GC1(IY, 1))/CGYNN/GPRNXGCP(1, IY+D+GYN/CGPRPx 
&GCPC1, IY) ))-(GCl(IY, 1 )-GCl(IY-l, 1) )/( GYS/ ( GPRPXGCP ( 1 , IY) )+GYSS/( 
8GPRSXGCPC 1 , IY-1 ) ) ) 

GDH2=CGC2C IY+1 , 1 )-GC2( IY, 1 ) )/(GYNN/(GPRNxGCP( 2, IY+1 ) )+GYN/(GPRP 
&XGCPC2, IY) ))-(GC2CIY,l)-GC2CIY-l, 1))/CGYS/(GPRPXGCPC2, IY) )+ 

8 GYSS/(GPRS*GCP(2,IY-1))) 

GDH20=CGC3C IY+1 , 1 )-GC3( IY, 1 ) )/ C GYNN/ ( GPRNXGCP C 3 , IY+1 ) )+GYN/(GPRPX 
8GCPC 3 , IY) ) ) - ( GC3 ( I Y, 1 ) -GC3 ( IY-1, l))/( GYS/ ( GPRPXGCP ( 3, IY))+GYSS/ 
a ( GPRSXGCP C 3 IY— 1 ) ) ) 

GH20 = C GC3(IY+1,1) -GC3 ( I Y, 1 ) )/ C GYNN/GPRN+GYN/GPRP ) - 
a ( GC3C IY, 1 ) -GC3 C IY-1 , 1 ) )/ C GYS/GPRP+GYSS/GPRS ) 

VVARCIY, 1)=VVAR(IY,1)+(GDHM+ GDH2 +GDH20 ) x( GT-298 )+ GH20QXGH20 
GOTO 450 

AT SOURCE IN METAL . NO HEAT GNERATEDAT THIS RIGION IY<JL. 

0 VVARCIY, 1)=0.0 
0 VVARCIY, 1)=VVARCIY, 1)X1 .E10 
0 CONTINUE 

CALL ADDCH1,1,NX,1,NY,V0LUME,CM,VM,CVAR,VVAR,NY,NX) 

WRIT EC 6 , X ) 1 xxxHl 1 ,GH1 

WRITE C 6, X) ' XXTEMP' , GTEMP 

WRITEC 6 , x ) 1 XX HI SOURE',VVAR 

WRITEC 6 , x ) »xxxxxx END SOURCE 1 , INDVAR 

RETURN 

AT SOURCES ON THE OXIDE PHASE H2 



MEL 087 9 0 
MEL 088 0 0 
MEL 0881 0 
MEL 08820 
MEL 088 30 
MEL 08840 
MEL 08850 
MEL 0886 0 
MEL 0887 0 
MEL 0888 0 
MEL 0889 0 
MEL 08 90 0 
MEL0891 0 
MEL 08920 
MEL 089 30 
MEL 0894 0 
MEL 08 950 
MEL 0896 0 
MEL 0897 0 
MEL 0898 0 
MEL 0899 0 
MEL09000 
MEL 09010 
MEL 09 020 
MEL 09030 
MEL 09040 
MEL 09 050 
MEL 09 06 0 
MEL 09 07 0 
MEL 0908 0 
MEL09090 
MEL 091 0 0 
MEL 091 1 0 
MEL09120 
MEL 0 9130 
MEL09140 
MEL 09 150 
MEL09160 
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501 DO 510 IY=1,NY MEL09170 

KK=19 MEL09180 

I F ( GCPOXXGCPI ! IY, 1 ) XGVOL I IY, 1 ) .GT. l.E-27) GOTO 9900 MEL09190 

PRINT 9000,KK,IY,GCPOX,GCPI!IY,1),GVOL!IY,1) MEL09200 

9900 CM(IY,1) = 0. MEL 09210 

VM! IY, 1 ) =0 . MEL 09220 

GH0XS = 0 . 0 MEL 09230 

CVAR!IY,1)=1.E-10 MEL 09290 

WAR! IY, 1 ) = 0 . MEL 09250 

510 CONTINUE MEL09260 

DO 520 IY=JL,NY MEL09270 

GAS=3.*GR2(IY,l)/2.E-6 MEL09280 

C HEAT LOSES BY CONVECTION TO THE GAS PHASE MEL09290 

VVAR(IY,1) = -(5.23E6)xgAS*(GH2(IY, 1)XGCP! IY, 1 ) -GH1 I IY, 1 ) XGCP0X+ MEL 09 300 
8 GCP0XX!GC3! IY,1)XGH200) -GCPI!IY,1)XGH0X0)/ ! GCPOXXGCPI !IY, 1) ) MEL 09 310 

IF!!!GTEMPO!IY,l)-GTMOX) .LT. 3) .OR. (GTEMP(IY,1) .LE. GTMOX) ) MEL09320 

8 GH0XS=GR00XXGR2(IY,1)*GV0LCIY,1)*GLS0X MEL 09330 

C SOLDIFICATION OF OXIDE - HEAT LOSES TO GAS PHASE MEL09390 

WAR! IY,1)=VVAR! IY,l)-GHOXS MEL 0935 0 

C RADIATION HEAT LOSES MEL09360 

IF( ( IY .GE.JCI) .AND. C IY .LE. JCO)) WAR! IY, 1 ) =VVAR( IY, 1 ) -GRCI MEL09370 
C CONTRIBUTION OF HEAT GENERATED BY NEW OXIDE COMBUSTION PRODUCT MEL09380 

WAR! IY,1)=VVAR!IY, 1)+GMDT!IY, l)XGCPOXX!GTOXV-GTEMPO!IY, 1))/ MEL 09390 

SGVOL(IY,l) MEL09900 

VVARCIY, 1)=VVAR! IY,1)X1.E10 MEL 09910 

520 CONTINUE MEL09920 

CALL ADD(H2,1,NX,1,NY,V0LUME,CM,VM,CVAR,VVAR,NY,NX) MEL09930 

WRITE! 6 , X) • XXXH2 ' , GH2 MEL09990 

WRITE!6,X) 'XX TEMP 0X',GTEMP0 MEL09950 

WRITE!6,x) 'xx H2 SOURE',VVAR MEL09960 

WRITE! 6 , X ) *XXXXXX END SOURCE ' , INDVAR MEL09970 

105 RETURN MEL09980 

CXXXXXXXXXXXXXXXEND OF CHAPTER 5 XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXMEL 09990 

MEL 0950 0 

CHAPTER 6: CALLED AT THE END OF EACH VARIABLE-RECALCULATION MEL09510 

CYCLE COMMENCED AT CHAPTER 9. ITNO = ITERATION NUMBER. MEL09520 

MEL 09530 

600 CONTINUE MEL09590 

RETURN MEL09550 

MEL 09560 

CHAPTER 7: CALLED AT END OF EACH SLAB-WISE CALCULATION. MEL09570 

MEL09580 

700 CONTINUE MEL09590 

CXXXXXXXXXXXXXXX CHAPER 7 XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXMEL 09600 
CALL SET! U 2, 1 , NX, 1 , NY, GTEMP, NY, NX) MEL09610 

CALL SET!U1,1,NX,1,NY,GTEMP0,NY,NX) MEL09620 

DO 106 1=1, NY MEL09630 

IFCMSLAB) GOTO 107 MEL09690 

JC1=C1H MEL 09650 

JC2=C2H MEL 0966 0 

JC3=C3H MEL 0967 0 

GOTO 108 MEL 09680 

107 JC1=C1 MEL 0969 0 

JC3=C3 MEL09700 

JC2 = C2 MEL 0971 0 

108 CALL GET! JC1 , GC1 , NY, NX) MEL09720 

CALL GET! JC2, GC2, NY, NX) MEL09730 

CALL GET! JC3, GC3 , NY, NX) MEL09790 

IF! I .GT.JL) GOTO 109 MEL09750 

GC1 ! 1 , 1 ) =1 . 0 MEL09760 

GC2! 1 , 1 ) =0 . 0 MEL 09770 

GC3! 1 , 1 ) =0 . 0 MEL 09780 

GOTO 106 MEL09790 

109 IF ! GC1 ! I , 1 ) . GT.1.0)GC1!I,1)=1.0 MEL09800 

IF ! GC2! 1,1) . GT.1.0)GC2!I,1)=1.0 MEL09810 

IF ! GC3! 1,1) .GT.1.0)GC3!I,1)=1.0 MEL09820 

IF !GC1!I,1) .LT. 0 . 0 )GC1! I , 1 ) =0 . 0 MEL09830 

IF ! GC2! 1,1) .LT. 0 . 0 )GC2! I , 1 ) =0 . 0 MEL09890 

IF !GC3!I,1) .LT. 0 . 0 )GC3! 1 , 1 ) = 0 . 0 MEL09850 

106 CONTINUE MEL09860 

CALL SET! JC1 , 1 , NX, 1 , NY, GC1 , NY, NX) MEL09870 

CALL SET!JC2,1,NX,1,NY,GC2,NY,NX) MEL09880 
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CALL SET(JC3,1,NX,1,NY,GC3,NY,NX) MEL09890 

RETURN MEL09900 

MEL09910 

CHAPTER 8* CALLED AT THE END OF EACH SWEEP; MEL09920 

NOT ACCESSED IF PARABOLIC. MEL09930 

MEL 09990 

iOO CONTINUE MEL 09950 

RETURN MEL09960 

MEL 0997 0 

CHAPTER 9= CALLED AT THE END OF EACH TIME STEP; MEL09980 

NOT ACCESSED IF PARABOLIC. MEL09990 

MEL10000 

100 CONTINUE MEL10010 

RETURN MEL10020 

MEL10030 

CHAPTER 1 0 : SET PHASE 1 DENSITY HERE WHEN IRH01=-1 IN DATA. MEL10090 

SET CURRENT-Z 'SLAB* DENSITY, Dl, IF MSLAB= . T . , MEL10050 

EG. IF(MSLAB) CALL SET( Dl , 1 , NX , 1 , NY, GDI , NY, NX) . MEL10060 

SET NEXT LARGER-Z ’SLAB* DENSITY, D1H, IF HSLAB= . T . & PARAB=F MEL10070 

EG. IF(HSLAB) CALL SET( D1H, 1 , NX, 1 , NY, GD1H, NY, NX) . MEL10080 

SET D(LN(D1))/DP (IE. D1DP) FOR UNSTEADY FLOW, MEL10090 

EG. IF(MSLAB) CALL SET ( Dl DP , 1 , NX, 1 , NY, GDI DP , NY, NX) . MEL10100 

MEL10110 

00 CONTINUE MEL 10120 

xxxxxxxxxxxxxCHAPTER 1 0 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxMEL10130 

IF (MSLAB) GOTO 101 MEL10190 

JD1=D1H MEL10150 

GOTO 102 MEL 1016 0 

01 JD1 = D1 MEL 1 0 17 0 

02 CALL SET( JD1,1,NX, 1 , NY, GDI , NY, NX) MEL10180 

RETURN MEL 1 0 1 9 0 

MEL 10200 

CHAPTER 11: SET PHASE 2 DENSITY HERE WHEN IRH02=-1 IN DATA. MEL10210 

SET CURRENT-Z ’SLAB’ DENSITY, D2, IF MSL AB= . T . , MEL10220 

EG. IF(MSLAB) CALL SET( D2, 1 , NX, 1 , NY, GD2 , NY, NX) . MEL10230 

SET NEXT LARGER-Z ’SLAB’ DENSITY, D2H, IF HSL AB= . T . & PARAB=F MEL10290 

EG. IF(HSLAB) CALL S ET( D2H , 1 , NX, 1 , NY, GD2H , NY, NX ) . MEL10250 

SET D(LN(D2))/DP FOR UNSTEADY FLOW, MEL10260 

EG. IF(MSLAB) CALL SET ( D2DP , 1 , NX, 1 , NY, GD2DP, NY, NX ) . MEL10270 

MEL 1 028 0 

00 CONTINUE MEL10290 

;xxxxxxxxxxx CHAPTER 11 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx MEL 10300 

IF (MSLAB) GOTO 103 MEL10310 

JD2=D2H MEL 10320 

GOTO 109 MEL10330 

03 JD2=D2 MEL10390 

09 CALL SET ( J D2 , 1 , NX, 1 , NY, GD2 , NY, NX ) MEL10350 

RETURN MEL10360 

MEL 1 037 0 

CHAPTER 12: SET PHASE 1 VISCOSITY HERE WHEN IEMU1=-1 IN DATA. MEL10380 

SET CURRENT-Z ’SLAB’ VISCOSITY (MU1), IF MSLAB= . T . , MEL10390 

EG. IF(MSLAB) CALL SET(MU1 , 1 , NX, 1 , NY, GVISC, NY, NX) . MEL10900 

SET NEXT LARGER-Z ’SLAB’ VISC. (MU1H), IF HSLAB= . T . 8 PARAB = F MEL10910 

EG. IF(HSLAB) CALL SET(MU1H, 1 , NX, 1 , NY, GVSCH, NY, NX) . MEL10920 

MEL 10930 

CHAPTER ALSO ACCESSED WHEN EMULAM=-1.0 IN DATA, SO THAT THE MEL10990 

LAMINAR VISCOSITY WHICH APPEARS IN WALL FUNCTIONS & IN THE MEL10950 

KE-EP TURBULENCE MODEL (IEMU1=2) MAY BE SET NON-CONSTANT. MEL10960 

SET CURRENT-Z ’SLAB’ VALUE (MU1LAM) WHEN LAMMU= . T . , MEL10970 

EG. IF(LAMMU) CALL SET( MU1 L AM, 1 , NX, 1 , NY, GVSCL , NY , NX) . MEL10980 

MEL10990 

!00 CONTINUE MEL10500 

exxxxxxxxxxx CHAOTER 12 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx MEL 1 051 0 

IF(MSLAB) CALL SET(MU1 , 1 , NX, 1 , NY, GVISC, NY, NX) MEL10520 

IF ( HSLAB ) CALL SET(MU1H, 1 , NX, 1 , NY, GVISC, NY, NX) MEL10530 

IF (LAMMU) CALL SET( MU1LAM, 1 , NX, 1 , NY, GVISC, NY, NX) MEL10590 

RETURN MEL10550 

MEL 1 056 0 

CHAPTER 13: SET EXCHANGE COEFFICIENT (E.C.) FOR VARIABLE MEL10570 

INDVAR WHEN SIGMA ( INDVAR ) =-l . 0 IN DATA. MEL10580 

SET CURRENT-Z ’SLAB’ E.C. (EXCO) IF MSL AB= . T . , MEL10590 

EG. IF(MSLAB) CALL SET( EXCO, 1 , NX, 1 , NY, GEXCO, NY, NX) . MEL10600 
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C SET NEXT SMALLER-Z 'SLAB' E.C. (EXCOL) IF LSLAB=.T., MEL10610 

C EG. IF(LSLAB) CALL SETC EXCOL , 1 , NX, 1 , NY, GEXCOL , NY, NX) . MEL10620 

C SET NEXT LARGER-Z ’SLAB' E.C. (EXCOH) IF HSLAB=.T., MEL10630 

C EG. IF(HSLAB) CALL SETC EXCOH, 1 , NX, 1 , NY, GEXCOH , NY, NX) . MEL10640 

C NOTE: FOR MSLAB, INDVAR=U1 , . . C4; FOR LSLAB, INDVAR = U1L , . . CAL MEL10650 

C a FOR HSLAB, INDVAR=U1H, . . C4H . IF PARAB= . T . SET MSLAB ONLY. MEL10660 

c MEL 1 067 0 

1300 CONTINUE MEL10630 

C ************ CHAPTER 13 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx MEL10690 

IFC INDVAR .NE.H1) GOTO 131 MEL10700 

IF (LSLAB) CALL SETC EXCOL , 1 , NX, 1 , NY, GGAMH1 , NY, NX) MEL10710 

IF (MSLAB) CALL SET ( EXCO , 1 , NX, 1 , NY , GGAMH1 , NY, NX ) MEL10720 

IF (HSLAB) CALL SET( EXCOH, 1 , NX, 1 , NY, GGAMH1 ,NY, NX) MEL10730 

131 IF (INDVAR .NE.C1) GOTO 132 MEL10740 

IF (LSLAB) CALL SET( EXCOL , 1 , NX, 1 , NY, GD( 1 , NY, 1 ) , NY, NX) MEL10750 

IF (MSLAB) CALL SETC EXCO, 1 , NX, 1 , NY, GDC 1 , NY, 1 ) , NY, NX) MEL10760 

IF (HSLAB) CALL SET ( EXCOH , 1 , NX, 1 , NY , GDC 1 , NY, 1 ) , NY, NX) MEL10770 

132 IFC INDVAR .NE. C2) GOTO 133 MEL10780 

IF (LSLAB) CALL SET ( EXCOL , 1 , NX, 1 , NY , GD( 2 , NY , 1 ) , NY, NX) MEL10790 

IF (MSLAB) CALL SETC EXCO , 1 , NX, 1 , NY, GDC 2, NY, 1 ) , NY, NX) MEL10800 

IF (HSLAB) CALL SET( EXCOH , 1 , NX, 1 , NY, GD( 2, NY, 1 ) , NY, NX) MEL10810 

133 IF (INDVAR .NE.C3 ) GOTO 134 MEL10820 

IF (LSLAB) CALL SETC EXCOL , 1 , NX, 1 , NY, GD( 3, NY, 1 ) , NY, NX) MEL10830 

IF (MSLAB) CALL SETC EXCO, 1 , NX, 1 , NY, GDC 3, NY, 1 ), NY, NX) MEL10840 

IF (HSLAB) CALL SET ( EXCOH , 1 , NX, 1 , NY, GD( 3 , NY , 1 ) , NY, NX) MEL10850 

134 RETURN - MEL10860 

C MEL1087 0 

C CHAPTER 14: SET INTER-PHASE FRICTION COEFFICIENT (CFP) HERE MEL10380 

C WHEN ICFIP = -1 IN DATA; ITS UNITS = FORCE / (CELL * RELATIVE MEL10390 

C SPEED OF PHASES). MEL10900 

C MEL10910 

1400 CONTINUE MEL10920 

RETURN MEL10930 

C MEL 10940 

C CHAPTER 15: SET INTER-PHASE MASS-TRANSFER RATE PER CELL (MDT ) MEL10950 

C HERE WHEN IMDOT = -1 IN DATA. MEL10960 

C MEL 1 097 0 

1500 CONTINUE MEL10980 

C ************ CHAPTER 15 *************************************** MEL10990 

CALL SET (MDT , 1 , NX, 1 , NY,GMDT , NY, NX) MEL11000 

RETURN MEL11010 

c MEL11020 

C CHAPTER 16: SET HERE PHASE 182 SATURATION ENTHALPIES MEL11030 

C ( HST1 8 HST2) WHEN IHSAT = -1 IN DATA. MEL11040 

c MEL11050 

1600 CONTINUE MEL11060 

RETURN MEL 11070 

END MEL11080 

C$DIRECTIVE**STRIDE MEL11090 

SUBROUTINE STRI DE( IZSTEP , IGOTO , IRN ) MEL11100 

c MEL11110 

C USE THIS SUBROUTINE TO SPECIFY THE GEOMETRY MEL11120 

C OF THE FORWARD STEP IN PARABOLIC CALCULATIONS. MEL11130 

C IZSTEP IS THE CURRENT FORWARD STEP, a NZSTP IS THE LAST MEL11140 

C FORWARD STEP (FOR PARAB= . T . EARTH SETS NZ=1 ). MEL11150 

C THE COMMON VARIABLE 'ZWL 1 GIVES THE DISTANCE OF THE MEL11160 

C PREVIOUS STEP FROM THE ORIGIN. MEL11170 

c MEL11180 

C$INCLUDE 9 , CMNGUSSI . FTN/G MEL11190 

LOGICAL LOGIC1, LOGIC MEL11200 

DIMENSION LOGIC(IOO) MEL11210 

COMMON/ L DATA/ LOGI Cl ( 309 ) MEL 11220 

EQUIVALENCE (LOGIC(l) , LOGICK210) ) MEL11230 

DIMENSION INTGR(IOO) MEL11240 

COMMON/ IDATA/INTGRK 194) MEL 11250 

EQUIVALENCE ( INTGR( 1 ) , INTGR1 ( 95) ) MEL11260 

DIMENSION RE(100) MEL11270 

COMMON/ RDAT A/ RE1 (421 ) MEL11280 

EQUIVALENCE ( RE( 1 ) , RE1 ( 322 ) ) MEL11290 

COMMON/BOUND/LOCREG(60), MEL 11 30 0 

aTRl,CPlRl(7),VPlRl(7),CP2Rl(5),VP2Rl(5),CPNRl(5),VPNRl(5), MEL 11310 

8TR2,CP1R2(7),VP1R2(7),CP2R2(5),VP2R2(5),CPNR2(5),VPNR2(5), MEL 11 320 
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aTR3,CPlR3(7),VPlR3(7),CP2R3(5), VP2R3( 5 ) , CPNR3( 5 ) , VPNR3( 5 ) , 
8TRA,CP1RA(7), VP1RA(7 ) , CP2RA( 5 ) , VP2RA( 5 ) , CPNRA( 5 ) , VPNRA( 5 ) , 
aTR5,CPlR5(7),VPlR5(7),CP2R5(5),VP2R5(5),CPNR5(5),VPNR5(5), 
aTR6,CPlR6(7),VPlR6(7),CP2R6(5),VP2R6(5),CPNR6(5),VPNR6(5), 

3TR7 , CP1R7 ( 7 ) , VP1R7 ( 7 ) , CP2R7 ( 5 ) , VP2R7 ( 5 ) , CPNR7 ( 5 ) , VPNR7 ( 5 ) , 
aTR8,CPlR8(7),VPlR8(7),CP2R8(5),VP2R8(5) , CPNR8 ( 5 ) , VPNR8 ( 5 ) , 
aTR9,CPlR9(7),VPlR9(7),CP2R9(5),VP2R9(5),CPNR9(5),VPNR9(5), 
aTR10,CPlR10(7),VPlR10(7),CP2R10(5),VP2R10(5),CPNR10(5),VPNR10(5) 
DIMENSION TCVREG( 350 ) 

EQUIVALENCE (TCVREGC 1) , TR1 ) 

NCLUDE 9, GUSSEQUI . FTN/G 

LOGICAL CARTES, POL AR,SPDATA, SKEW, TWODYZ, ONEDZ, STOVAR( 25), 

& S0LVAR(25),PRINT(25),RESID(25), CMPRSS , CON EMU, LSP1, 

& CONRHO , EMDOT , ONEPHS , INCORE( 1 0 ) , I SAVED , SAVEI , SAVEM, 

& RESTRT,XCYCLE, MONITR, REGION ( 10), STEADY, WHOL EP , SL ABPP , 

& RAIN, BLOCKZ, PWSTAG ,RUN(30),PLOT, RESMAP , FLAG, BLOCK, 

& TEST, CONCICA) , DISTIL 
LOGICAL SPRESS , PARAB , DONACC , OVERLY , SACC, 

& GUSSIE, CATHY, CONNIE, CORA, ESTER, FLASH, FLORA, 

& FOCS, GENMIX, HESTER, PICALO, PL ANT, SPLASH, HELGA, 

& TACT,TIBALT,TOPSI,PAMELA,TABLES,WSTAG, 

& CONMOD,GROSTA,SUBPST, SUBWGR 
EQUIVALENCE (CARTES, LOGICl(l) ) , (POLAR, LOG IC1(2)), 

& (SPDATA,L0GIC1(3) ), ( SKEW , LOGICl ( A )),( TWODYZ, LOGICl ( 5 )) , 

& (0NEDZ,L0GIC1(6)),(ST0VAR(1),L0GIC1(7)), 

8 ( SOL VAR(1),L0GIC1( 32)), (PRINT (1) , LOGICl ( 57 ) ) , 

a ( RESI D(l), LOGICl ( 32)), (CMPRSS, LOGICl (107)), 
a (CONEMU, LOGICl (108) ), (LSP1, LOGICl (109) ), 
a ( CONRHO, LOGICl ( 110)), (EMDOT, L OGIC1 ( 1 1 1 ) ) , 
a ( ONEPHS, LOGICl (1 12) ), (INCORE(l) , LOGICl (113)) 

EQUIVALENCE( I SAVED, LOGICl (123) ), (SAVEI, L OGIC1 ( 12A ) ) , 
a (SAVEM,L0GIC1(125)),(RESTRT,LGGIC1(126)), 
a (XCYCLE, LOGICl (127)) , (MONITR, LOGICl ( 128 )) , 
a (REGION(l) , LOGICl ( 129 ) ) ,( STEADY , LOGICl ( 1 39 )) , 
a (WHOL EP, LOGICl (1 AO) ) , ( SL A3P° , L 0GIC1 ( 1 A1 ) ) , 
a (RAIN, LOGICl (1A2)), ( BLOCKZ, L OGIC1 ( 1 A3 ) ) , 
a (PWSTAG, LOGICl(lAA) ), (RUN(l) , LOGICl (1A5) ), 
a ( PLOT, LOGICl (17 5 ) ), (RESMAP, LOGICl ( 176) ), 

a (FLAG, LOGICl (177)), (BLOCK, LOGICl ( 178 )), (TEST, LOGICl ( 179 ) ) 
EQUIVALENCE(CONCKl), LOGICl (206) ) 

EQUIVALENCE ( SPRESS , LOGIC( 100 )),( PARAB, LOGIC( 99 )) , 

8 (DONACC, L0GIC(98)) 

EQUIVALENCE ( OVERLY, LOGIC( 96 )),( SACC, LOGIC( 95) ) , 
a (WSTAG, L0GIC(9A) ) 

EQUIVALENCE ( DISTI L , LOGIC( 90 ) ) 

EQUIVALENCE (PAMELA, LOGIC (67 ) ) , (TOPSI , LOGIC (66 )) , 
a (TIBALT,L0GIC(65)),(TACT,L0GIC(6A)),(HELGA,L0GIC(63)), 
a (SPLASH, L0GIC(62)), (PLANT, L0GIC(61)), (PICALO, LOGIC(60) ), 
a (HESTER, L0GIC(59)), (GENMIX, LOGIC (58) ), (FOCS, L0GIC(57) ), 
a (FLORA, LOGIC (56 )), (FLASH, L0GIC(55) ), (ESTER, L0GIC(5A)) , 
a (C0RA,L0GIC(53) ), (CONNIE, L0GIC( 52)), (CATHY, LOGIC( 51)), 
a (GUSSIE, L0GIC(50)) 

EQUIVALENCE (TABLES, L0GIC(A8)),( CONMOD, L0GIC( A7 ) ) , 
a (GROSTA,LOGIC(A6)),(SUBPST,LOGIC(A5)) 

EQUIVALENCE (SUBWGR, LOGIC( A3) ) 

INTEGER FSTEP,FSWEEP,TSTSWP,ITAB(8),MTABVR(8) 

DIMENSION ISPCS0(25),LITER(25) 

EQUIVALENCE (NX, INTGR1 ( 1 ) ) , ( NY, INTGR1 ( 2) ) , 
a (NZ,INTGR1(3)),(ISPCS0(1),1NTGR1(A)), 
a ( NREGN, INTGR1 ( 29 ) ) , ( NPHI , INTGR1 ( 30 ) ) , 
a (LITKE, INTGRK31)) , (LITHYD, INTGRK32) ) , 
a (LITH,INTGR1(33)), ( LITCNC, INTGR1 ( 3 A ) ) , 

a (LITSLB,INTGR1(35)), ( NRUN, INTGR1 ( 36 ) ) , 

8 (LITER(1),INTGR1(37)),(FSTEP,INTGR1(62)), 

a (FSWEEP,INTGR1(6 3) ), (LSTEP, INTGRK6A)) 

EQUIVALENCE ( L SWEEP, INTGR1 ( 65 ) ) , (NPRINT, INTGR1 (66 ) ) , 
a (IERRP, INTGRK67)), (IMAXP, INTGRK68) ), 
a ( IEMU1 , INTGRK69)), (IXMON, INTGR1 (70 ) ) , 

a ( IYMON, INTGRK71)), ( IZMON ,INTGR1(72)), 

a ( KEMU >INTGR1(73)),( KMAIN , INTGRK7A)), 

a ( KINDEX, INTGR1 (7 5 ) ) , ( KGEOM, INTGRK76)), 

8 (KINPUT,INTGR1(77)),(KS0DAT,INTGR1(78)) 



MEL11330 
MEL113A0 
MEL11350 
MEL 1136 0 
MEL11370 
MEL 11380 
MEL11390 
MEL11A00 
MEL11A10 
MEL11A20 
MEL11A30 
MEL11AA0 
MEL11A50 
MEL11A60 
MEL 1 1 A7 0 
MEL 1 1 A3 0 
MEL11A90 
MEL 11500 
MEL 11510 
MEL 11520 
MEL 11530 
MEL115A0 
MEL 11550 
MEL 1 1 56 0 
MEL 11570 
MEL11580 
MEL 11590 
MEL 116 00 
MEL 116 10 
MEL 11620 
MEL11630 
MEL116A0 
MEL11650 
MEL11660 
MEL 1167 0 
MEL 1168 0 
MEL 11690 
MEL11700 
MEL11710 
MEL 11720 
MEL11730 
MEL117A0 
MEL11750 
MEL 1176 0 
MEL 1177 0 
MEL 1 178 0 
MEL11790 
MEL11800 
MEL11810 
MEL 11820 
MEL11830 
MEL118A0 
MEL11850 
MEL11860 
MEL 1187 0 
MEL 1188 0 
MEL11890 
MEL 119 00 
MEL11910 
MEL 11920 
MEL11930 
MEL119A0 
MEL11950 
MEL11960 
MEL11970 
MEL11980 
MEL 11990 
MEL 12000 
MEL 1201 0 
MEL 12020 
MEL12030 
MEL120A0 
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a 
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a 

a 

a 

a 

a 

a 



EQUIVALENCECKFL AG, INTGRH 79) ), CKCOMPF, INTGR1C80) ) , 

CKSORCE, INTGRK81)), ( KL ES1D, INTGRH 82 ) ), 

CKLES2D, INTGRH 83)), ( KL ES3D, INTGR1 C 89 ) ), 
CKC0MPP,INrGRlC85)),CKADJST,INTGRlC86)), 

CKFLUX, INTGRH 87)), C KSHI FT , I NTGR1 C 88 ) ) , 

(KOUTPT, INTGRH 89 )),CKDIF, INTGRH 90)), 

(KCOMPU, INTGRH 91) ),(KCOMPV, INTGRH 92)), 

( KCOMPW, I NTGR1 ( 93 ) ) , ( KCOMPR, I NTGRK 9 A ) ) 

EQUIVALENCE ( KEL I N , I NTGR( 10 0 ) ) , ( MEANDF, I NTGR( 99 ) ) , 

( NUMCLS , INTGRC98)), 

(IRH01,INTGR(97 )), (IRH02,INTGR(96) ), 
CIZW1,INTGRC95)),CIZW2,INTGRC99)), 

(MGRI D, INTGR(93)),( KWALL , INTGRC 92) ) , 
(IZPR1,INTGR(91)),(IZPR2,INTGRC90)), 

(I STPR1, INTGRC 89) ), (ISTPR2, INTGRC 88) ), 

C NT PR IN, INTGRC 87 ) ) 

EQUIVALENCE C IMDOT , INTGRC 86 )), C IHSAT , INTGRC85) ) , 

CICFIP, INTGRC 84) ), ( NTABL E, I NTGRC 83 ) ) , C NTABVR , INTGRC 82 ) ) , 

C L INTAB, INTGRC 81 )),CNPRTAB, INTGRC 80) ), CNMON, INTGRC 79)), 

C ITABC1), INTGRC 71 )), ( MTABVR C 1 ) , I NTGRC 6 3 ) ) , 
CNZPRIN,INTGRC62)), 

CNXPRIN, INTGRC61) ) , CNYPRIN, INTGRC 60) ), 

CLDISTL, INTGRC 59)) , CTSTSWP, INTGRC 58) ), 

C KDB EXP, INTGRC 57 )), ( KDBRHO, INTGRC56 ) ) , 

C KDBMDT, INTGRC 55) ) 

EQUIVALENCE C I LOOP1 , INTGRC 53) ) , C ILOOPN, I NTGRC 52) ) , 
CIPBP,INTGRC51)),CIZPRF,INTGRC50)), 

CIZPRL, INTGRC 99)), C ISTPRF, INTGRC 98 ) ) , 

C ISTPRL , INTGRC 97 ) ) , C KDBGEN , INTGRC 96 ) ) , 

C I VELF, INTGRC 95 )), Cl VELL, INTGRC 99)), 
CIKEF,INTGRC93)),CIKEL,INTGRC92)) , 

CIENTF, INTGRC91)), Cl ENTL , INTGRC 90 )) , 

C ICNCF, I NTGRC 39 )), (ICNCL, INTGRC 38) ), 

C NV EL, INTGRC 37 )),CNKE,I NTGRC 36)), 

CNENT, INTGRC 35)), C NCNC, INTGRC 39 ) ) , 

C NZSTP, INTGRC 33 )), CNPRMNT, INTGRC 32)) 

DIMENSION SPARE1 C 20),XUDISTC30),YVDIST(30), ZWDIST C 50 ) , 
SIGMAC25),CRITC25),DTFALSC25),RESREFC25), 
TITLE(25),FIINITC25), TFRACC 30 ) 

EQUIVALENCE C SPARE1 C 1 ) , RE1 C 1 ) ) , 

C TFRACC 1 ), RE1 C 21 ) ) , C XULAST, RE1C 51)), 

C YV LAST, RE1C52)),CZW LAST, RE1C53)), 
CXUDISTC1),RE1C59)),CYVDISTC1),RE1C89)), 
CZWDISTC1),RE1C119)), C SIGMAC 1 ) , RE1 C 169 ) ) , 
CCRITC1),RE1(189)),(DTFALSC1),RE1C219)), 

CRESREFC1) ,RE1C239) ), C EMU1 , RE1 C 269 ) ) , 
CRH01,RE1C265)),CRH02,RE1C266)), 

C TLAST , RE1 C 267 ) ) , C CFIPS , RE1C 268 ) ) 

EQUIVALENCE C AMDOT , RE1 C 269 ) ) , C FI INI TC 1 ) , RE1 C 270 ) ) , 
CRELAXP,RE1C295)),(TITLEC1),RE1(296)), 

C DT , RE1C321 ) ) , 

CRI NNER,ZWDISTC50)),CSNALFA, ZWDIST C 99)), 

CPBAR,ZWDIST(98) ) 

EQUIVALENCE C SLOEMU , REC 10 0 ) ) , C SLORHO , REC 99 ) ) , C RLXRHO , REC 98 ) ) 
CRHOMAX,REC97 )), (RHOMIN,REC96) ), 

C EMUMAX, REC 95 ) ) , C EMUMIN, REC 99 ) ) , 

C TKEMAX, REC93)),C TKEMIN, REC 92 ) ) , 

C RLXPZ, REC 9 1 ) ) , 

C ABUOY , REC 90 ) ) , C HREF, REC 89 ) ) , 

C AGRAVX, REC 88 ) ) , C AGRAVY, REC 87 ) ) , 

C AGRAVZ, REC86)),CATIME,REC85)) 

EQUIVALENCE C RLXPXY, REC 89 ) ) , 

C EPSMAX,REC83)),CEPSMIN,REC82)), 
CARH01,REC81)),CBRH01,REC80)),CCRH01,REC79)) 

EQUIVALENCE C AZW2, REC 73 ) ) , C BZW2 , REC 79 ) ) , 

CPINT , REC7 5 ) ) , CCZW2, REC76 ) ) , 

CDZ,REC78) ) , C ZWL , REC77) ) , CZW2M1T, REC72) ) , 

CVELMIN,REC70)), CVELMAX,REC69) ), 

CFALC0M,REC68)), CCMDOT,REC67)), 

CRLXMDT , REC66 ) ) , CAMDTMX, REC 6 5 ) ) , C AMDTMN , REC 6 9) ) , 

C RADMAX, REC63)),C RADMIN, REC 62) ) , 

( PRESS, REC61) ), CCP1,REC60)), 



MEL 12050 
MEL 1206 0 
MEL 12070 
MEL 12080 
MEL 12090 
MEL 12100 
MEL12110 
MEL 12120 
MEL12130 
MEL 12190 
MEL12150 
MEL12160 
MEL 12170 
MEL12180 
MEL12190 
MEL12200 
MEL12210 
MEL12220 
MEL12230 
MEL12290 
MEL12250 
MEL 1226 0 
MEL 1227 0 
MEL 1228 0 
MEL12290 
MEL12300 
MEL12310 
MEL12320 
MEL12330 
MEL 12390 
MEL 12350 
MEL 1236 0 
MEL12370 
MEL12380 
MEL12390 
MEL 1290 0 
MEL 1291 0 
MEL 12920 
MEL 12930 
MEL 12990 
MEL 12950 
MEL12960 
MEL 1297 0 
MEL12980 
MEL12990 
MEL12500 
MEL12510 
MEL12520 
MEL12530 
MEL12590 
MEL 12550 
MEL12560 
MEL1257 0 
MEL12580 
MEL12590 
MEL12600 
MEL12610 
MEL12620 
MEL 12630 
MEL12690 
MEL12650 
MEL12660 
MEL12670 
MEL12680 
MEL12690 
MEL 127 0 0 
MEL12710 
MEL12720 
MEL127 30 
MEL 127 90 
MEL12750 
MEL 127 6 0 
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& (CP2,RE(59) ), (CP3,RE(58)), ( FSTOIC, REC 57 ) ) , 

& CARRC0N,RE(56)),(PREEXP,RE(55)), CCA1 I , REC 59 ) ) 
EQUIVALENCE (CA2I , RE( 53 ) ) , (CA13, RE< 52 ) ) , (CA23 , RE(51 ) ) , 
8 (H1SAT,RE(50)),(H2SAT,RE(99)), 

8 ( EMULAM, RE( 48 ) ) 

EQUIVALENCE ( TWPRCN, RE( 42) ) , C PSATEX, RE( 91 ) ) 

EQUIVAL ENCEC EL EXP , REC 39 ) ) , ( EWAL L , REC 38 ) ) , 

8 (SIGE, REC 37 ) ) , ( S IGK , REC 36 ) ) , ( TAUDK , REC 35 ) ) 

DIMENSION L0C1(6),L0C2(6) , L0C3C 6 ) , L0C9C 6 ) , L0C5C6 ) , 

8 LOC6C6 ), LOC7C6), L0C8C6) , L0C9C6), L0C10C6) 

EQUIVALENCE ( LOCREGC 1 ), LOCI ( 1 )),( LOCREGC 7 ), L0C2C 1 ) ) 

( LOCREG( 13),L0C3(1)),( LOCREGC 19), LOCAL 1) ) 

( LOCREGC 25) , L0C5( 1 ) ) , ( LOCREG( 31 ) , LOC6 ( 1 ) ) 

( L OCREG( 37 ) , LOC7 ( 1 ) ) , ( LOCREGC 93 ) , L0C8 ( 1 ) ) 
(LOCREGC 49), L0C9C 1 ) ) , ( L0CREGC55) , LOCIOC 1 ) ) 
XUFRAC(30),YVFRAC(30) .ZWFRACC30), ENDIT(25) 
XFRACC30) ,YFRACC30) ,2FRAC(30) 



EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
DIMENSION 
DIMENSION 



EQUIVAL ENCEC XUDISTC 1 ) , XUFRACC 1 ) , XFRAC( 1 ) ) , 

8 (YVDIST(1),YVFRAC(1),YFRAC(1) ), 

8 ( ZWDISTC 1 ) , ZWFRACC 1 ) , ZFRACC 1 ) ) , ( LITCNC, LITC) , 

8 ( ENDITC 1 ) , CRITC 1 ) ) , ( RLXP , RELAXP ) 

EQUIVALENCE (NPRMON, LITER CD) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 1 STARTS: 
SATLIT-EQUIVALENT IRUN; 

EQUIVALENCE ( IRUN, INTGRC 11 ) ) 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

USER PLACES HIS VARIABLES, ARRAYS, 

USER PLACES HIS DATA STATEMENTS HERE. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

GO TO (10,20,30,90,50,60), IGOTO 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



STANDARD SECTION 1 ENDS. 
USER SECTION 1 STARTS: 
EQUIVALENCES ETC. HERE. 

SECTION 1 ENDS. 
STANDARD SECTION 2 STARTS; 



STANDARD SECTION 2 ENDS. 
USER SECTION 2 STARTS: 



SECTION 1: SET FORWARD STEP SIZE DZ FOR IZSTEP>1 WHEN 
ZWL AST <0. IN DATA. 

AT IZSTEP=1 EARTH SETS DZ = ABSCZWDISTC 1 ) *ZWLAST) 



10 CONTINUE 

SXXXUSER SETS DZ HERE. . . 
RETURN 



SECTION 2: SET X-WIDTH (XULAST) OF GRID FOR IZSTEP > 1 
WHEN XULAST < 0 . IN DATA. 

AT IZSTEP=1 EARTH SETS XULAST = ABSCXULAST) 



20 CONTINUE 

EXXXUSER SETS XULAST HERE... 
RETURN 



SECTION 3: SET Y-WIDTH (YVLAST ) OF GRID FOR IZSTEP > 1 
WHEN YVLAST < 0 IN DATA. 

AT IZSTEP=1 EARTH SETS YVLAST = ABSCYVLAST) 



30 CONTINUE 

EXXXUSER SETS YVLAST HERE... 
RETURN 



SECTION 9: SET INNER RADIUS (RINNER) OF GRID FOR IZSTEP > 1 
WHEN RINNER < 0 IN DATA. 

AT IZSTEP=1 EARTH SETS RINNER = ABSCRINNER) 



90 CONTINUE 

EXXXUSER SETS RINNER HERE... 
RETURN 



SECTION 5: SET SLOPE CSNALFA) OF INNER EDGE OF GRID FOR 
IZSTEP > 1 WHEN SNAL FA> 1, & CARTES = .FALSE. 

AT IZSTEP = 1 EARTH SETS SNALFA = SNAL FA-2 . 



50 CONTINUE 



MEL 1277 0 
MEL12780 
MEL 12790 
MEL12800 
MEL 1281 0 
MEL 1 2820 
MEL12830 
MEL 12890 
MEL 12850 
MEL 1286 0 
MEL 1287 0 
MEL12880 
MEL12890 
MEL12900 
MEL 1291 0 
MEL 12920 
MEL 12930 
MEL 12990 
MEL 12950 
MEL 1 296 0 
MEL 1297 0 
MEL12980 
MEL 12990 
MEL 13000 
MEL13010 
MEL 13020 
MEL13030 
MEL 13090 
MEL 13050 
MEL13060 
MEL 1307 0 
MEL 1 308 0 
MEL13090 
MEL 13100 
MEL 131 1 0 
MEL13120 
MEL13130 
MEL13190 
MEL 13150 
MEL13160 
MEL 1317 0 
MEL13180 
MEL13190 
MEL13200 
MEL13210 
MEL 13220 
MEL 13230 
MEL 13290 
MEL 13250 
MEL 1326 0 
MEL 1327 0 
MEL 1328 0 
MEL13290 
MEL 13300 
MEL13310 
MEL 1332 0 
MEL 13330 
MEL 1 3390 
MEL13350 
MEL13360 
MEL 1337 0 
MEL13380 
MEL 1 3390 
MEL 1 3900 
MEL13910 
MEL13920 
MEL13930 
MEL13990 
MEL13950 
MEL 1 396 0 
MEL1397 0 
MEL13980 
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CXXXXXUSER SETS SNALFA HERE... MEL13490 

RETURN MEL13500 

MEL13510 

SECTION 6: SET MEAN PRESSURE ( PBAR) AT NEXT FORWARD STEP MEL13520 

WHEN PBAR < 0. IN DATA. MEL13530 

FOR UNCONFINED FLOWS WITH IMPRESSED NON-ZERO MEL13540 

PRESSURE GRADIENTS SET PBAR HERE; FOR CONFINED MEL13550 

FLOWS EARTH AUTOMATICALLY COMPUTES PRESSURE REQUIRED. MEL13560 

MEL 1357 0 

60 CONTINUE MEL13580 

CXXXXXUSER SETS PBAR HERE... MEL13590 

RETURN MEL 136 00 

END MEL136 10 



34 



Es MEL3SAT FORTRAN A1 



[ RECTI VEX*S AT LIT 

*FILE NAME: MODSTL.FTN 

^ABSTRACT: SATELLITE MODEL MAIN PROGRAM. 

^DOCUMENTATION: PHOENICS INSTRUCTION MANUAL (SPRING 1983). 
^AUXILIARY SUBROUTINES (TAPES, ETC.) ARE IN SATELLITE LIBRARY 
SERVICEU, WHICH MUST BE INCLUDED IN LINK EDIT TO RUN. 
(XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 1 STARTS: 



5 TER 1 COMMON BLOCKS AND USER'S DATA. 



‘JCLUDE 9 , CMNGUSSI . FTN/G 
LOGICAL LOGIC1, LOGIC 
DIMENSION LOGIC(IOO) 

COMMON/ L DATA/ L OG IC 1 ( 309 ) 

EQUIVALENCE ( LOGIC( 1 ) , L0GIC1 (210 ) ) 

DIMENSION INTGR(IOO) 

COMMON/IDATA/INTGR1 (199) 

EQUIVALENCE ( INTGR( 1 ) , INTGR1 ( 95) ) 

DIMENSION RE(100) 

COMMON/RDATA/REK 421 ) 

EQUIVALENCE ( RE( 1 ) , RE1 ( 322) ) 

COMMON/ BOUND/ L OCREG ( 6 0 ) , 

&TR1,CP1R1(7),VP1R1(7),CP2R1(5) , VP2R1 ( 5) , CPNR1 ( 5) , VPNR1 ( 5) , 
8TR2,CP1R2(7 ) ,VP1R2(7),CP2R2(5) ,VP2R2(5) ,CPNR2(5) ,VPNR2(5) , 
8TR3,CP1R3(7), VP1R3(7 ) , CP2R3 ( 5 ) , VP2R3( 5 ) , CPNR3( 5 ) , VPNR3( 5 ) , 
&TR4,CP1R4(7), VP1R4(7) , CP2R4( 5 ) , VP2R4( 5 ) ,CPNR4( 5) , VPNR4(5 ) , 
8TR5,CP1R5(7), VP1R5 ( 7 ) , CP2R5( 5 ) , VP2R5( 5 ) , CPNR5( 5 ) ,VPNR5(5), 
&TR6,CP1R6(7),VP1R6(7),CP2R6(5),VP2R6(5),CPNR6(5),VPNR6(5), 
8TR7,CP1R7(7),VP1R7(7),CP2R7(5),VP2R7(5),CPNR7(5),VPNR7(5), 
&TR8,CP1R8(7),VP1R3(7),CP2R8(5),VP2R8(5),CPNR8(5),VPNR8(5), 
8TR9,CP1R9(7),VP1R9(7),CP2R9(5),VP2R9(5),CPNR9(5),VPNR9(5), 
8TR10,CP1R10(7) ,VP1R10(7) ,CP2R10(5), VP2R10(5) ,CPNR10(5),VPNR10(5) 
DIMENSION TCVREG( 350 ) 

EQUIVALENCE (TCVREG( 1 ) , TR1 ) 

ICLUDE 9,GUSSEQUI . FTN/G 

LOGICAL CARTES, POL AR,SPDATA, SKEW, TWODYZ, ONEDZ, STOVAR( 25 ) , 

& S0LVAR(25),PRINT(25),RESID(25), CMPRSS , CON EMU, LSP1 , 

& CONRHO, EMDOT, ONEPHS, INCORE( 10 ) , I SAVED , SAVEI , SAVEM, 

8 RESTRT , XCYCL E, MONITR, REGI ON( 1 0 ) , STEADY, WHOL EP, SLABPP, 

8 RAIN, BLOCKZ, PWSTAG, RUN(30 ), PLOT, RESMAP, FLAG, BLOCK, 

a TEST, CONCK4) , DISTIL 
LOGICAL SPRESS , PARAB , DON ACC, OVERLY , SACC, 
a GUSSIE, CATHY, CONNIE, CORA, ESTER, FLASH, FLORA, 
a FOCS, GENMIX, HESTER, PI CAL 0, PL ANT, SPL ASH, HELGA, 
a T ACT, TIBALT,TOPSI, PAMELA, TABLES, WSTAG, 
a CONMOD,GROSTA, SUBPST, SUBWGR 
EQUIVALENCE (CARTES, LOGICl(l) ), (POLAR, L0GICK2) ) , 
a (SPDATA,L0GIC1(3)),(SKEW,L0GIC1(4)),( TWODYZ, L0GICK5)), 
a (0NEDZ,L0GIC1(6)),(ST0VAR(1),L0GIC1(7)), 
a (S0LVAR(1),L0GIC1(32)),(PRINT(1),L0GIC1(57)), 
a ( RESI D(1 ), LOGICK 82)), (CMPRSS, L0GICK1 07)), 

8 ( C0NEMU,L0GIC1(1 08 )), (LSP1, LOGICK 109)) , 

8 (CONRHO, LOGIC1(110)),( EM DOT, LOGICl(lll)), 

a ( ONEPHS, LOGICK 1 12) ), (INCORE(l), LOGICK 113) ) 

EQUIVALENCE( I SAVED, LOGICK 123)), (SAVEI, LOGICK 124)), 
a ( SAVEM, LOGICK 1 25 )), (RESTRT, LOGICK 1 26)), 

a ( XCYCL E, LOGICK 127)), (MONITR, LOGICK 128)), 
a ( REGION (1), LOGICK 129) ), (STEADY, LOGICK 139) ), 

a ( WHOL EP, LOGICK 140)), (SLABPP, LOGICK 141)), 
a (RAIN, LOGICK 142)), (BLOCKZ, LOGICK 143)), 
a ( PWSTAG, LOGICK 144) ),(RUN(1), LOGICK 145)), 
a (PLOT, LOGICK 175 )), (RESMAP, LOGICK 176) ), 

a (FLAG, LOGICK 177)), (BLOCK, LOGICK 178)), (TEST, LOGICK 179)) 
EQUIVALENCE(C0NC1( 1 ), LOGICK 206 ) ) 

EQUIVALENCE ( SPRESS , LOGIC( 1 00 )),( PARAB , LOGIC( 99 )) , 

8 ( DON ACC, L0GIC(98)) 

EQUIVALENCE ( OVERLY, LOGIC( 96 )),( SACC, LOGIC( 95 )) , 

8 (WSTAG, L0GIC(94)) 

EQUIVALENCE ( DISTI L , LOGIC( 90 ) ) 

EQUIVALENCE ( PAMELA, LOGIC( 67 )),( TOPSI , LOGIC( 66 )) , 

8 (TIBALT, LOGIC( 65)), (TACT, L0GIC(64)), (HELGA, LOGIC(63)), 
a ( SPLASH, LOGIC( 62)), (PLANT, LOGIC( 61)), (PICALO, LOGIC(60 ) ) , 



MEL00010 
MEL00020 
MEL00030 
MEL00040 
MEL 00050 
MEL 0 0 060 
MEL 00 07 0 
MEL 00080 
MEL00090 
MEL00100 
MEL 00110 
MEL00120 
MEL00130 
MEL00140 
MEL00150 
MEL 0 016 0 
MEL 0017 0 
MEL 0013 0 
MEL00190 
MEL 00200 
MEL00210 
MEL 00220 
MEL00230 
MEL00240 
MEL00250 
MEL00260 
MEL 0027 0 
MEL00230 
MEL00290 
MEL 00300 
MEL00310 
MEL00320 
MEL00330 
MEL00340 
MEL00350 
MEL00360 
MEL00370 
MEL00380 
MEL00390 
MEL00400 
MEL0041 0 
MEL00420 
MEL00430 
MEL00440 
MEL00450 
MEL00460 
MEL00470 
MEL00480 
MEL00490 
MEL00500 
MEL00510 
MEL00520 
MEL 00530 
MEL00540 
MEL00550 
MEL00560 
MEL 0057 0 
MEL 00580 
MEL 0059 0 
MEL00600 
MEL 0 0610 
MEL00620 
MEL00630 
MEL00640 
MEL 00650 
MEL 0066 0 
MEL00670 
MEL00680 
MEL00690 
MEL 0 0700 
MEL00710 
MEL00720 
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a (HESTER, L0GIC(59)), ( GENMIX, LOGIC( 58 ) ) , ( FOCS, LOGIC( 57 ) ) , 

a ( FLORA, LOGIC( 56)), (FLASH, L0GIC(55)),( ESTER, L0GIC(54)), 
a ( CORA, LOGIC(53 )), (CONNI E,LOGIC( 52)), (CATHY, LOGIC( 51)), 

a (GUSSIE,LOGIC(50)) 

EQUIVALENCE (TABLES, L OGIC( 48 )),( CONMOD, LOGIC( 47 ) ) , 
a (GR0STA,L0GIC(46)),(SUBPST,L0GIC(45) ) 

EQUIVALENCE (SUBWGR, LOGIC( A3) ) 

INTEGER FSTEP,FSWEEP,TSTSWP, ITAB ( 8 ) , MTABVR( 8 ) 

DIMENSION ISPCS0(25),LITER(25) 

EQUIVALENCE (NX, INTGR1 ( 1 ) ) , ( NY, INTGR1 ( 2 ) ) , 
a (NZ,INTGR1(3) ), (ISPCSO(l), INTGRK A) ), 
a (NREGN, INTGRK29) ) , ( NPHI , INTGR1 ( 30 ) ) , 
a (LITKE, INTGRK31)) , (LITHYD, INTGR1 ( 32 ) ) , 

a (LITH,INTGR1(33)),(LITCNC,INTGR1(34)), 
a (LITSLB, INTGRK35) ) , (NRUN, INTGRK36) ) , 
a ( LITER( 1 ) , INTGR1 ( 37 ) ) , ( FSTEP , INTGR1 (62) ) , 

a (FSHEEP, INTGR1 (63)) , ( L STEP , INTGR1 ( 6 A ) ) 

EQUIVALENCE (L SHEEP, INTGR1 (65) ) , ( NPRINT , INTGR1 (66 ) ) , 
a (IERRP, INTGRK67) ), (IMAXP, INTGRK68)), 
a (IEMU1, INTGRK69) ), (IXMON, INTGRI (70 ) ) , 
a (IYMON, INTGRI (7 l)),(IZMON, INTGRI (72)), 
a (KEMU, INTGRI (73)), (KMAIN, INTGRI (79) ), 
a (KINDEX, INTGRI ( 75 )), (KGEOM, INTGRI (76)), 

a ( KINPUT, INTGRI ( 77 )), (KSODAT, INTGRI (78) ) 
EQUIVALENCE(KFLAG, INTGRI (79) ), ( KCOMPF, INTGRI (80 ) ), 
a (KSORCE, INTGRI (81 )),(KLES1D, INTGRI (82)), 
a ( KLES2D, INTGRI ( 83 ) ) , (KLES3D, INTGRI ( 84) ), 

8 (KCOMPP, INTGRI (85) ), (KADJST, INTGRI (86) ) , 

8 ( KFLUX, INTGRI (87 ) ) , (KSHI FT, INTGRI (83 ) ) , 

a ( KOUTPT, INTGRI ( 89) ), (KDIF, INTGRI (90)), 
a (KCOMPU, INTGRI (91) ), (KCOMPV, INTGRI (92) ) , 
a ( KCOMPW, INTGRI ( 93 ) ),(KCOMPR, INTGRI (94) ) 

EQUIVALENCE (KELIN, INTGR( 100 ) ) , (MEANDF, INTGRC99) ), 
a (NUMCLS,INTGRC93)), 

a ( IRHOl, INTGR(97) ) , (IRH02, INTGR(96) ) , 
a ( IZH1 , INTGR( 95) ) , ( IZW2, INTGR( 94) ) , 
a (MGRID,INTGR(93)),( KHALL ,INTGR(92)), 

8 (IZPR1,INTGR(91) ), (IZPR2, INTGR(90) ), 
a ( ISTPR1 , INTGR(89 ) ) , (ISTPR2, INTGR(88) ) , 
a (NTPRIN,INTGR(87)) 

EQUIVALENCE ( IMDOT , INTGR( 86 ) ) , ( IHSAT , INTGR(85) ) , 
a (ICFIP,INTGR(84)),(NTABLE,INTGR(83)),( NTABVR, INTGR(82)) 
a (LINTAB,INTGR(81)),(NPRTAB,INTGR(80)), (NMON, INTGR(79)), 
a (ITAB(1),INTGR(71)),(MTABVR(1),INTGR(63)), 
a (NZPRIN,INTGR(62)), 

a (NXPRIN,INTGR(61)),(NYPRIN,INTGR(60)), 
a ( LDISTL , INTGR( 59 ) ) , ( TSTSWP , INTGR ( 58 ) ) , 

a ( KDBEXP , INTGR ( 57 ) ) , (KDBRHO, INTGR (56) ) , 

a ( KDBMDT , INTGR ( 55 ) ) 

EQUIVALENCE ( I LOOP1 , I NTGR ( 53 ) ) , ( I LOOPN , INTGR( 52 ) ) , 
a (IPBP,INTGR(51)),(IZPRF,INTGR(50)), 
a (IZPRL,INTGR(49)),(ISTPRF,INTGR(48)), 
a ( ISTPRL , INTGR( 47 ) ) , ( KDBGEN, INTGR( 46 ) ) , 
a (IVELF,INTGR(45)),(IVELL,INTGR(44)), 
a (IKEF,INTGR(43)),(IKEL,INTGR(42)), 
a ( IENTF, INTGR( 41 ) ) , ( IENTL , INTGR( 40) ) , 
a (ICNCF,INTGR(39)),( ICNCL ,INTGR(38)), 
a (NVEL, INTGR(37) ), (NKE, INTGR(36) ), 
a (NENT, INTGR(35) ), (NCNC, INTGR(34) ), 
a ( NZSTP , INTGR( 33 ) ) , ( NPRMNT , INTGR ( 32 ) ) 

DIMENSION SPARE1(20),XUDIST(30),YVDIST(30), ZWDIST ( 50 ) , 
a SIGMA(25),CRIT(25) , DTFALS( 25 ) , RESREF( 25 ) , 
a TITLE(25),FIINIT(25) ,TFRAC(30) 

EQUIVALENCE ( SPARE1 ( 1 ) , RE1 ( 1 ) ) , 
a (TFRAC(1),RE1(21)), (XULAST, RE1 (51 ) ) , 

a ( YVLAST , RE1(52) ), (ZHLAST ,RE1(53) ), 

a (XUDIST(1),RE1(54)),(YVDIST(1),RE1(84)), 
a ( ZWDIST( 1 ),RE1 ( 114 )), (SIGMA (1),RE1( 164)), 

a (CRIT(1),RE1(189)),(DTFALS(1),RE1(214)), 
a (RESREF(1),RE1(239) ), ( EMU1 , RE1 ( 264) ) , 
a (RH01,RE1(265) ), ( RH02, RE1 ( 266 ) ) , 
a (TLAST,RE1(267)),(CFIPS,RE1(268)) 
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: MEL3SAT FORTRAN A1 



EQUIVALENCE ( AMDOT, RE1 ( 269 ) ) , < FIINIT( 1 ) , RE1 ( 270) ) , 

8 (RELAXP,RE1(295) ), (TITLE(l) , REK296) ), 

& (DT,RE1(321)), 

a (RINNER,ZWDIST(50)), ( SNALFA, ZWDIST( 99 ) ) , 

& (PBAR,ZWDIST(98) ) 

EQUIVALENCE ( SL OEMU, RE( 10 0 ) ) , ( SLORHO, RE( 99 ) ) , ( RLXRHO, RE( 98 ) ) 
& (RHOMAX,RE(97)), ( RHOMIN , RE( 96 ) ) , 

& ( EMUMAX, RE( 95 ) ) , ( EMUMIN, RE( 99 ) ) , 

8 ( TKEMAX, REC93) ), ( TKEMI N , RE( 92 ) ), 

a ( RLXPZ, RE( 91 ) ) , 
a (ABUOY,RE(90)), ( HREF, RE( 89 ) ) , 

8 ( AGRAVX, RE( 88 ) ) , ( AGRAVY , RE( 87 ) ) , 

8 ( AGRAVZ, RE (86)), (ATIME,RE(85)) 

EQUIVALENCE ( RLXPXY, RE( 89 ) ) , 
a ( EPSMAX, REC83) ) , ( EPSMIN , RE( 82 ) ) , 
a (ARH01,RE(81)),(BRH01,RE(30)),(CRH01,RE(79)) 

EQUIVALENCE ( AZ.W2 , RE( 7 3 ) ) , ( BZW2 , REC 7 A ) ) , 
a (PINT,RE(75)),(CZW2,RE(76)), 
a (DZ,RE(73)),(ZWL,RE(77)), ( ZW2M1T, RE( 72) ), 
a (VELMIN,RE(70) ), C VELMAX, REC 6 9 ) ) , 
a ( FALCOM, REC 63 ) ) , (CMDOT , REC 67 ) ) , 

a (RLXMDT,RE(66) ), ( AMDTMX, REC 65 ) ) , C AMDTMN, REC 6 A ) ) , 
a (RADMAX,RE(63)), ( RADMI N , REC 6 2 ) ) , 

a (PRESS, RE(61)), (CPI, REC60)), 
a (CP2,RE(59)),(CP3,RE(58)),(FSTOIC,RE(57)), 
a (ARRC0N,RE(56)), ( PREEXP ,RE(55)),(CA1I,RE(54)) 

EQUIVALENCE (CA2I , RE( 53 ) ) , ( CA13, RE( 52 ) ) , ( CA23, RE( 51 ) ) , 
a (H1SAT,RE(50)),(H2SAT,RE(<49)), 
a ( EMULAM, RE( 98 ) ) 

EQUIVALENCE ( TWPRCN, RE( 92 ) ) , ( PSATEX, RE( 91 ) ) 

EQUIVAL ENCE( EL EXP, RE( 39 ) ) , ( EWALL , RE( 38 ) ) , 

8 (SIGE,RE(37)),(SIGK,RE(36)),(TAUDK,RE(35)) 

DIMENSION LOCI ( 6 ) , LCC2( 6),L0C3(6),L0C9(6),L0C5(6), 
a LOC6(6),LOC7(6),LOC3(6),LOC9(6),LOC10(6) 

EQUIVALENCE ( LOCREG( 1 ) , LOCI ( 1 ) ) , ( LOCREG(7 ) , L0C2( 1 ) ) 

( LOCREG( 13 ) , L0C3( 1 ) ) , ( LOCREG( 19 ) , L0C9( 1 ) ) 
(L0CREG(25),L0C5(1>>, ( LOCREG( 31 ) , LOC6 ( 1 ) ) 

( LOCREG( 37 ) , L0C7 ( 1 ) ) , ( L OCREG( 93 ) , L0C8 ( 1 ) ) 
(L0CREG(99),L0C9(1)),(L0CREG(55), LOCI 0(1)) 
XUFRAC(30),YVFRAC(30) , ZWFRAC( 30 ) , ENDIT( 25 ) 
XFRAC(30) , YFRAC( 30),ZFRAC(30) 



EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
DIMENSION 
DIMENSION 



EQUIVALENCE( XUDIST ( 1 ) , XUFRAC( 1 ) , XFRAC( 1 ) ) , 
a (YVDIST(1),YVFRAC(1),YFRAC(1)), 

a ( ZWDIST ( 1 ) , ZWFRAC( 1 ) , ZFRAC( 1 ) ) , ( LITCNC, L ITC) , 

a (ENDIT(1),CRIT(1)),(RLXP,RELAXP) 

EQUIVALENCE ( NPRMON, LITER(l)) 

CLUDE 9 , CMNGRFIC . FTN/G 

LOGICAL GRAPHS, ORTHOG, ANTSYM 

COMMON / GRAF3 / GRAPHS, ORTHOG, ANTSYM, NPRT 

REAL NAME 

COMMON / GRAF9 / NAME (25),ITITL (5) 

COMMON/CPI/IPWRIT, IDUM(293) 

DIMENSION GDTAPE(3),DFAULT(9) 

DIMENSION ARRAY1 ( 309 ) , ARRAY2( 199 ) , ARRAY3( 921 ) 

LOGICAL ARRAY1 , LSPDA, WRT, RD, NAMLST 
INTEGER ARRAY2, XPL AN E, YP LANE, ZP LANE 

INTEGER P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,EP,H1,H2,H3,C1,C2, 
aC3,C9 

REAL NORTH, LOW 

EQUIVALENCE ( ARRAY1 ( 1 ), CARTES ),( ARRAY2( 1 ), NX) 

EQUIVALENCE ( ARRAY3( 1 ) , SPARE1 ( 1 ) ) , (Ml , R1 ) , ( M2, R2) 
EQUIVALENCE ( LSTRUN, INTGR( 12 )),( NAMLST, LOGIC( 88 ) ) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 1 ENDS. 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 1 STARTS: 
GRAFFIC ARRAYS DIMENSIONED AS NEEDED... 

COMMON/GRAF1/PHIK1) /GRAF2/PHI2( 1 ) 

POROSITY a SPECIAL DATA ARRAYS DIMENSIONED AS NEEDED... 
DIMENSION PE( 1, 1, 1 ), PN( 1,1,1), PH(1,1,1),PC( 1,1,1) 

DIMENSION LSPDA(1),ISPDA(1),RSPDA(1) 

USER PLACES HIS VARIABLES, ARRAYS, EQUIVALENCES ETC. HERE. 
DATA TITEM0/9HTEM0/ 

DATA TITEMP/9HTEMP/ 



MEL 01950 
MEL 0196 0 
MEL 0197 0 
MEL01980 
MEL01990 
MEL01500 
MEL 01510 
MEL 01520 
MEL01530 
MEL 01590 
MEL01550 
MEL0156 0 
MEL0157 0 
MEL 01580 
MEL01590 
MEL01600 
MEL01610 
MEL 01620 
MEL01630 
MEL01690 
MEL01650 
MEL 0166 0 
MEL 0167 0 
MEL 01680 
MEL01690 
MEL 017 00 
MEL01710 
MEL 017 20 
MEL01730 
MEL 01790 
MEL01750 
MEL 017 6 0 
MEL 01770 
MEL01780 
MEL01790 
MEL 018 0 0 
MEL01810 
MEL 01820 
MEL 01830 
MEL 01890 
MEL 01850 
MEL 0186 0 
MEL 0187 0 
MEL01880 
MEL 01890 
MEL 01900 
MEL 01 91 0 
MEL 01920 
MEL 01930 
MEL 01990 
MEL01950 
MEL 0196 0 
MEL 0197 0 
MEL 01980 
MEL01990 
MEL 020 0 0 
MEL02010 
MEL 02020 
MEL 02030 
MEL 02090 
MEL 02050 
MEL 0206 0 
MEL 0207 0 
MEL02080 
MEL02090 
MEL02100 
MEL02110 
MEL 02120 
MEL02130 
MEL02190 
MEL02150 
MEL 0216 0 



37 



FILE: MEL3SAT FORTRAN A1 



DATA TITDEN/4HRH01/ 

DATA NLSP,NISP,NRSP/1,1,1/ 

C USER PLACES HIS DATA STATEMENTS HERE. 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 1 ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 2 STARTS: 

C 

CHAPTER 2 SET CONSTANTS, AND ARRANGE FILE MANIPULATIONS. 

C 

C PLEASE DO NOT ALTER, OR RE-SET, ANY OF THE REMAINING 
C STATEMENTS OF THIS CHAPTER. 

DATA CELL, EAST, WEST, NORTH, SOUTH, HIGH, LOW, VOLUME/ 

8 0.,1.,2.,3.,4.,5.,6.,7. / 

DATA P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,KE,EP,H1,H2,H3,C1,C2, 
&C3,C 4/ 1,2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20/ 

DATA FIXFLU,FIXVAL,0NLYMS,WALL/1.E- 10,1. El 0,0. 0,-1 0.0/ 

DATA IPLANE,XPLANE,YPLANE,ZPLANE/0,1,2, 3/ 

DATA WRT , RD, DFAULT/ .TRUE. , . FALSE ., 4HDEFA, 4HULT ., 4HDTA/ , 1HG/ 
DATA GDTAPE/4HGUSI,4HE1 . D.2HTA/ 

DATA NL DATA, NI DATA, NRDATA/3 09, 194,421/ 

DATA NLCREG . NTCVRG/60 , 350/ 

CD CALL SETGRF 

CALL TAPES( 1 0 , GDT APE, 3,1, 4XNRDATA) 

c READ DEFAULT FILE IF BLOCKDATA ABSENT 

IF( INTGR1 (29) .NE.10) GO TO 2 

CALL WRIT40(40HDATA ESTABLISHED IN BLOCK DATA. ) 

GO TO 3 

2 CALL DEFLT 

CD 2 CALL TAPESC1,DFAULT,4,2,4*NRDATA) 

CD CALL DATAIO( RD, 1 ) 

CALL WRIT40(40HDATA TAKEN FROM DEFAULT. DTA ON GROUP A/C) 

3 CALL WRIT40C40HFILE MODSTL.FTN IS THE SATLIT USED. ) 

C 

CHAPTER 3 DEFINE DATA FOR NRUN RUNS. 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 2 ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 2 STARTS: 

C— GROUP 41 MULTI-RUNS : RUN( 1-30 ) < . T . , 29* . F . > 

C 

DO 410 11=1,1 
410 RUNCII )= .TRUE. 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 2 ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 3 STARTS: 
DO 10 IRUN=1 , 30 

I F C .NOT.RUN(IRUN)) GO TO 10 
NRUN=NRUN+1 
LSTRUN=IRUN 
10 CONTINUE 

DO 999 IRUN=1 , LSTRUN 

IF( .NOT.RUN(IRUN)) GO TO 999 
INTGRC 1 1 ) = IRUN 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 3 ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 3 STARTS: 

C ALL INTEGER VARIABLES ARE DEFAULTED TO 0, AND REAL VARIABLES 

C TO 0.0, UNLESS OTHERWISE INDICATED. 

C E.G. BY VARI ABL E<1 0> , OR <10. 0> AS APPROPRIATE. 

C THE DEFAULT SETTINGS OF ALL LOGICAL VARIABLES ARE ALWAYS 
C INDICATED, E.G. VARI ABL E< .T . >, OR VARIABL E< . F . > . 

C 

C RUN1 

C 

C GROUP 1 . FLOW TYPE : 

C PARAB<.F.>,CARTES<.T.>,ONEPHS<.T.> 

ONEPHS=. FALSE. 

C 

C GROUP 2. TRANSIENCE : 

C STEADY<.T.>,ATIME,LSTEP<1>,FSTEP<1> 

C TLAST<1 . El 0>, TFRACC 1-30 ) <30*1 .> 

STEADY=. FALSE. 

C LSTEP=30 

C ATIME=0 . 

C TLAST=50E-4 

C SERVICE SUBROUTINE FOR ’NT’ POWER-LAW TIME STEPS: 
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CALL GRDPWRC0,30,50E-4, 1 . ) 



- GROUP 3. X-DIRECTION : 

NX<1 > , XULAST<1 . 0> , XFRACC 1-30 ) 

SERVICE SUBROUTINE FOR POWER-LAW GRID: 
CALL GRDPWRC 1 , NX, XULAST , 1 . ) 



- GROUP 4. Y-DIRECTION : 

NY<1 > , YVLAST <1 . 0 > , YFRACC 1 -30 ) , RI NNER , SNAL FA 
SERVICE SUBROUTINE FOR POWER-LAW GRID: 

NY=1 00 

CALL GRDPWRC 2, 100, A . E-3 , 1 . ) 



- GROUP 5. Z-DI RECTI ON : 

NZ<1 > , ZWL AST <1 . 0>, ZFRACC 1-30 ) 

SERVICE SUBROUTINE FOR POWER-LAW GRID: 
CALL GRDPWRC 3, NZ, ZWL AST, POWER) 



- GROUP 6. MOVING GRID : 

MGRID, IZW1, IZW2, AZW2, 3ZW2, CZW2, PINT , ZW2M1T 



- GROUP 7. BLOCKAGE: BLOCK< . F . > , IPLANE, IPWRIT 

XSET CONSTANT POROSITIES OVER SUB-DOMAINS USING: 

CALL CONPOR(IR,TYPE,VALUE,IXF,IXL,IYF,IYL,IZF,IZL), WHERE: 

IR = RUN SECTION NUMBER, E.G. 1 FOR RUN1 SECTION; 1 TYPE* = EAST, 
WEST, NORTH, SOUTH, HIGH, LOW & CELL. ’ VALUE ' =WANTED POROSITY 
OVER REGION IXF, . . . IZL . 

^DIMENSION ARRAYS PEC NX, NY, NZ) , PNC NX, NY, NZ) , PHC NX, NY, NZ) , 8 
PCC NX , NY , NZ ) ABOVE. 

*FOR FULLY-BLOCKED CELLS CIE. T VALUE T = 0.0) USER NEED SET ONLY 
THE 1 CELL 1 POROSITY CTO ZERO), AS CELL-FACE AREAS ARE THEN 
AUTOMATICALLY ZEROED. 

*FOR SATELLITE PRINTOUT OF ALL POROSITIES IN DOMAIN, 1 IPLANE 1 = 
XPLANE YPLANE OR ZPLANE, FOR DESIRED CROSS-SECTI ON DIRECTION. 

*FOR EACH 1 TYPE 1 A MAXIMUM OF 10 CALLS TO CONPOR IS ALLOWED, 
BUT IF REQUIREMENTS EXCEED THIS PROVISION SET BLOCK=.T. 8 
IPWRIT = -l , AND SET POROSITY ARRAYS EXPLICITLY HERE AS WANTED. 
IN THIS CASE, THE USER MUST SET A L L ELEMENTS OF 
ARRAYS PE, PN, PH, PC C MANY MAY BE 0.0 OR 1.0). HE MAY USE: 
CALL CRCPARRAY, VALUE, IXF, IXL , IYF, IYL , IZF, IZL , NX, NY, NZ) 

ANY NUMBER OF TIMES, TO SET 1 PARRAY* C= PE, ETC.) TO 
1 VALUE 1 OVER RANGE IXF TO IXL, IYF TO IYL, IZF TO IZL. 

XCONPOR MUST NOT BE USED IN CONJUNCTION WITH EXPLICIT 
SETTINGS OF THE ARRAYS CINCLUDING SETTINGS VIA CR) . 



- GROUP 8. DEPENDENT VARIABLES TO BE SOLVED FOR OR STORED : 
S0LVARC1-25X25X.F.>,ST0VARC1-25X25X.F.>,C0NC1C1-4X4X.T.> 
USE FOLLOWING NAMED INTEGERS FOR ARRAY ELEMENTS 1-20: 
P1,PP,U1,U2,V1,V2,W1,W2,M1,M2,RS,KE,EP,H1,H2,H3,C1,C2,C3,C4 . 
SOLVAR C P 1 ) = . TRUE . 

SOLVAR C PP ) = . TRUE . 

SOLVAR CV1X.TRUE. 

SOLVAR C V2 ) = . TRUE . 

SOLVAR C Rl ) = . TRUE . 

SOLVAR C R2 ) = . TRUE . 

SOLVAR C RS ) = . TRUE « 

SOLVAR C HI ) = . TRUE . 

SOLVAR C H2 ) = . TRUE . 

SOLVAR C Cl ) = . TRUE . 

SOLVAR C C2 ) = . TRUE . 

SOLVAR C C3 ) = . TRUE . 

STOVAR C U1 ) = . TRUE . 

STOVAR C U2 ) = . TRUE . 

STOVAR C 24 ) = .TRUE . 



- GROUP 9. VARIABLE LABELS « 

TITLEC1-25X2HP1,2HPP, 2HU1 , 2HU2, 2HV1 , 2HV2, 2HW1 , 2HW2, 2HR1 , 
2HR2, 2HRS , 2HKE, 2HEP , 2HH1 , 2HH2 , 2HH3 , 2HC1 , 2HC2, 

2HC3 , 2HC4 , 2HRX, 2HRY, 2HRZ, 2*4H*xx*> 

TITLECU1) =TITEMO 
TITLECU2XTITEMP 
TITLECPPXTITDEN 
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GROUP 10 PROPERTIES: 

: IRH01<1>,IRH02<1>,RH01<1 .0>,RH02<1 . 0>, 

l ARH0K1 . 0>, BRH0X1 . 0>,CRH01<1 . 0> 

: IEMU1<1>, EMUK1 . 0>, EMULAM<1 . E-10> 

: IHSAT,H1SAT,H2SAT,PSATEX<1. 0> 

; SIGMAC 1-25 XI. 0,2. 0,1. , 1 . E10, 1 . , 1 . E10, 1 . , 1 . E10, 

; 4X1. 0,1. 314, 1.0,1.E10, 10*1. 0> 

IRH01=-1 
IEMU1=-1 
I RH02=-1 
SIGMAC 6 ) = 1 . El 0 
SIGMAC 1 5 ) =1 . El 0 
SIGMA C HI ) =~1 
SIGMAC Cl ) = -l 
SIGMAC C2 ) =-l 
SIGMAC C3 ) =-l 



GROUP 11 INTER-PHASE TRANSFER PROCESSES : 

; ICFIP,CFIPS, IMD0T,CMD0T,CA1I<1 . E6>, CA2I<1 . E6> 
ICFIP= 1 
CFIPS=5E07 
IMD0T=-1 



GROUP 12 SPECIAL SOURCES : 

: ISPCSOC 1 -25 ),AGRAVX, AGRA VY, AGRAVZ, ABUOY, HREF 



GROUP 13 INITIAL FIELDS : 
: FI I NIT C 1-25X25x1 . E-l 0> 

FI INITC 6 ) = 1 . E-10 
FIINITC15)=1 .E-10 
FIINITC 1 ) =100000. 
FIINITC5)=10101 
FIINITC9)=10101 
FIINIT C 14 ) =101 01 
FIINITC 17 )=10101 
FIINIT C 18 ) =101 01 
FIINITC19)=10101 
FIINIT C 1 0 ) =0 . 0 



GROUP 1 A BOUNDARY/INTERNAL CONDITIONS : 

: ILOOP1,ILOOPN,XCYCLE<.F.>,PBAR,REGIONC1-10X10X.T.> 

: xn.B. ALL 10 REGIONS ARE DEFAULTED .TRUE.. THE USER SHOULD 

: SET REGIONCIX. FALSE. FOR UNUSED REGIONS 'I'. 

DO 140 1=1,10 
140 REGIONCIX. FALSE. 

GROUP 15 TO 24; REGIONS 1 TO 10 

: ONLY THOSE REGIONS ARE ACTIVE WHICH ARE SPECIFIED BY THE 

: USER, PREFERABLY BY WAY OF:- 

: CALL PLACEC IREGN, TYPE, IXF,IXL,IYF,IYL,IZF,IZL) & 

: CALL COVALC IREGN, VARBLE, COE FF, VALUE) 

CALL PLACEC1,S0UTH,1,NX,1,1,1,NZ) 

CALL PLACEC3,NORTH,1,NX,100,100,1,NZ) 

CALL C0VALC3, Ml, FIXVAL, 100000. ) 

CALL C0VALC3, M2, FIXVAL, 100000. ) 



: GROUP 25 GROUND STATION : 

: GROSTA<.F.>, N AML ST < . F . > 

GROSTA= . TRUE . 

: XNAMLST ACTIVATES NAMELIST IN GROUND. 



GROUP 26 SOLUTION TYPE AND RELATED PARAMETERS : 
: WHOLEP< . F. >, SUBPST< . F . >, DONACC< . F . > 



: GROUP 27 SWEEP AND ITERATION NUMBERS : 

) FSWEEP<1>,LSWEEP<1>,LITHYD<1>,LITC<1>,LITKE<1>,LITH<1>, 

LSWEEP=20 

: LITERC 1-25X9X1, -1,1 5xi> 

LITERC2X 50 
LITERC 17 )=20 
L ITERC 18 ) =20 



MEL03610 
MEL03620 
MEL03630 
MEL03640 
MEL03650 
MEL03660 
MEL03670 
MEL03680 
MEL 036 9 0 
MEL03700 
MEL 037 1 0 
MEL03720 
MEL03730 
MEL03740 
MEL03750 
MEL 037 6 0 
MEL 0377 0 
MEL03780 
MEL 037 90 
MEL 03800 
MEL03810 
MEL 03820 
MEL 03830 
MEL03840 
MEL03850 
MEL03860 
MEL 0387 0 
MEL03880 
MEL03890 
MEL03900 
MEL 0391 0 
MEL 03920 
MEL03930 
MEL 0394 0 
MEL03950 
MEL 0396 0 
MEL0397 0 
MEL03980 
MEL03990 
MEL04000 
MEL 0401 0 
MEL04020 
MEL04030 
MEL 04040 
MEL04050 
MEL 0406 0 
MEL 0407 0 
MEL04080 
MEL04090 
MEL04100 
MEL04110 
MEL04120 
MEL04130 
MEL04140 
MEL 04 150 
MEL04160 
MEL 0417 0 
MEL04180 
MEL04190 
MEL04200 
MEL04210 
MEL 04220 
MEL04230 
MEL 04240 
MEL04250 
MEL04260 
MEL04270 
MEL04280 
MEL04290 
MEL04300 
MEL04310 
MEL 04320 
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. E : MEL3SAT FORTRAN A1 



L ITER( 1 9 ) =20 

IVELF<1>,NVEL<1>,IVELL<10000>, 
IKEF<1>,NKE<1>,IKEL<10000>, 
IENTF<1>,NENT<1>,IENTL<10000>, 
ICNCF<1 > , NCNC<1 >, ICNCL<1 00Q0> , 
IRH01F<1>,NRH01<1>, IRHOl L<1 000Q>, 
IRH02F<1>,NRH02<1>,IRH02L<10000> 



- GROUP 28 TERMINATION CRITERIA : 

ENDITC 1-25X9X1 . E-l 0 , 0 . 5 , 15*1 . E-l 0> 



- GROUP 29 RELAXATION : 

RLXP<1 . >,RLXPXY<1 . >,RLXPZ<1 . >,RLXRH0<1 . >,RLXMDT<1 . >, 

RLXRHO=0 . 2 

DTFALSCC1X4.E-5 

DTFAL S ( C2 ) =4 . E-5 

DTFALSCC2) =4 . E-5 

DT FA LS( 3-25X23X1 . E10> 



- GROUP 30 LIMITS : 

VELMAXCl . El 0> , VELMIN<-1 . El 0> , RHOMAX< 1 . El 0>, RH0MIN<1 . E-10>, 
TKEMAX<1 . El Q> , TKEMI N<1 . E-l 0> , EMUMAX<1 . El 0> , EMUMI N<1 . E-10>, 
EPSMAX<1 . E10>, EPSMIN<1 . E-10> , AMDTMX<1 . El 0> , AMDTMN<-1 . E10> 



- GROUP 31 SLOWING DEVICES : SL0RH0<1 . > , SL0EMU<1 . > 



- GROUP 32 PRINT-OUT OF VARIABLES : 

PRINT ( 1-25X . T . , . F. ,23X.T.>,SUBWGR<.F.> 
PRINT(U1)=.TRUE. 

PRINT(U2)= .TRUE. 

PRINT(PP)=.TRUE. 

PRINT(24)= .TRUE. 



- GROUP 33 MONITOR PRINT-OUT : 

IYM0N=25 

IXM0N<1 > , IYM0N<1> / IZM0N<1 > , NPRM0N<1 > , NPRMNT<1 > 



- GROUP 34 FIELD PRINT-OUT CONTROL : 

NPRINT<100>,NTPRIN<100>,NXPRIN<1>/NYPRIN<1>, NZPRIN<1 > , 

IZPRF<1>,ISTPRF<1>,IZPRL<10000>,ISTPRL<10000> 

NUMCLS<10>,KOUTPT 

KOUTPT=-l 

NTPRI N=1 



- GROUP 35 TABLE CONTROL : 

TABLES<.F.>,NTABLE,NTABVR,LINTAB,NPRTAB,NMON, 

ITAB(l-8),MTARVR(l-8) 



GROUP 36-38 ARE NOT DOCUMENTED IN THE INSTRUCTION 
MANUAL AND ARE INTENDED FOR MAINTENANCE PURPOSES ONLY 
- GROUP 36 DEBUG PRINT-OUT SLAB AND TIME-STEP * 
IZPRK1>,IZPR2<1>,ISTPR1<1>,ISTPR2<1> 



- GROUP 37 DEBUG SWEEP AND SUBROUTINES : 

KEMU,KMAIN,KINDEX,KGEOM,KINPUT,KSODAT,KCOMPF,KSORCE, 
KS0LV1,KS0LV2,KS0LV3,KC0MPP,KADJST,KFLUX,KSHIFT,KDIF, 
KCOMPU, KCOMPV, KCOMPW, KCOMPR , KWALL , KDBRH0<-1> , KDBEXP , KDBMDT 
KDBGEN 



- GROUP 38 MONITOR, TEST, AND FLAG : 

MONITR< . F . >, FL AG< . F . >, TEST <.T.>,KFLAG<1> 
END OF MAINTENANCE-ONLY SECTION 



- GROUP 39 ERROR AND RESIDUAL PRINT-OUT : 

I ERRP<10 00>,RESREF( 1,3-24X25X1 . > , RESMAP< . F , >, 

RESIDC1-25X2X.F. , 23X .T . >, KOUTPT 

I ERRP=LSWEEP 

RESREFC 1 ) =5 . E-5 

RESREF(Vl)=5.E-5 

RESREF(V2)=5.E-5 

RESREF ( HI ) =40 . 

RESREF C H2 ) =40 . 



MEL 04330 
MEL 04340 
MEL 04350 
MEL 0436 0 
MEL 0437 0 
MEL 0438 0 
MEL 0439 0 
MEL 04400 
MEL 0441 0 
MEL 04420 
MEL 04430 
MEL 04440 
MEL 04450 
MEL 0446 0 
MEL 0447 0 
MEL 0443 0 
MEL 04490 
MEL 0450 0 
MEL 04510 
MEL 04 520 
MEL 04 53 0 
MEL 04540 
MEL 04550 
MEL 0456 0 
MEL 04 57 0 
MEL 0453 0 
MEL 0459 0 
MEL 046 0 0 
MEL 046 1 0 
MEL 046 20 
MEL04630 
MEL 046 40 
MEL 046 50 
MEL 046 6 0 
MEL 046 7 0 
MEL 046 3 0 
MEL 046 9 0 
MEL 047 0 0 
MEL 047 1 0 
MEL 047 20 
MEL 047 30 
MEL 047 40 
MEL 047 50 
MEL 0476 0 
MEL 0477 0 
MEL 047 8 0 
MEL 047 9 0 
MEL 04800 
MEL 048 1 0 
MEL 04820 
MEL 0483 0 
MEL 04840 
MEL 04850 
MEL 0486 0 
MEL 0487 0 
MEL 0488 0 
MEL 0489 0 
MEL 049 00 
MEL 04910 
MEL 04920 
MEL 04930 
MEL 04940 
MEL 04950 
MEL 0496 0 
MEL 0497 0 
MEL 0498 0 
MEL 0499 0 
MEL05000 
MEL 050 1 0 
MEL 05020 
MEL05030 
MEL05040 
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FILE: MEL3SAT FORTRAN A1 



RESREF(Cl)=5.E-5 

RESREF(C2)=5.E-5 

RESREF(C3)=5.E-5 



GROUP 40 SPECIAL DATA : L0GIC(1. .10), INTGR(1. .10),RE(21. .30), 
NLSP<1> > NISP<1>» NRSP<1>, SPDATA< . F. >, LSPDAd ) » ISPDAd ) , RSPDAd ) 
USE FIRST 10 ELEMENTS OF ARRAYS LOGIC S INTGR AND 21ST 
TO 30TH OF ARRAY RE FOR TRANSFERRING SPECIAL DATA FROM 
SATELLITE TO GROUND, BUT IF REQUIREMENTS EXCEED THIS 
PROVISION SET SPDATA = .T., AND DIMENSION ARRAYS LSPDA , 

ISPDA , RSPDA ABOVE AND IN GROUND AS NEEDED, AND SET HERE 
NLSP , NISP, NRSP TO DIMENSION VALUES. 



— GROUP 42 RESTARTS AND DUMPS : SAVEM< . F . >, RESTRT < . F . >, KINPUT 



GROUP 43 GRAFFIC : 

GRAPHS< .F.>,ORTHOG<.T. >, ANTSYM , NPRT<1 > , ITITL<5*4H****> 
FOR A GRAFFIC RUN, DIMENSION PHI1 3 PHI2 AS FOLLOWS: 
PHI1(NX*NY*NZ*NM) 

PHI2((NX+2)*(NY+2)*(NZ+2)*(NM+IBLK)) , WHERE 
NM=NO . OF VARIABLES STORED + DENSITY ( -I ES ) 

IBLK=0 IF BL0CK=.FALSE.,=4 IF A 3D RUN, 

=3 IF A 2D.YZ RUN. 



MEL05050 
MEL 0506 0 
MEL 0507 0 
MEL05080 
MEL05090 
MEL05100 
MEL05110 
MEL 05120 
MEL05130 
MEL05140 
MEL 05150 
MEL05160 
MEL05170 
MEL05180 
MEL05190 
MEL05200 
MEL0521 0 
MEL05220 
MEL 05230 
MEL05240 
MEL 05250 
MEL05260 
MEL 0527 0 
MEL05280 





IFdRUN. EQ.l) GO TO 900 


MEL05290 


c — 


RUN2 


MEL05300 




IFdRUN. EQ. 2) GO TO 900 


MEL 0531 0 


c — 


RUN3 


MEL05320 




IFdRUN . EQ . 3 ) GO TO 900 


MEL05330 


c — 


RUN4 


MEL 05340 




IFdRUN. EQ. 4) GO TO 900 


MEL 05350 


c — 


RUN5 


MEL 0536 0 




IFdRUN. EQ. 5) GO TO 900 


MEL0537 0 


c — 


RUN6 


MEL05380 




IFdRUN. EQ. 6) GO TO 900 


MEL05390 


c — 


RUN7 


MEL 05400 




IFdRUN. EQ. 7) GO TO 900 


MEL05410 


c— - 


RUN8 


MEL 05420 




IFdRUN. EQ. 8) GO TO 900 


MEL 05430 


c — 


RUN9 


MEL05440 




IFdRUN. EQ. 9) GO TO 900 


MEL05450 


c — 


RUN10 


MEL 0546 0 




IFdRUN. EQ. 10) GO TO 900 


MEL0547 0 


c — 


RUN11 


MEL05480 




IFdRUN.EQ.il) GO TO 900 


MEL05490 


c — 


RUN12 


MEL05500 




IFdRUN. EQ. 12) GO TO 900 


MEL05510 


c — 


RUN13 


MEL05520 




IFdRUN. EQ. 13) GO TO 900 


MEL 05530 


c — 


RUN14 


MEL05540 




IFdRUN. EQ. 14) GO TO 900 


MEL05550 


c— 


RUN15 


MEL05560 




IFdRUN. EQ. 15) GO TO 900 


MEL 0557 0 


c 


RUN16 


MEL05580 




IFdRUN. EQ. 16) GO TO 900 


MEL05590 


c— 


RUN17 


MEL05600 




IFdRUN. EQ. 17) GO TO 900 


MEL05610 


c— 


RUN18 


MEL05620 




IFdRUN. EQ. 18) GO TO 900 


MEL05630 


c— 


RUN19 


MEL05640 




IFdRUN. EQ. 19) GO TO 900 


MEL05650 


c— 


RUN20 


MEL05660 




IFdRUN. EQ. 20) GO TO 900 


MEL05670 


c— 


RUN21 


MEL05680 




IFdRUN. EQ. 21) GO TO 900 


MEL05690 


c— 


RUN22 


MEL05700 




IFdRUN. EQ. 22) GO TO 900 


MEL05710 


c— 


RUN23 


MEL05720 




IFdRUN. EQ. 23) GO TO 900 


MEL 057 30 


c— 


RUN24 


MEL05740 




IFdRUN. EQ. 24) GO TO 900 


MEL05750 


c— - 


RUN25 


MEL05760 
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I F( IRUN . EQ . 25) GO TO 900 

■ RUN26 

I F( IRUN . EQ . 26 ) GO TO 900 

- RUN27 

I F C IRUN . EQ . 27 ) GO TO 900 

■ RUN28 

IF( IRUN . EQ . 28 ) GO TO 900 

* RUN29 

I F( IRUN . EQ . 29 ) GO TO 900 

- RUN30 

IFC IRUN . EQ . 30) GO TO 900 
10 CONTINUE 

• ALL RUNS 

:XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX USER SECTION 3 ENDS. 
:XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX STANDARD SECTION 4 STARTS 

IRITE GENERAL DATA ON TO THE GUSIE1.DTA TAPE, ETC... 

IF(SPDATA) CALL WRTSPCC L SPD4 , NL3P , ISPDA , NISP , RSPDA , NRSP ) 
IF(BLOCK) CALL WRTPQR( PE, PN , PH , PC, NX, NY, NZ, IPLANE) 

OLD PRACTICES RETAINED FOR REFERENCE: 

IF(SPDATA) CALL SPCDAT(IRUN) 

I F ( BLOCK) CALL PORDAT ( IRUN) 

I FC GRAPHS ) CALL SORT(IRUN) 

I F ( RESTRT ) GO TO 902 
DO 901 INDVAR= 1 , 25 

IF(IFIX(FIINIT(INDVAR)+0.1) .NE. 10101) GO TO 901 
CALL FLDDAT ( IRUN) 

GO TO 902 

1 CONTINUE 

2 CALL DATAIOCWRT , 10 ) 

IF(MONITR) CALL DATAIO ( WRT , -6 ) 

CONTINUE 

STOP 

END 

RECTI VEX* FLDDAT 
XFILE NAME: MODFLD.FTN 

^ABSTRACT: SATELLITE MODEL FLDDAT SUBROUTINE. 

^DOCUMENTATION: PHOENICS INSTRUCTION MANUAL (SPRING 1983). 
SUBROUTINE FLDDAT ( IRUN ) 

FIELDS IS USED TO SPECIFY NON-UNIFORM INITIAL FIELDS. 

EARTH SETS UNIFORM INITIAL FIELDS TO FI INITCMPHI ) , EXCEPT 
WHEN FINIIT(MPHI)=10101 . 0 WHICH IS THE SIGNAL THAT EARTH 
READS THE NON-UNIFORM FIELDS TAPE THE CONTENTS OF WHICH 
ARE SET HERE. 

IT IS ESSENTIAL TO PROVIDE SETTINGS FOR ALL THOSE MPHI’S 
FOR WHICH FI I NIT ( MPHI ) HAS BEEN SET TO 10101. THE SETTING 
ORDER MUST FOLLOW THE STANDARD ORDER OF BLOCK DATA. 
NOTE: EARTH PRINTS OUT THE INITIAL FIELDS IF KOUTPT=-l . 

PLEASE NOTE THAT THE PHI ARRAY M U S T BE SET AS 
PHI(IY,IX)=. . .,NOT AS PHI ( IX, IY) = . . . 

ALSO, YOU MUST DIMENSION PHI TO THE MAXIMUM 
DIMENSIONS USED IN THE RUNS DEFINED IN SATLIT, VIZ. 

DIMENSION PHI ( NY MAX., NX MAX.) 



TER 1 



PRELIMINARIES 



CLUDE 9 , CMNGUSSI . FTN/G 
LOGICAL LOGIC1, LOGIC 
DIMENSION L OGICC 100) 

COMMON/ LDATA/LOGI Cl (309) 

EQUIVALENCE ( LOGIC( 1 ) , LOGIC1 ( 210 ) ) 

DIMENSION INTGR(IOO) 

COMMON/IDATA/INTGR1 (194) 

EQUIVALENCE ( INTGR( 1 ) , INTGRK95) ) 

DIMENSION RE(100) 

COMMON/RDATA/RE1 ( 421 ) 

EQUIVALENCE ( RE( 1 ) , RE1 ( 322 ) ) 

COMMON/ BOUND/ L OCR EG ( 60), 

&TR1 , CP1R1 ( 7 ) , VP1R1 (7 ) , CP2R1 (5),VP2R1(5),CPNR1(5),VPNR1(5), 
&TR2,CP1R2(7),VP1R2(7),CP2R2(5),VP2R2(5),CPNR2(5),VPNR2(5), 



MEL 0 577 0 
MEL 05780 
MEL 057 90 
MEL 058 0 0 
MEL 05810 
MEL 05820 
MEL 058 3 0 
MEL 058 4 0 
MEL 05850 
MEL 0586 0 
MEL 0 587 0 
MEL 0 588 0 
MEL 0589 0 
MEL 0590 0 
MEL 059 1 0 
MEL 05920 
MEL 05930 
MEL 05940 
MEL 059 5 0 
MEL 0596 0 
MEL 0597 0 
MEL 0598 0 
MEL 0599 0 
MEL 06000 
MEL 06010 
MEL 06 020 
MEL 06 03 0 
MEL 06 040 
MEL 06 050 
MEL 06 06 0 
MEL 06 07 0 
MEL 06 03 0 
MEL 06 09 0 
MEL 06100 
MEL 06 1 1 0 
MEL 06 1 20 
MEL 06130 
MEL 06140 
MEL 06150 
MEL 06 16 0 
MEL 06 17 0 
MEL 06 180 
MEL 06190 
MEL 0620 0 
MEL 06210 
MEL 06220 
MEL 06230 
MEL 06240 
MEL 06 250 
MEL 0626 0 
MEL 0627 0 
MEL 0628 0 
MEL 0629 0 
MEL 06 30 0 
MEL 06310 
MEL 06320 
MEL 06330 
MEL 06 340 
MEL 06 350 
MEL 06 36 0 
MEL 0637 0 
MEL 06380 
MEL 06 390 
MEL 0640 0 
MEL 0641 0 
MEL 06420 
MEL 06430 
MEL 06 440 
MEL 06 450 
MEL 0646 0 
MEL 06 47 0 
MEL 06480 
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aTR3,CPlR3(7),VPlR3(7),CP2R3(5),VP2R3(5),CPNR3(5),VPNR3(5), 
aTR4,CPlR4(7),VPlR4(7),CP2R4(5),VP2R4(5),CPNR4(5),VPNR4(5), 
aTR5,CPlR5(7),VPlR5(7),CP2R5(5),VP2R5(5),CPNR5(5),VPNR5(5), 
8TR6,CP1R6(7),VP1R6(7),CP2R6(5),VP2R6(5),CPNR6(5),VPNR6(5), 
8TR7,CP1R7(7),VP1R7(7),CP2R7(5), VP2R7 ( 5 ) , CPNR7 ( 5) , VPNR7 ( 5) , 
aTR8,CPlR8(7),VPlR8(7) , CP2R8 ( 5 ) , VP2R8 ( 5 ) , CPNR8 ( 5) , VPNR8( 5) , 
aTR9,CPlR9(7),VPlR9(7),CP2R9(5),VP2R9(5) ,CPNR9(5),VPNR9(5), 
aTR10,CPlR10(7),VPlR10(7),CP2R10(5),VP2R10(5),CPNR10(5),VPNR10(5) 
DIMENSION TCVREG( 350 ) 

EQUIVALENCE ( TCVREG( 1 ) , TR1 ) 

C$INCLUDE 9 , GUSSEQUI . FTN/G 

LOGICAL CARTES, POLAR, SPDATA, SKEW, TWODYZ, ONEDZ, STOVAR( 25), 

& S0LVAR(25),PRINT(25),RESID(25),CMPRSS, CON EMU , LSP1 , 

& CONRHO, EMDOT, ONEPHS, INCOREC 10) , ISAVED,SAVEI , SAVEM, 

& RESTRT,XCYCLE,MONITR, REGION! 10), STEADY, WHOLEP, SLABPP, 

8 RAIN, BLOCKZ,PWSTAG, RUN( 30) , PLOT, RESMAP, FLAG, BLOCK, 

8 TEST, C0NC1C4), DISTIL 
LOGICAL SPRESS,PARAB, DON ACC, OVERLY, SACC, 

& GUSSIE, CATHY, CONNIE, CORA, ESTER, FLASH, FLORA, 

& FOCS,GENMIX, HESTER, PICAL 0, PL ANT, SPL ASH, HELGA, 

8 TACT, TI BALT, TOPS I, PAMELA , TABLES , WSTAG, 
a CONMOD, GROSTA,SUBPST, SUBWGR 
EQUIVALENCE ( CARTES , L0GIC1 ( 1 )),( POLAR, L0GIC1 ( 2) ) , 

8 (SPDATA, LOGICK 3)), (SKEW, LOGICK 4)), (TWODYZ, LOGICK 5)), 
a (0NEDZ,L0GIC1(6)),(ST0VAR(1),L0GIC1(7)), 
a ( SOL VAR(l), LOGICK 32) ) , ( PRINT (1 ) , L0GIC1 ( 57 ) ) , 
a (RESID(l), LOGICK 82)) , (CMPRSS , LOGICK 1 07 )) , 
a ( CONEMU, LOGICK 1 08) ), ( LSP1 , LOGICK 1 09 )) , 

8 (CONRHO, LOGICK 110) ), (EMDOT, LOGICK 111)), 
a ( ONEPHS, LOGICK 1 12 ) ) , ( INCOREC 1 ) , LOGICK 113)) 

EQUIVALENCE! I SAVED, LOGICK 123) ) , (SAVEI , LOGICK 124) ) , 
a (SAVEM, LOGICK 125)), (RESTRT, LOGICK 126) ) , 
a (XCYCLE, LOGICK 127)), (MONITR, LOGICK 128)), 
a (REGION! 1), LOGICK 129)), (STEADY, LOGICK 139)), 
a (WHOLEP, LOGICK 140) ), ( SLABPP, LOGICK 141 )) , 
a (RAIN, LOGICK 142)), (BLOCKZ, LOGICK 143)), 
a (PWSTAG, LOGICK 144) ),(RUN(1 ), LOGICK 145)), 
a (PLOT, LOGICK 175)) , (RESMAP, LOGICK 176) ) , 

a (FLAG, LOGICK 177)), (BLOCK, LOGICK 178)), (TEST, LOGICK 179)) 
EQUIVALENCE(CONCKl), LOGICK 206)) 

EQUIVALENCE (SPRESS, LOGIC! 100)) , (PARAB, LOGIC! 99) ), 
a ( DON ACC, L0GIC(98)) 

EQUIVALENCE ( OVERLY, LOGIC( 96 )),( SACC, LOGIC( 95) ) , 

8 (WSTAG, L0GIC(94)) 

EQUIVALENCE ( DISTI L , LOGIC( 90 ) ) 

EQUIVALENCE ( PAMELA , LOGIC( 67 )), (TOPSI , LOGIC( 66 )) , 
a (TIBALT,L0GIC(65) ), (TACT, LOGIC( 64)), (HELGA, LOGIC ( 63)), 

a ( SPLASH, LOGIC! 62) ),( PLANT, L0GIC(6 1) ) , ( PICALO , LOGIC(60 ) ) , 

a (HESTER,L0GIC(59)),(GENMIX,L0GIC(58)),(F0CS,L0GIC(57)), 
a (FLORA, LOGIC (56)), (FLASH, LOGIC( 55)), (ESTER, L0GIC(54)), 

8 ( CORA, L OGIC( 53)), (CONNIE, LOGIC( 52) ),( CATHY, L0GIC(51)), 

a (GUSSIE, LOGIC(50) ) 

EQUIVALENCE (TABLES , LOGIC( 48 )), (CONMOD, LOGIC( 47 )) , 

8 (GROSTA, L0GIC(46 ) ) , ( SUBPST , LOGIC( 45 ) ) 

EQUIVALENCE (SUBWGR, L0GIC(43) ) 

INTEGER FSTEP,FSWEEP,TSTSWP,ITAB(8),MTABVR(8) 

DIMENSION I S PCSO ( 25 ), LITER ( 25) 

EQUIVALENCE ( NX, INTGRK 1 ) ) , ( NY, INTGRK 2) ) , 
a (NZ, INTGRK 3 )),(ISPCS0(1), INTGRK 4)), 
a ( NREGN, INTGR1(29)),(NPHI,INTGR1(30)), 
a ( L ITKE, INTGRK 31 )),(LITHYD, INTGRK 32)), 
a (LITH, INTGRK 33)), (LITCNC, INTGRK34)), 
a (LITSLB, INTGRK 35 )),(NRUN, INTGRK 36)), 
a ( L ITER( 1 ), INTGRK 37 )),( FSTEP, INTGRK 62) ) , 
a ( FSWEEP, INTGRK 63 )), (LSTEP, INTGRK 64)) 

EQUIVALENCE ( LSWEEP , INTGRK 65) ) , ( NPRINT, INTGRK 66 )) , 
a (I ERRP, INTGRK 67) ),(IMAXP, INTGRK 68)), 
a ( IEMU1, INTGRK 69 )), ( IXMON, INTGRK7 0 ) ) , 
a (IYM0N,INTGR1(71)),(IZM0N, INTGRK 72)), 
a (KEMU, INTGRK 73 )),(KMAIN, INTGRK 74)), 
a (KINDEX, INTGRK 75) ), (KGEOM, INTGRK 76) ), 
a (KINPUT , INTGRK 7 7 ) ) , (KSODAT ,INTGR1(78)) 



MEL06490 
MEL06500 
MEL 0651 0 
MEL06520 
MEL06530 
MEL06540 
MEL 06550 
MEL 0656 0 
MEL 0657 0 
MEL06580 
MEL 06 59 0 
MEL06600 
MEL06610 
MEL 06620 
MEL06630 
MEL 06640 
MEL06650 
MEL 0666 0 
MEL0667 0 
MEL06680 
MEL06690 
MEL 067 00 
MEL067 1 0 
MEL 067 20 
MEL 067 30 
MEL06740 
MEL067 50 
MEL06760 
MEL 0677 0 
MEL06780 
MEL 0679 0 
MEL06800 
MEL 068 1 0 
MEL 06820 
MEL06330 
MEL 06840 
MEL 06850 
MEL06860 
MEL0687 0 
MEL06880 
MEL06890 
MEL06900 
MEL06910 
MEL06920 
MEL06930 
MEL 06940 
MEL06950 
MEL0696 0 
MEL0697 0 
MEL06980 
MEL06990 
MEL07000 
MEL07010 
MEL07020 
MEL07030 
MEL07040 
MEL07050 
MEL07060 
MEL07 07 0 
MEL 07 080 
MEL07090 
MEL 07 100 
MEL 07110 
MEL07120 
MEL07130 
MEL07140 
MEL 07 150 
MEL07160 
MEL07170 
MEL07180 
MEL07190 
MEL07200 
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EQUIVALENCECKFLAG, INTGR1 (79 ) ) , ( KCOMPF, INTGR1 (80) ) , 



(KS0RCE,INTGR1(81)), (KLES1D, INTGRK82) ) , 
(KLES2D, INTGR1 (83) ) , ( KLES3D, INTGR1 (84) ) 

(KADJST,INTGR1(86)), 
(KSHIFT/INTGRK88))/ 



(KCOMPP, INTGRK85)), 
(KFLUX/INTGRK87)) 



(KOUTPT,INTGR1(89)),(KDIF,INTGR1(90)), 



(KCOMPU, INTGR1 ( 91 ) ) , 
(KCOMPW, INTGRK93) ) 



(KCOMPV, INTGRK92) ) . 
( KCOMPR, INTGR1 ( 94 ) ) 



EQUIVALENCE ( KEL IN, INTGR( 100)) , (MEANDF, INTGR( 99 ) ) , 
(NUMCLS,INTGR(98)), 

( IRHOl , INTGR(97 ) ) , ( IRH02, INTGR(96)), 

( IZW1, INTGR(95)), ( IZN2 , INTGR( 94 ) ), 

(MGRID, INTGR(93) ), ( KUAL L , INTGR( 92 ) ), 
(IZPR1,INTGR(91)),(IZPR2,INTGR(90)) , 

(ISTPR1, INTGR(89) ) , ( ISTPR2 , INTGR( 88 ) ) , 
(NTPRIN,INTGR(87)) 

EQUIVALENCE ( IMLOT , INTGR(86 ) ) , ( IHSAT, INTGR(85) ) , 

(ICFIP, INTGR(84) ), (NTABLE, INTGR(83) ) , ( NTABVR, INTGR(82)) 
(L INTAB, INTGR(81) ), (NPRTAB, I NTGR ( 8 0 ) ) , ( NMON , INTGR( 7 9 ) ) , 
( ITAB(1 ), INTGR(71) ) , (MTABVR( 1) , INTGR(63) ) , 

(NZPRIN, INTGR( 6 2 ) ) , 

( NXPRIN, INTGR(61) ) , (NYPRIN, INTGR(60) ) , 

INTGR(59) ), (TSTSNP, INTGR(58 ) ), 
(KDBRH0,INTGR(56)), 



(LDISTL 

( KDBEXP, INTGR(57)) 

( KDBMDT ,INTGR(55)) 
EQUIVALENCE ( IL OOP1 , INTGR( 53 ) ) 



(IL00PN,INTGR(52)), 

( IZPRF, INTGR(50) ), 
(ISTPRF,INTGR(48)), 

( ISTPRL / INTGR( 47 ) ) , ( KDBGEN, INTGR( 46 ) ) , 
(IVELF,INTGR(45)), (IVELL, INTGRC44)), 

( IKEL , INTGR( 42) ) , 

( I ENTL , I NTGR ( 40 ) ) , 

( ICNCL ,INTGR(33)), 

( NKE, INTGR( 36 ) ) , 

( NCNC, INTGRC 34 ) ) , 

(NPRMNT, INTGR(32)) 

DIMENSION SPAREK20) , XUDIST( 30 ) , YVDIST( 30 ) ,ZWDIST(50) 
8 SIGMA (25),CRIT(25),DTFALS(25),RESREF(25), 

8 TITLE(25),FIINIT(25),TFRAC(30) 

EQUIVALENCE ( SPARE1 ( 1 ) , RE1 ( 1 ) ) , 



( IPBP , INTGR( 51 ) ) 

( IZPRL , INTGR( 49 ) ) , 



( IKEF, INTGR( 43) ) 

( IENTF, INTGR( 41 ) ) , 
(ICNCF,INTGR(39)) J 
( NVEL , INTGR( 37 ) ) , 
(NENT , INTGR(35) ) , 

( NZSTP, INTGR(33)) 



(TFRAC(1),RE1(21) ), ( XULAST , RE1 ( 51 ) ) , 
(YVLAST,RE1(52)),(ZWLAST,RE1(53)), 

(XUDIST(l) ,RE1(54) ) , (YVDIST( 1 ),RE1(84) ) , 
(ZWDIST(1),RE1(114)), (SIGMA(l) ,RE1(164)), 
(CRIT(1),RE1(189)),(DTFALS(1),RE1(214)), 

(RESREF(l) ,RE1(239) ) , ( EMU1 , RE1 ( 264 ) ) , 

(RH01,RE1(265)), (RH02,RE1(266)), 

(TLAST,RE1(267 )) , ( CFIPS , RE1 ( 268 ) ) 

EQUIVALENCE ( AMDOT , RE1 ( 26 9 ) ) , ( FI INIT ( 1 ) , RE1 ( 27 0 ) ) , 

8 (RELAXP,RE1(295) ), (TITLE(l) ,RE1(296) ), 

8 ( DT , RE1 ( 321 ) ) , 

8 (RINNER,ZWDIST(50)),(SNALFA,ZWDIST(49))> 

8 (PBAR,ZWDIST(48)) 

EQUIVALENCE ( SLOEMU , RE( 1 0 0 ) ) , ( SLORHO , RE( 99 ) ) , ( RLXRHO , RE( 98 ) ) , 



( RHOMAX, RE(97)),(RHOMIN,RE(96))j 
( EMUMAX, RE( 95) ) , ( EMUMIN , RE( 94 ) ) , 

( TKEMAX, RE( 93 ) ) , (TKEMIN, RE( 92) ) , 

( RLXPZ , RE( 9 1 ) ) , 

(ABUOY,RE(90)),(HREF,RE(89)), 

( AGRAVX, RE( 88 ) ) , ( AGRAVY , RE(87 ) ) , 

( AGRAVZ , RE( 86 ) ) , ( ATIME, RE( 85) ) 

EQUIVALENCE ( RLXPXY, RE( 84 ) ) , 

8 ( EPSMAX,RE(83)) , (EPSMIN,RE(82) ) , 

8 (ARH01,RE(8D) , ( BRHOl , RE( 80 ) ) , ( CRHOl , RE( 7 9 ) ) 

EQUIVALENCE ( AZW2, RE(73)),( BZH2, RE(74) ) , 



(PINT,RE(75)),( CZW2 , RE( 76 ) ) , 

(DZ,RE(78)), (ZWL,RE(77) ), (ZW2M1T, RE(7 2) ) , 
(VELMIN,RE(70)), ( VELMAX, RE( 69 ) ) , 

( FA L COM, RE( 68 ) ) , (CMDOT,RE(67 ) ), 

( RLXMDT , RE( 66 ) ) , ( AMDTMX , RE( 6 5 ) ) , ( AMDTMN, RE( 64 ) ) 
( RADMAX, RE(63)),(RADMIN,RE(62)), 
(PRESS,RE(61)),(CP1,RE(60))> 



MEL 07 21 0 
MEL 07 220 
MEL 07 230 
MEL 07 240 
MEL07250 
MEL 07 26 0 
MEL 07 27 0 
MEL 07 280 
MEL 07290 
MEL 07 30 0 
MEL 07 310 
MEL 07 320 
MEL 07 330 
MEL 07 340 
MEL 07 350 
MEL 07 36 0 
MEL 07 37 0 
MEL 07 38 0 
MEL 07 39 0 
MEL 07 400 
MEL07410 
MEL 07 420 
MEL 07 430 
MEL07440 
MEL 07 450 
MEL 07 46 0 
MEL 0747 0 
MEL 07 48 0 
MEL 07 49 0 
MEL 07 50 0 
MEL 07 510 
MEL 07 520 
MEL 07 530 
MEL 07 540 
MEL 07 550 
MEL 07 56 0 
MEL 07 57 0 
MEL 07 58 0 
MEL 07 59 0 
MEL 076 00 
MEL 07 6 1 0 
MEL 076 20 
MEL07630 
MEL 07640 
MEL076 50 
MEL 07 66 0 
MEL 07 67 0 
MEL07680 
MEL 07690 
MEL077 00 
MEL 07710 
MEL 077 20 
MEL077 30 
MEL 077 40 
MEL07750 
MEL 077 6 0 
MEL0777 0 
MEL07780 
MEL 077 90 
MEL 078 00 
MEL 0781 0 
MEL 07 820 
MEL 07830 
MEL 07840 
MEL 07 850 
MEL 0786 0 
MEL 0787 0 
MEL07880 
MEL 0789 0 
MEL 07 9 00 
MEL 07 9 1 0 
MEL 07 920 
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& CCP2,REC59)), CCP3,REC58)), C FSTOIC, REC 57 ) ), 

a CARRC0N,REC56)), C PREEXP , REC 55 ) ) , CCA1I ,REC54) ) 

EQUIVALENCE CCA2I , RE( 53 ) ) , ( CA1 3, RE( 52) ) , ( CA23, RE( 51 ) ) , 
a CH1SAT,REC50)),CH2SAT,REC49)), 
a ( EMULAM, REC48 ) ) 

EQUIVALENCE C TWPRCN , RE( 42 ) ) , ( PSATEX , RE( 41 ) ) 
EQUIVALENCECELEXP,REC39)),CEWALL,REC38)), 
a CSIGE,REC37)),CSIGK,REC36)),CTAUDK,REC35)) 

DIMENSION L0C1C6),L0C2C6),L0C3C6),L0C4C6),L0C5C6) , 
a LOC6C6),LOC7C6),LOC8C6),LOC9C6),LOC10C6) 

EQUIVALENCE C LOCREGC 1 ) , LOCI C 1 ) ) , ( L 0CREGC7 ) , L0C2C 1 ) ) 

EQUIVALENCE C lOCREG( 13) , L0C3C 1 ) ) , ( LOCREG( 19 ) , L0C4C 1 ) ) 
EQUIVALENCE C L OCREGC 25 ) , L0C5 (1 ) ) , ( LOCREGC 31 ) , L0C6 (1 ) ) 
EQUIVALENCE C LOCREGC 37 ) , LOC7 ( 1 ) ) , ( LOCREGC 43 ) , L 0C8 ( 1 ) ) 
EQUIVALENCE C LOCREGC 49 ) , LOC9 ( 1 ) ) , ( L OCREG( 55 ) , LOCI 0 ( 1 ) ) 
DIMENSION XUFRACC30) ,YVFRACC30) , ZWFRACC 30) , ENDITC 25 ) 

DIMENSION XFRACC30),YFRACC30) , ZFRACC30) 

EQUIVALENCECXUDIST Cl), XU FR ACC 1 ) , XFRACC 1 ) ) , 
a CYVDISTC1),YVFRACC1) ,YFRACC1) ), 

a C ZWDIST C 1), ZWFRACC 1),ZFRACC1)),CLITCNC,LITC), 

a CENDITC1),CRITC1)),CRLXP,RELAXP) 

EQUIVALENCE CNPRMON, LITERC 1 ) ) 

DIMENSION PHIC100,1),FDTAPEC3) 

INTEGER FTAPE 
LOGICAL FIRST 

INTEGER P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,EP,H1,H2,H3,C1,C2, 
aC3,C4 

DATA P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,KE,EP,H1,H2,H3,C1,C2, 
aC3,C4/l, 2 , 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20/ 

DATA FIRST, FTAPE, FDTAPE/ . TRUE . , 9 , 4HFI EL , 4HDS . D,2HTA/ 

IF(FIRST) CALL TAP ESC FTAPE, FDTAPE, 3 , 1 , 4*NX*NY) 

FIRST=. FALSE. 

C 

CHAPTER 2 USER CHAPTER FOR SETTING INITIAL FIELDS. 

C 

NNZ=NZ 

IFCPARAB) NNZ=1 
DO 2000 IZ=1,NNZ 
DO 2000 MPHI=1 , 25 

IFCMPHI . EQ . 2) GO TO 2000 

IFCIFIXCFIINIT (MPHI )+0.1).NE. 10101) GO TO 2000 
IFC . NOT . STOVARCMPHI ) . AND . . NOT . SOLVARCMPHI ) ) GO TO 2000 
C*XXXXX*XXXXXXUSER SECTION STARTS HEREXXXXXXXXXXXX 

FOR EACH VARIABLE, MVAR, FOR WHICH FI INITCMVAR) IS SET TO 
10101.0 IN BLOCK DATA, THE USER MUST INSERT HERE FORTRAN 
STATEMENTS FOR SETTING HIS WANTED NON-UNIFORM INITIAL FIELD. 
THE FORTRAN-STATEMENT MODEL TO BE FOLLOWED BY THE USER FOR 
EACH MVAR IS AS FOLLOWS... 

IFCMPHI .NE. MVAR) GO TO 205 
DO 201 IX=1 , NX 
DO 201 IY=1 , NY 

201 PHICIY, IX) = USER INSERTS FIELD VALUE AT IX,IY,IZ HERE. 

WRITEC FTAPE) C C PHI C IY, IX) , IY=1 , NY) , IX=1 , NX) 

205 CONTINUE 

INSERT HERE ADDITIONAL FORTRAN SEQUENCES, SIMILAR TO THAT 
GIVEN ABOVE, FOR EACH MVAR FOR WHICH FIINITCMVAR)=10101 . 0 
CxxxxxxxxxxxxxUSER SECTION ENDS HERE xxxxxxxxxxxx 
WRITEC 6 , X ) 'xxxxxx BEGINNING ’,MPHI 
IF CMPHI.NE. 5) GOTO 210 
DO 201 IY=1 , 25 

201 PHICIY, 1)=0.0 
DO 202 IY=26 , 28 

202 PHICIY, 1)=0. 0375 
DO 203 IY=29, 100 

203 PHICIY, 1)=0. 00001 

WRITEC FTAPE) C PHI CIY,1) ,IY=1,100) 

WRITEC6 ,*) ' WRITTEN *,MPHI 

GOTO 270 

210 IF CMPHI.NE. 9) GOTO 230 
DO 211 IY=1 , 100 



MEL 07 930 
MEL07940 
MEL07950 
MEL07960 
MEL07970 
MEL07980 
MEL 07 99 0 
MEL08000 
MEL 08 0 1 0 
MEL08020 
MEL08030 
MEL 08040 
MEL08050 
MEL08060 
MEL 03 07 0 
MEL03080 
MEL08090 
MEL 031 00 
MEL 08110 
MEL08120 
MEL08130 
MEL08140 
MEL 08150 
MEL 0816 0 
MEL08170 
MEL08180 
MEL08190 
MEL 08200 
MEL 08210 
MEL 08220 
MEL 08230 
MEL 08240 
MEL 08250 
MEL08260 
MEL 0827 0 
MEL 08280 
MEL 08290 
MEL 08 30 0 
MEL08310 
MEL 08 320 
MEL08330 
MEL 08340 
MEL08350 
MEL08360 
MEL0837 0 
MEL08380 
MEL08390 
MEL08400 
MEL08410 
MEL 08420 
MEL 08430 
MEL08440 
MEL08450 
MEL08460 
MEL 0847 0 
MEL 08480 
MEL 08490 
MEL08500 
MEL08510 
MEL 08520 
MEL 08530 
MEL08540 
MEL08550 
MEL08560 
MEL0857 0 
MEL08580 
MEL08590 
MEL 086 00 
MEL08610 
MEL08620 
MEL08630 
MEL08640 
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>11 PHKIY, 1) = 1 .0 

WRITEC FTAP E)C PHI CIY,1), IY=1, 100) 

WRITEC 6 , X ) 1 WRITTEN 

GOTO 270 

>30 IF C MPHI .NE. 14 ) GOTO 240 
DO 213 IY=1,25 
>13 PHKIY, 1) = 1.58E06 
PHI C 26 , 1)=3.486E06 
PHI C 27 , 1 ) =6 . 67 0 E06 
PHIC28,1)=6 . 77 0 E06 
PHI C 29 , 1 ) =6 . 3E06 
PHI C 30 , 1)=4.225E06 
PHI( 31 , 1 ) =2 . 22E06 
PHI (32, 1 ) = 2 . 37 E05 
PHI(33, 1)=-1 .72E06 
PHIC34,1)=-3.2E06 
PHIC35,1)=-6.57E06 
PHI ( 36 , 1 ) = -7 . 03E06 
DO 216 IY=37,100 
>16 PHKIY, 1)=-1.28E07 

WRITEC FTAPE)( PHI ( IY,1), IY=1 ,100) 

WRITEC 6 , x ) * WRITTEN 

GOTO 270 

!40 IF C MPHI .NE. 17 ) GOTO 250 
DO 217 IY=1 , 25 
>17 PHI C IY, 1 ) =1 . 0 
PHI C 26 , 1 ) = 0 . 9 
PHI C 27 , 1 ) = 0 . 8 
PHI C 28 , 1 ) = 0 . 7 
PHIC29,1)=0.6 
PHIC30,1)=0.5 
PHIC31, 1)=0.4 
PHIC32,1)=0.3 
PHIC33, 1 ) = 0 . 2 
PHI C 34 , 1 ) =0 . 1 
PHI C 35 , 1 ) = 1 . E“1 0 
DO 219 IY=36 , 100 
!19 PHI C IY, 1 ) =1 .E-10 

WRITEC FTAPE)( PHI C IY, 1 ), IY=1, 100) 

WRITEC 6 , X ) 1 WRITTEN 

GOTO 270 

!50 IF C MPHI .NE. 18 ) GOTO 260 
DO 251 IY= 1 , 25 
!51 PHI C IY, 1 ) = 1 . E-10 
PHIC26, 1 )=0 . 05 
PHI C 27 , 1 ) = 0 . 1 55 
PHI C 28 , 1 ) = 0 . 26 
PHIC29,1)=0 .365 
PHIC30,1)=0.47 
PHI C31,l)=0.575 
PHI C32,l)=0.68 
PHIC33,1)=0 .785 
PHIC34,1)=0.89 
PHI C 35 , 1 ) =1 .0 
DO 253 IY=36 ,100 
!53 PHI C IY, 1 ) =1 . 0 

WRITEC FTAPE)C PHI C IY, 1 ), IY= 1,1 00) 

WRITEC 6 , ^ ) 1 WRITTEN 

GOTO 270 

!6 0 I FC MPHI .NE. 19 ) GOTO 270 
DO 261 IY=1 , 25 
16 1 PHI C IY, 1 ) =1 . E-l 0 
PHIC26,1)=0.05 
PHI C 27 , 1 ) =0 . 045 
PHI C28,l)=0.04 
PHI C 29 , 1 ) =0 . 035 
PHI C30,l)=0.03 
PHI C 31 , 1 ) = 0 . 025 
PHIC32,1)=0.02 
PHIC33,1)=0.015 
PHI C 34 , 1 ) = 0 . 0 1 
PHI C 35, 1 ) = 1 . E-l 0 



MEL 086 50 
MEL 0866 0 

1 > MPHI MEL 0867 0 

MEL 08680 
MEL 0869 0 
MEL 087 0 0 
MEL 087 1 0 
MEL 087 20 
MEL 087 30 
MEL 08740 
MEL 087 5 0 
MEL 0876 0 
MEL 0877 0 
MEL 08780 
MEL 087 90 
MEL 083 00 
MEL 033 1 0 
MEL 0882 0 
MEL 0333 0 
MEL 088 40 
MEL 033 5 0 

f , MPHI MEL 0886 0 

MEL 0887 0 
MEL 08830 
MEL 0889 0 
MEL 089 0 0 
MEL 03 9 1 0 
MEL 089 20 
MEL 089 30 
MEL 03940 
MEL 03 950 
MEL 03 96 0 
MEL 0397 0 
MEL 0398 0 
MEL08990 
MEL 09000 
MEL 09 0 1 0 
MEL 09 020 
MEL 09 03 0 

1 , MPHI MEL 09 040 

MEL 09050 
MEL 09 06 0 
MEL 09 07 0 
MEL 09 08 0 
MEL 09 09 0 
MEL 09 1 0 0 
MEL 09 1 1 0 
MEL09120 
MEL 09 1 30 
MEL 091 40 
MEL 09 1 50 
MEL 09 16 0 
MEL 09 17 0 
MEL 09 18 0 
MEL 09 1 9 0 
MEL 09 20 0 
MEL 0921 0 

1 , MPHI MEL 09 220 

MEL 09 230 
MEL 09240 
MEL 09 250 
MEL 0926 0 
MEL 0927 0 
MEL 0928 0 
MEL09290 
MEL 093 0 0 
MEL0931 0 
MEL 09320 
MEL 09330 
MEL 09 340 
MEL 09 350 
MEL 0936 0 
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262 

270 

2000 



DO 262 IY=36 , 1 00 
PHI ( IY, 1 ) =1 . E-l 0 
WRITE(FTAPE)(PHI(IY,1 
WRITEC6 , X) 

CONTINUE 

CONTINUE 

RETURN 

END 



)»IY=1,100) 

» WRITTEN 



' > MPHI 



MEL 0937 0 
MEL0938 0 
MEL09390 
MEL 09400 
MEL 09410 
MEL09420 
MEL09430 
MEL09440 
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